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Inertial  Instability  of  Nonparallel  Flow  on  an  Equatorial  0  Plane 

Timothy  J.  Dunkerton 
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(Manuscript  received  22  June  1992.  in  final  form  7  December  1992) 

ABSTRACT 

A  simple  theoretical  model  was  developed  to  investigate  the  inertial  instability  of  zonallv  nonuniform,  non¬ 
parallel  flow  near  the  equator.  The  basic  state  was  independent  of  height  and  time  but  included  cross-equatorial 
shear  with  longitudinal  variation,  as  observed  in  the  tropical  mesosphere  and  elsewhere.  Numerical  solutions 
were  obtained  for  the  most  unstable  modes. 

It  is  shown  that,  in  addition  to  previously  known  “global”  (symmetric  and  nonsymmetric)  modes  of  inertial 
instability,  there  exist  “local”  modes  within  regions  of  anomalous  potential  vorticity.  Local  modes  may  be 
exactly  stationary  or  display  zonal  phase  propagation,  but  are  distinguished  from  global  modes  by  their  zero 
group  velocity  and  concentration  of  amplitude  within,  or  downstream  from,  the  region  of  most  unstable  flow. 
Local  stationary  instability  has  the  largest  growth  rate  and  occurs  in  strong  inhomogeneous  shear  when  the  in 
situ  mean  flow  is  near  zero,  that  is,  quasi-stationary  with  respect  to  the  (stationary)  basic-state  pattern.  This 
situation  is  expected  in  an  equatorial  Rossby  wave  critical  layer. 

The  local  mode  has  properties  similar  to  those  of  “absolute”  instability  of  nonparallel  flow  as  discussed 
elsewhere  in  fluid  dynamics. 


1.  Introduction 

Inertial  instability  may  arise  in  conservative  axisym- 
metric  flow  near  the  equator  when  there  is  nonzero 
latitudinal  shear  y  =  uy.  In  this  case  potential  vorticity 
(PV)  on  one  side  of  the  equator  is  anomalous,  satisfying 
a  necessary  condition  for  centrifugal  parcel  instability 
(Dunkerton  1981 ).  Analytic  eigenmodes  of  symmetric 
instability  in  zonallv  uniform  shear  on  an  equatorial 
0  plane  were  obtained  by  Dunkerton  and  indepen¬ 
dently  by  Stevens  (1983).  Subsequently  Boyd  and 
Christidis  ( 1982)  and  Dunkerton  (  1983)  found  that 
low- wavenumber,  zonallv  nonsymmetric  instabilities 
have  larger  growth  rates  than  symmetric  instability 
when  (i)  the  vertical  wavenumber  m  is  near  or  below 
a  neutral  point  of  symmetric  instability 


(N  is  static  stability  and  0  =  212 /a  =  2.29 
X  10"11  m_l  s~!),  and  (ii)  nonzero  integer  zonal 
wavenumbers  s  *  ka  exist  below  a  short-wave  cutoff 
km**  -  0/y-  Although  symmetric  instability  has  max¬ 
imum  inviscid  growth  approaching  y/2  as  |  m  |  oc , 
the  nonsymmetric  mode  is  preferred  when  second-or¬ 
der  viscosity  and  diffusion  are  added  and  conditions 
(i)  and  (ii)  are  met  (Dunkerton  1983).  Similar  con¬ 
clusions  pertain  to  midlatitude  flow  containing  non- 
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symmetric  inertial  and  barotropic  instabilities,  as 
shown  by  Stevens  and  Ciesielski  ( 1986). 

Anomalous  PV  exists  in  the  tropical  winter  meso¬ 
sphere  (Dunkerton  1981;  Hitchman  andLeovv  1986), 
upper  troposphere  between  South  America  and  the  In¬ 
dian  Ocean  ( Liebmann  1 987 ) ,  lower  troposphere  near 
Indonesia  ( Krishnamurti  et  al.  1985.  1988),  and  Pacific 
south  equatorial  current  (Philander  1989,  p.  62).  In¬ 
ertial  instabilities  were  found  in  a  nearly  inviscid  2D 
axisymmetric  model  of  the  troposphere  ( Held  and  Hou 
1980)  and  in  3D  middle  atmosphere  GCMs  (Hunt 
1981;  O’Sullivan  and  Hitchman  1992).  This  instability 
was  suggested  as  a  possible  explanation  of  layered 
structures  in  the  tropical  mesosphere  ( Hitchman  et  al. 
1987;  Fritts  et  al.  1992)  and  mesoscale  anomalies  in 
the  midlatitude  upper  troposphere  (Ciesielski  et  al. 
1989).  Conditional  symmetric  instability  is  thought  to 
explain  frontal  rainbands  according  to  Bennetts  and 
Hoskins  ( 1979)  and  several  others. 

Unfortunately  for  the  theory,  none  of  the  regions  of 
anomalous  PV  in  atmosphere  or  ocean  is  zonallv  sym¬ 
metric.  Observations  and  numerical  models  suggest, 
for  example,  that  inertial  instabilities  in  the  mesosphere 
coincide  with  tropical  penetration  of  planetary  Rossby 
waves  from  midlatitudes  (Hitchman  et  al.  1987; 
O’Sullivan  and  Hitchman  1992).  For  the  theory  to  be 
relevant,  it  must  be  generalized  to  a  nonparallel  or 
zonally  nonuniform  basic  state. 

The  purpose  of  this  paper  is  to  develop  a  simple 
theoretical  formalism  and  thereby  to  document  two 
types  of  inertial  instability  in  nonparallel  flow:  local 
instability  (stationary  and  zonally  propagating)  and 


^  1993  American  Meteorological  Society 


2 


15  AUGUST  1993 


DUNK.ERTON 


2745 


global  instability  (symmetric  and  nonsymmetric).  Of 
these,  the  local  instabilities  are  new  and  probably  most 
relevant  to  middle  atmosphere  models  in  which  break¬ 
ing  planetary  waves  produce  regions  of  significantly 
anomalous  PV  near  the  equator. 

The  theoretical  model  is  described  in  section  2,  in¬ 
cluding  a  brief  review  of  symmetric  and  nonsymmetric 
instability.  Growth  rates  and  structure  of  nonparallel 
instabilities  are  illustrated  in  section  3  for  a  wide  range 
of  model  parameters  using  a  simple  analytic  basic  state 
resembling  the  onset  of  planetary  wave  breaking.  The 
theoretical  interpretation  of  local  instability  is  discussed 
further  in  section  4  and  the  Appendix. 

2.  Background 

To  isolate  the  effect  of  horizontal  basic-state  varia¬ 
tions,  a  simple  two-dimensional  ( x-y )  model  is  desir¬ 
able,  analogous  to  the  shallow-water  system.  In  a  ver¬ 
tically  continuous  atmosphere,  these  equations  describe 
perturbations  oscillating  in  height  with  constant  si¬ 
nusoidal  variation.  Vertical  wavelength  is  then  con¬ 
tained  in  Lamb’s  parameter  «  (Andrews  et  ai.  1987) 
and  is  specified  a  priori.  In  this  section,  the  simple 
model  is  developed,  its  solution  method  described,  and 
a  brief  review  given  of  parallel  flow  instabilities  (several 
features  of  which  carry  over  to  the  nonparallel  case). 

a.  Theoretical  model 

The  hydrostatic  primitive  equations  on  an  equatorial 
0  plane  are 

u,  +  uux  +  v(uy  -  0y)  +  wu:  +  =  X  (2. la) 

v,  +  u(vx  +  0y)  +  vvy  +  wv.  +  4>y=Y  (2.1b) 

<\>zt  +  +  v</>;y  +  w(N2  +  =  Q  (2.1c) 

ux  +  vy -i - (pow)r  =  0  (2.1d) 

Po 

where  u ,  v,  w  are  zonal,  meridional,  and  vertical  ve¬ 
locity,  <j>  is  geopotential,  p0  is  basic-state  density  ps 
exp(  —z/H),  and  X,  Y,  Q  are  dissipative  terms.  Equa¬ 
tions  (2.1 )  are  essentially  those  of  Andrews  and  Mc¬ 
Intyre  ( 1976)  written  in  log-pressure  notation  of  Hol¬ 
ton  (1975). 

To  formulate  an  idealized  model  of  instability  the 
dependent  variables  may  be  expanded  as 

m(x,  y,  z,  t )  =  U(x,  y)  +  u’{x,  y,  z,  t)  (2.2) 

and  similarly  for  v,  w,  <t>.  Here  U(x,  y)  is  a  barotropic 
basic  state  independent  of  height  and  time  for  the  pur¬ 
pose  of  linear  stability  analysis.  Linearized  perturbation 
equations  are  then 

«;  +  uu'x  +  v'(  uy  -  / 3y)  +  <p'x  =  A  +  KVyu'  -  wwV 

(2.3a) 


v'i  +  iiv’y  +  Byir  +  < j>’y  =  B  +  KTyv'  -  vmr v'  ( 2.3b ) 
<t>,  +  ii<px  +  t  ( ux  +  u'y )  =  C  +  KTjf(j)'  —  vmz<j> 

(2.3c) 

under  the  foilowing  assumptions: 

1 )  The  basic-state  velocity  field  is  barotropic.  tem¬ 
porally  constant,  and  horizontal  ( IV m  0).  More  gen¬ 
erally  we  could  regard  the  mean  flow  as  slowly  varying 
in  height  and  time  ( Boyd  1978). 

2)  Advection  of  basic-state  temperature  is  neglected, 

that  is,  u'  •  =  o.  This  term  has  a  counterpart 

u'  •  V„4>  normally  retained  in  the  shallow-water  sys¬ 
tem.1 

3 )  Perturbations  vary  as  exp imz,  and  nr  >  1  /  (4 H2) 
whereupon  s  =  m2/N 2  as  in  ( 1.1). 

4)  \<t>::\  <N3. 

5)  Dissipative  terms  are  written  as  second-order 
viscosity  and  diffusion. 

In  (2.3a-c)  the  zonally  symmetric  component  of  U 
was  kept  on  the  lhs,  that  is. 


U-»u(y)  +  U(x,  y) 

(2.4a) 

V(x,y) 

(2.4b) 

on  the  rhs 

-A  =  ( Uux  +  Vu’y)  +  (u'Ux  +  V'Uy) 

(2.5a) 

-B  =  ( Uv'x  +  Vv'y)  +  (u'Vx  +  V'Vy) 

(2.5b) 

~C  a  (  U<t>X  +  V  ()>'y). 

(2.5c) 

b.  Solution  method 


To  solve  (2.3a-c)  a  semi-implicit  method  was  used 
in  which  variables  were  expanded  in  zonal  harmonics 
(letting  k  be  integers  for  the  moment). 


f 

]  iok 

u*. 

uk 

Vk 

dk 
•  Bk 

ia 


exp 


ikx 


exp 


exp  • 


a 

ikx 

a 

ikx 


a 


(2.6a) 

(2.6b) 

(2.6c) 


Terms  on  the  lhs  of  ( 2.3a-c)  were  treated  implicitly  in 
time,  that  is. 


1  In  real  atmospheres  the  basic  state  will  include  variations  of  static 
stability  induced  by  potential  vorticity  anomalies  (Hoskins  et  al. 
1985).  For  simplicity,  only  the  variation  of  induced  horizontal  velocity 
U.  V  was  retained  in  the  definition  of  A,  B,  C.  This  assumption  is 
reasonable  in  the  tropical  middle  atmosphere,  where  incident  Rossbv 
waves  have  approximately  horizontal  ray  paths. 
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Cik  *  \(unk  +  «r‘),  etc.,  (2.7) 

so  that  for  each  harmonic 
k 

+  Sjt(wr  ~  0y)  4-  -  <£*  =  -tig  -  will  (2.3a) 

'•  +uvk  +  $yuk  +  —  0*  =  Bl  +  uvnk  (2.8b) 

.  -  .  (  k  „  d  \ 

-u<j>k  +  e  (^-  w*  +  I  =  -iCnk  -  w0Z,  (2.8c) 

where  co  m  u  -  ku/ a,  co  »  2i/  At,  and 

Ac  -  At  +  diffusion  terms,  etc.  (2.9) 

After  some  manipulation  (2.8a-c)  reduce  to  a  geo¬ 
potential  equation 


Zeuiky 

0\v  +  — 7T 

ew  -  k~  • 


-  u 


w  ecj-  —  k~ 


=  0 


(2.14) 

which  is  nonlinear  in  e  unless,  for  example,  ear  >  Ar 
in  which  case 


2ky 

vyv  +  — r-  vy  -  v 
oj 


to 


(0  -  Ty)  +  eA 


0.  (2.15) 


If  k  =  0 

vyy  -  ev[fiy((3y  -  y)  -  or]  =  0  (2.16) 

and  the  eigencondition  for  symmetric  instability  is 


A 


d_  1 0 y\ 

dy  \  A  / 


—  eA<£  4-  2) 


or  +  72/4  = 


(2/2  4~  QJVff 

|w| 


(2.17) 


(2.10)  (Dunkerton  1981 ).  Solutions  for  v  are  Hermite  funo 
for  each  harmonic,  where  tions  (polynomials  times  a  Gaussian)  centered  about 

*  a  shifted  equator”  ys  m  y/2f3 ,  the  center  of  anomalous 

d  53  py(py  ~  y)  -  (2. 1 1 )  vorticity  region.  For  n  =  0  the  neutral  point  of  sym- 

In  (2.10)  X),  whose  definition  is  omitted  for  brevity,  metric  mstabihty  (hereafter,  “neutral  point”)  is  given 

includes  the  explicit  forcing  and  solutions  from  the  by  ^ 1-1 ''  For  sma  *  and  Ty  ~ 
previous  time  step.  Equation  (2. 10)  was  solved  by  dis-  ,  N0  kB 

crctizing  in  y  with  64  grid  points  and  using  a  tridiagonal  w  +  7  /4  s»  - — -  +  —  .  (2.18a) 

algorithm.  The  domain  was  periodic  in  x  and  extended  m  we 


from  -4ys  to  +6ys  where  ys  a  7/2/3.  There  are  no 
critical  latitude  («  =  0)  or  inertial  latitude  (A  =  0) 
singularities  when  a>  is  complex.  A  spectral  rather  than 
pseudospectral  method  was  used  to  evaluate  D  since 
in  this  case  Uk,  Vk  =  0  for  any  |  k\  ¥=  1;  that  is,  a  purely 
wave  1  contribution  to  basic-state  variation  in  x  was 
assumed.  Product  terms  Ak,  etc.,  could  be  determined 
efficiently  without  transforms.  In  the  numerical  inte¬ 
grations 

ii  =  y{y  -  yc)  +  Uq  (2.12) 

where  y  and  z70  are  constants.  Fields  of  U,  V are  spec¬ 
ified  in  section  3. 

c.  Instability  of  uniform  parallel  flow 

When  U  =  V  =  0,  harmonic  perturbations  propor¬ 
tional  to  exp/(Ar.v  —  (^t)  satisfy 

-ww  +  t>(7  -fly)  +  k<t>  =  0  (2.13a) 

+wv  +  0yu  +  4>y  =  0  (2.13b) 

-we0  +  ku  +  vy  =  0  (2.13c) 


^  eneut  , 

3 

—  .  (2.18b) 

Equations  (2.18a,b)  were  derived  in  nondimensional 
form  by  Boyd  and  Christidis  (1982). 

At  finite  k  the  eigenproblem  (2.10)  or  (2.14)  must 
be  solved  numerically.  Growth  rates  for  y  =  10-3  s-1 
are  shown  in  Fig.  1  illustrating  the  relative  importance 
of  symmetric  and  nonsymmetric  modes.  Curves  were 
obtained  from  the  geopotential  equation  by  a  shooting 
method;  symbols  correspond  to  values  from  the  time- 
dependent  model.  They  are  in  excellent  agreement  ex¬ 
cept  at  very  small  growth  rates,  where  the  time-depen¬ 
dent  model  could  not  determine  growth  rate  accurately, 
and  at  large  e,  where  latitudinal  resolution  was  inade¬ 
quate. 

An  approximate  model  derived  from  (2.15)  contains 
eigenfrequencies  similar  to  the  exact  values  of  Fig.  1 
over  a  large  portion  of  parameter  space.  This  model 
begins  with  constant  y  and  assumes  that  y  variations 
of  intrinsic  frequency  can  be  neglected  outside  A: 


where  co  =  o>  —  ku.  The  geopotential  equation  derived  2ky 

from  (2. 13a— c)  is  just  (2.10)  with  !D  —  0  and  o>  rede-  ~  —  ® 

fined  as  intrinsic  frequency.  This  equation  is  linear  in 
the  eigenvalue  t.  Boyd  (1978)  derived  the  meridional  where  co0  is  a  constant  (complex)  intrinsic  frequencv 
velocity  equation  within  the  unstable  region.  It  follows  that 


kff 

—  +  (A 

to0 


0  (2.19) 
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Fig.  1.  Growth  rate  as  a  function  of  zonal  wavenumber  in  parallel 
flow,  for  t  -  10”5  s’*;  €  in  units  m**2  s2.  Symbols  denote  results  of 
time-dependent  numerical  model;  solid  and  dashed  curves  show  an¬ 
alytic  results  from  a  shooting  method.  (The  dashed  curve  had  no 
symmetric  instability  because  e  < 


(o>  -/cy2/#)2  1  , 

1  -k2y2/02  47’ 


N0 


k& 


Vl-*V/0a  +  - 

m  |  o)0e 


(2.20) 


where 


flow.  Some  characteristics  of  these  modes  are  now  de¬ 
scribed.  The  basic  state,  independent  of  height  and 
time,  was  designed  as  a  simple  analog  of  incipient 
planetary  wave  breaking  in  the  tropical  mesosphere. 
The  evolution  of  inertial  instability  in  a  time-dependent 
Rossby  wave  critical  layer  was  recently  discussed  bv 
O’Sullivan  and  Hitchman  ( 1992)  in  a  three-dimen¬ 
sional  model.  In  order  to  understand  the  selection  of 
unstable  modes,  it  will  be  worthwhile  to  explore  the 
parameter  dependence  of  inertial  instability  in  a  sim¬ 
pler  two-dimensional  model. 

a.  Basic  state 

The  basic  state  was  defined  by 
ii(y)  +  U(x,y)  =  y(y-yc) 

cos (x/a)  +  Uo.  (3.1) 

In  (3.1),  the  latitude  of  maximum  shear  is  yc,  \p  is 
dimensionless  amplitude  of  basic-state  variation,  ys 
=  7/2/8  is  the  center  of  anomalous  zonal-mean  vortic- 
ity  as  defined  previously,  and  ilo  is  a  constant  mean 
flow  independent  of.v,  v.  The  profile  (3.1)  is  shown 
in  Fig.  2  for  yc  =  1.8  vy,  \p  =  1,  and  Uo  =  0. 

In  Fig.  2,  latitudinal  shear  at  y  =  yc  was  increased 
(decreased)  relative  to  the  zonal  mean  by  50%  at  180° 
(0°)  longitude.  The  flow  was  inertially  unstable  at  all 
x,  but  more  so  at  the  center  of  figure.  This  region  could 
represent  a  zone  of  enhanced  latitudinal  shear  set  up 
by  Rossby  waves  penetrating  from  latitudes  north  of 
the  equator.  In  reality,  there  would  be  some  latitudinal 
phase  tilt  during  critical  layer  development  (O’Sullivan 
and  Hitchman  1992) — an  unnecessary  complication 
for  stability  analysis. 

The  meridional  component  of  basic  state  V(x,  y) 
was  included  in  all  simulations,  that  is. 


rPy* 


+  tanh(^ 

\  ys 


coo  —  w  —  ky2 


0-2uk 
2(fi2  —  k2yz)  ■ 


(2.21) 


A  formal  justification  of  approximations  leading  to 
(2.20)  will  not  be  given  except  to  note  that  terms  of 
0(k2y2)  arise  from  A  at  large  y  and  cannot  be  ignored 
for  latitudinal  trapping  and  quantization  of  eigenfunc¬ 
tions  near  the  short-wave  cutoff.  The  approximate  dis¬ 
persion  relation  ( 2.20 )  has  the  desired  behavior  at  small 
and  large  k.  In  between,  agreement  with  exact  growth 
rates  and  phase  speeds  is  good  for  e »  eneut  ( not  shown ). 
At  smaller  e  the  cubic  equation  fails  in  two  ways:  growth 
rates  are  overestimated  at  intermediate  k  and  there  is 
no  Kelvin  wave  instability  at  small  k  below  the  neutral 
point  (Boyd  and  Christidis  1982). 


3.  Instability  of  nonparailel  flow 

The  time-dependent  model  of  section  2  was  used  to 
determine  the  growth  of  unstable  modes  in  nonparallel 


U,  +  Vy  =  0.  (3.2) 

This  had  a  minor  effect  compared  to  the  zonal  com¬ 
ponent  (3.1).  Although  the  details  are  omitted,  it  can 
be  stated  that  all  of  the  results  shown  here  were  similar 
when  V  =  0.  Local  instabilities  were  also  found  when 
U  3  0  and  V  -  V(x)  (localization  due  to  dV/dx). 
The  effect  of  V  was  generally  to  disrupt  rather  than 
enhance  the  instability  in  comparison  to  parallel  flow. 
These  cases  were  felt  to  be  less  interesting  and  were 
not  investigated  further. 

b.  Latitude  of  maximum  shear 

The  latitude  of  maximum  latitudinal  shear  was  yc 
according  to  ( 3. 1 ).  On  the  other  hand,  the  maximum 
rate  of  “parcel  instability"  isi?  =  y—(Jy(f}y  +  f)  where 
f  is  relative  vorticity.  This  quantity  is  contoured  in  Fig. 
2:  it  did  not  always  maximize  at  y  ~  yc  but  was  a  func¬ 
tion  of  (and  other  parameters),  having  (for  \p  =  I ) 
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Fig.  2.  Profile  of  nonparallei  basic  state  (3.1)  with  =  1.  Contours  show  dimensionless  parcel 
displacement  growth  rate  R/y.  Contour  interval  0.1,  with  0.1  and  0.4  contours  omitted. 


an  overall  maximum  value  R  =  0.75y  when  y  =  yc 
=  1 .5ys.  Maximum  R  was  located  north  of  ys  because 
the  wave’s  contribution  to  basic  state  expanded  the 
zone  of  instability  at  180°  (recall  that  ys  was  the  center 
of  anomalous  zonal-mean  vorticity).  As  \ye\  -*■  oo, 
however,  the  relative  vorticity  at  |  yc  |  became  small 
relative  to  j3y,  so  the  region  of  anomalous  vorticity 
next  to  the  equator  returned  to  a  parallel  configuration. 
For  parallel  flow  the  maximum  R  =  y/2  ax  y  -  ys. 

Growth  rates  of  the  most  unstable  mode  as  a  func¬ 
tion  of  yc  are  shown  in  Fig.  3  for  y  =  4  X  10“5  s“!, 
=  256  nr  s~2,  =  l,  and  iio  -  0.  For  these  cal¬ 

culations  eight  zonal  harmonics  were  used,  and  It 
-  900  s.  Vertical  wavelength  km  for  this  choice  of 
e  when  N  =  0.02  s~! .  The  asymptotic  forms  of  mode 
and  growth  rate  were  realized  after  a  few  days  in  most 
cases,  although  as  will  be  seen  later,  this  time  scale 
depended  on  y.  The  results  appear  reasonable:  growth 
rate  peaked  at  yc  «s  1.7  ys  and  for  large  \yc\  instabilities 
developed  as  in  parallel  flow  (not  shown).  The  inter¬ 
pretation  of  Fig.  3  is  slightly  more  complicated,  how¬ 
ever,  because:  1 )  instabilities  near  the  center  of  the  plot 
were  exactly  stationary  “s”  and  grew  about  twice  as 
fast  as  in  parallel  flow — a  larger  variation  than  expected 
from  R  alone.  Peak  growth  occurred  a  little  to  the  right 
of  1.5 ys.  2)  Instabilities  adjacent  to  those  labeled  “s” 
were  zonallv  propagating  but  not  global  in  extent. 

The  basic-state  zonal  flow  [u{y)  4*  U(x,  y)  in  (3.1 )] 
at  location  of  maximum  R(  180°,  v)  also  varied  with 
yc;  it  was  zero  at  this  point  when  yc  ^  2.2  v*.  This 
variation  affected  the  growth  rate  in  addition  to  the 


variation  of  R .  Growth  was  optimized  for  a  combi¬ 
nation  that  maximized  R  and  minimized  the  in  situ 
mean  flow.  The  dependence  of  growth  rate  on  mean 
flow  speed,  or  breakdown  of  Galilean  invariance,  is 
characteristic  of  local  instability  of  nonparallei  flow  as 
discussed,  for  example,  by  Pierrehumbert  (1984). 
We  return  to  this  point  in  the  next  subsection  and  in 
section  4. 

Instability  structure  for  yc  =  l.Sys  is  shown  in  Fig. 
4.  The  disturbance  was  exactly  stationary  and  confined 


0  12  3 


yc  /  ys 

Fig.  3.  Growth  rate  as  a  function  of yc.  Local 
stationary  instability  denoted  by  **s/‘ 
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Fig.  4.  Structure  of  local  stationary  mode  for  y  =  4  X  10'3  s'1,  e'1  =  256  m2  s'2, 
4  ~  1,  «o  =  0,  and  yc  =  1.8y,.  Negative  o  contours  are  dashed. 


to  the  middle  half  of  the  domain.  Velocity  was  directed 
from  low  to  high  pressure  indicating  a  dynamical  source 
of  instability  (Coriolis  force).  Curiously  the  geopoten¬ 
tial  phase  tilted  slightly  although  no  tilt  was  imposed 
on  the  basic  state  (3.1).  Identical  results  were  obtained 
with  doubled  resolution  in  x,  y  and  At  =  450  s. 

Figure  4  exemplifies  “local  stationary”  instability. 
Note  that  1 )  although  localized  in  x  it  is  not  infinites¬ 
imal  in  size.  This  indicates  that  terms  in  addition  to 
v'Uy  are  important  in  the  dynamics,  as  expected  from 
section  2c.  An  oversimplified  model  could  be  formu¬ 
lated  in  which  v'Uy  was  the  only  term  in  A,  all  linear 
terms  x  d/dx  were  dropped,  and  B,  C  a  0.  Growth 
rate  in  this  case  could  be  obtained  from  the  symmetric 
formula  (2.18a  with  k  =  0)  but  was  a  function  of  x, 
implying  a  5-function  catastrophe.  (This  behavior  was 
verified  numerically.)  Further,  2)  absence  of  zonal 
propagation  over  a  finite  range  of  parameters  (in  this 
case,  yc )  apparently  requires  the  zonally  symmetric  (s 
=  0)  component  of  perturbation,  an  important  part  of 
the  total  wave  field.  Most  of  the  perturbation  energy 
was  contained  in  the  lowest  four  harmonics.  When  the 
s  =  0  component  was  artificially  excluded  from  the 
numerical  model,  real  phase  speed  was  nonzero  except 
where  its  trajectory  crossed  the  cr  =  0  axis  at  a  point 
near  yc  =  1.8.1'*. 

In  parallel  flow,  nonsymmetric  instability  propagates 
zonally  except  in  the  special  case  cr(k,  <)  »  -i?0 
(Dunkerton  1983;  Stevens  and  Ciesielski  1986).  Only 
symmetric  instability  is  trivially  “stationary”  regardless 


of  parameters.  In  nonparallel  flow,  local  stationary  in¬ 
stability  is  possible.  It  is  in  a  sense  locally  symmetric, 
although  this  description  ignores  the  spatial  variation 
of  instability  structure  (e.g..  Fig.  4). 

c.  Dependence  on  u0 

The  term  z70  in  (3.1)  added  a  constant  mean  flow 
independent  of  x,  y.  In  parallel  flow  do  would  have  no 
effect  other  than  to  shift  phase  speed  by  this  constant. 
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Fig.  5.  Growth  rate  as  a  function  of  w0. 
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For  local  modes  in  nonparallel  flow,  growth  rate  is  al¬ 
tered. 

Figure  5  shows  «,*/  y  for  the  same  parameters  as  Fig. 
2  but  with  yc  =  1.8  ys  and  i70  varied  over  a  large  range. 
Local  stationary  instability  occurred  near  the  point  of 
zero  mean  flow,  but  for  large  |  ii0\,  the  instability  was 
displaced  off  center  and  forced  to  propagate  in  the  di¬ 
rection  of  M0. 

An  example  of  this  zonallv  propagating  or  advected 
instability  is  shown  in  Fig.  6  for  u0  =  -30  ms-1.  Prop¬ 
agation  was  to  the  west,  and  the  pattern  was  observed 
to  repeat  every  six  days.  Obviously  this  was  insufficient 
time  for  global  traverse,  and  in  any  case,  the  disturbance 
reentering  from  the  right  had  little  role  to  play  in  the 
sequence. 

It  is  thought  that  local  instabilities  generally  belong 
to  a  class  defined  by  “absolute”  instability  (Pierre- 
humbert  1984).  Absolute  instability  is  sensitive  to  a 


mean  flow  that  displaces  the  center  of  mode  away  from 
the  source  of  instability  and  thereby  lowers  the  effective 
growth  rate.  When  the  mean  flow  is  stable  at  large  |  x\ , 
a  sufficiently  large  1 170 1  can  stabilize  the  problem.  This 
did  not  happen  at  the  edges  of  Fig.  5  because  of  the 
zonallv  symmetric  component  of  shear  and  periodic 
boundary  conditions;  advected  modes  simply  wrapped 
around  and  continued  to  exist  as  global  modes  in  x. 
with  smaller  growth  rate.  A  sponge  region  near  x  =  0, 
to  simulate  an  infinite  domain,  dramatically  altered 
the  advected  mode  as  discussed  at  the  end  of  this  sec¬ 
tion.  The  interpretation  of  local  modes  as  absolute  in¬ 
stability  is  deferred  to  the  Appendix. 

Of  the  two  subtypes  of  local  instability  (stationary 
and  advected),  the  stationary  variety  seems  more  rel¬ 
evant  to  equatorial  Rossby  wave  critical  layers,  for  in¬ 
stance,  when  a  quasi-stationary  Rossby  wave  induces 
a  region  of  significantly  anomalous  potential  vorticity 


LONGITUDE 

Fig.  6.  Structure  of  unstable  advected  mode  for  parameters  of  Fig.  4  but  with  Uq  =  -30  m  s"1 . 
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near  the  zero-wind  line.  Here  the  advecting  mean  flow 
is  small.  That  is  presumably  why  instabilities  simulated 
by  O’Sullivan  and  Hitchman  (1992)  were  quasi-sta- 
tionary  (locally  symmetric)  and  followed  the  anoma¬ 
lous  PV  during  critical  layer  development. 

d.  Transition  to  longitudinally  global  instability 

The  value  of  yp  in  ( 3. 1 )  determined  the  zonal  vari¬ 
ation  of  shear.  Not  surprisingly,  it  was  found  that  as  yp 
decreased  there  was  a  transition  to  global  instability, 
as  shown  in  Fig.  7  ( p  signifying  a  globally  propagating 
or  parallel-flow  instability).  Growth  rate  was  reduced 
as [yp  -*  0.  Parameters  were  the  same  as  before,  but 
with  yc  -  1.8y„  ii0  =  0,  and  variable  yp.  The  structure 
became  essentially  that  of  nonsymmetric  instabilitv  in 
parallel  flow  (Stevens  and  Ciesielski  1986),  with  en¬ 
hanced  amplitude  near  the  most  unstable  region  (see 
Fig.  8,  for  yp  -  0.2).  These  were  global  modes  that 
recycled  through  the  entire  domain.  (Following 
Pierrehumbert,  “global”  refers  only  to  the  longitudinal 
extent  of  eigenmodes.) 

Results  at  very  small  p  depended  on  initial  condi¬ 
tions.  A  wavenumber  1  initial  disturbance,  used  in  most 
simulations  here,  produced  wave  1  nonsymmetric  in¬ 
stability.  Initial  wave  2  led  to  wave  2  instability  at  small 
yp,  with  slightly  larger  growth  rate.  Wave  3  was  near 
the  short-wave  cutoff  and  less  unstable. 

The  time-dependent  model  simulated  only  the  most 
unstable  mode;  however,  we  expect  local  and  global 
modes  of  instability  to  coexist,  at  least  in  cases  such  as 
Fig.  2  where  the  flow  is  unstable  at  all  x.  Whether  global 
modes  of  inertial  instability  are  important  in  the  middle 
atmosphere  or  anywhere  else  is  uncertain.2 

e.  Variation  of  e 

The  value  of  e  1  in  Figs.  3—8  (256  m2  s-2)  was  near 
the  neutral  point  =  305  m2  s“2  for  y  =  4 
X  10  3  s  1 .  Vertical  wavelength  was  ~5  km  when  N 
-  0-02  s  1 .  [The  reader  may  refer  to  ( 1 . 1 )  for  other 
N ,  and  recall  the  transformation  of  variables  in  Dunk- 
erton  (1983)  allowing  generalization  to  arbitrary  values 
of  shear.] 

Recalling  the  discussion  of  section  2c,  the  neutral 
point  limits  symmetric  instability  in  parallel  flow.  There 
is  no  neutral  curve  at  small  t  ( except  at  small  k  adjacent 
to  the  neutral  point;  see  Boyd  and  Christidis  1982)  to 
prevent  instability  at  nonzero  s  =  ka.  For  y  =  4 
X  10  ~5  s-1 ,  several  unstable  integer  wavenumbers  are 


There  is  an  intriguing  visible  light  photograph  of  Jupiter  taken 
by  the  Hubble  Space  Telescope  on  28  May  1991  revealing  horizontally 
tilted,  banded  structures  next  to  the  equator  (e.g..  see  the  cover  page 
of  AGU  publication  Earth  in  Space.  January  1992  issue).  This  would 
suggest  horizontally  divergent  circulations,  perhaps  due  to  a  dvnamical 
instability. 


allowed  below  the  short-wave  cutoff  (s  =  0-3 ).  Growth 
rate  varies  monotonically,  increasing  with  e  to  a  limiting 
value  y/2  (cf.  Fig.  1 ). 

In  nonparallel  flow,  the  growth  of  local  instability 
depended  on  e  in  a  similar  wav  as  shown  in  Fig.  9.  The 
same  paraineters  were  used  in  (3.1)  with  yc  =  i.8.y„ 
P  -  1 .  and  u0  =  0.  Growth  rate  varied  monotonically, 
reaching  a  plateau  at  higher  «.  Over  most  of  this  range 
the  instability  was  stationary,  displaying  slow  eastward 
propagation  only  at  small  e.  It  seems  clear  from  this 
example  that  nonparallel  instability  exists  below  the 
neutral  point  whether  one  adopts  a  zonal-mean  flow 
value  («ne'ut  ~  305  m2  s-2)  or  hypothetical  local  value 
(  €  neut  ~  1545  m2  s  2)  at  center  of  domain. 

The  structure  of  instability  at  small  e  had  several 
interesting  features  (not  shown):  the  latitudinal  scale 
was  enlarged,  as  in  classical  equatorial  wave  theory,  so 
instability  extended  into  the  stable  part  of  the  domain. 
Its  structure  resembled  a  Rossby  wave  in  the  Northern 
Hemisphere  and  eastward  inertia-gravity  wave  in  the 
Southern  Hemisphere  (Matsuno  1966).  The  real  phase 
speed,  though  slow  eastward,  was  for  practical  purposes 
stationary  with  respect  to  the  ground.  The  disturbance 
propagated,  with  respect  to  the  fluid,  westward  north 
of  the  equator  and  eastward  south  of  the  equator.  Be¬ 
cause  of  latitudinal  shear,  a  fused  Rossby/ inertia- 
gravity  structure  could  be  excited  by  equatorial  insta¬ 
bility  that  acted  as  a  kind  of  “wavemaker”  in  this  in¬ 
stance.  [Frontal  excitation  of  inertia-gravity  waves  in 
a  localized  region  of  conditional  symmetric  instability 
was  noted  by  Thorpe  and  Rotunno  (1989)  and  Jones 
and  Thorpe  ( 1 992 ) .]  The  instability  was,  however,  local 
in  x. 

The  structure  of  instability  at  large  e  resembled  Fig. 

4  but  was  contracted  in  latitude,  remaining  local  in  x. 

Figure  9  suggests  that  as  in  parallel  flow,  scale  selec¬ 
tion  at  finite  \m\  is  unrealized  in  the  simple  model. 


7.  Growth  rare  ns  a  function  of  4*  Globally  propagating  or 
parallel-flow  instability  denoted  by  "p." 
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Fig.  8.  Structure  of  global  instability  for  parameters  of  lug.  4  but  with  =  0.2. 


Dunkerton  (1981,  1983 )  argued  that  vertical  diffusion  matters  inasmuch  as  the  basic  state  conserves  potential 
selects  finite  \m\.  With  diffusion  there  is  a  net  growth  vorticity  on  fluid  parcels  and  inertial  instability  is  to 
rate  w/ff  =  &>,-  —  veN1.  In  the  real  world  it  is  uncertain  some  extent  a  parcel  instability  dependent  on  anom- 
whether  diffusion  causes  scale  selection  or  is  a  posteriori  alous  PV.)  It  would  therefore  be  consistent  to  require 
caused  by  instability.  ( In  the  mesosphere,  breaking  in-  7  »  1  -2  ( X 1 0  '5  s  ' 1 )  for  instability  ( cf.  Hitchman  and 
temal  gravity  waves  will  help  stabilize  the  flow.)  Other  Leovy  1986;  Hitchman  et  al.  1987).  2)  Realistic  basic 
factors  may  be  important.  1 )  Stability  depends  on  basic-  states  contain  a  continuous  spectrum  of  horizontal  and 
state  transience;  instabilities  must  grow  rapidly  com-  vertical  motions.  The  initial  evolution  need  not  be 
pared  to  mean  flow  change.  In  the  tropical  mesosphere,  dominated  by  the  most  unstable  mode;  the  instability 
planetary  waves  will  modify  the  basic  state  on  a  time  at  finite  amplitude  need  not  be  modal.  3 )  The  merid- 
scale  of  1-5  days.  (A  Lagrangian  time  scale  is  what  ional  circulation  itself  has  a  weak  stabilizing  influence 

at  large  e  (section  3f). 


f  Effect  of  meridional  circulation 

Cross-equatorial  shear  in  the  middle  atmosphere  is 
due  to  diabatic  advection  by  a  mean  meridional  cir¬ 
culation.  Including  in  (2.3)  a  zonally  symmetric  com¬ 
ponent  S0,  independent  of  latitude,  had  a  minor  effect 
on  local  stationary  instability  as  shown  in  Fig.  10.  For 
this  series  7  =  4*  10~5  s'1,  e'1  =  256  m2  s'2,  yc  =  1.8y.s, 
^  =  1,  and  «o  =  0.  The  location  of  instability  was  shifted 
north  and  east  when  fi0  >  0  (not  shown).  Growth  rate 
was  slightly  faster  in  most  cases  with  e~‘  »  256  m2  s'2 
plotted  in  Figs.  3,  5,  7,  and  9  when  S0  =  5  m  s'1 .  At 
large  e  there  was  a  reduction  of  growth  suggesting  the 
possibility  of  weak  scale  selection  when  v0  >  0.  The 
physical  explanation  may  be  that  v0  advects  the  insta¬ 
bility  away  from  maximum  R,  the  effect  being  more 
dramatic  for  high-*  modes  with  short  meridional  scale. 


log2  s 

Fig  .  9.  Growth  rate  as  a  function  of  6. 
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Fig.  10.  Growth  rate  as  a  function  of  C0- 

g.  Results  for  weak  shear 

The  preceding  results  raise  two  questions  about  in¬ 
stability  in  nonparallel  flow.  1 )  Observe  from  Figs.  4 
and  6  that  local  instability  had  the  same  zonal  scale  as 
basic-state  variation.  There  was  at  most  one  stationary 
or  two  propagating  cells  in  the  ,x  direction.  In  general, 
we  expect  multicell  nonsymmetric  instability  in  locally 
parallel  flow.  Under  what  circumstances,  then,  can  a 
wave  train  of  many  cells  develop  in  the  unstable  part 
of  the  flow?  2)  Advected  instability  was  observed  at 


large  !  z70 1  because  the  flow  was  unstable  at  all  ,v  and 
periodic  boundary  conditions  were  assumed.  Real  flows 
could  be  stable  at  some  longitude.  Is  advected  instability 
possible  in  this  case?  The  first  question  will  be  addressed 
now  and  the  second  in  the  next  subsection. 

It  was  seen  in  parallel  flow  that  the  minimum  zonal 
scale  for  latitudinal  trapping  of  eigenfunctions  (decay 
at  large  \y\)  is  defined  by  a  short-wave  cutoff  kmMy 
=  j3.  The  argument  could  be  extended  to  nonparallel 
flow;  minimum  zonal  scale  depends  on  shear.  In  the 
case  illustrated  above,  it  was  impossible  to  accomodate 
several  cells  in  the  x  direction.  For  multiple  cells  y 
must  be  reduced.  [That  is.  in  the  configuration  (3.1 ). 
There  may  exist  situations  where  y  is  unimportant 
outside  the  region  of  anomalous  PV,  such  that  trapping 
of  eigenfunctions  is  ensured  by  the  basic  state  alone, 
apart  from  k.] 

Figure  1 1  shows  an  example  in  which  y  =  0.5 
X  lO'3  s'1,  e  =  4  m~2  s2,  yc  =  1.8ys,  f  =  1,  and  Mq 
=  0.  Numerical  integration  used  32  harmonics  (note 
change  of  axes);  results  are  shown  at  240  days.  For 
this  value  of  7  a  long  time  was  required  to  obtain  the 
mode,  due  to  the  assumed  7  and  wave  1  initial  con¬ 
dition.  (Less  time  would  have  been  required  from  high- 
k  initial  conditions.)  This  example  demonstrated  that 
multicell  instability  is  possible  when  the  scale  of  mean- 
flow  variation  is  much  greater  than  k„ Lx  =  y/0.  An¬ 
other  example  of  small  7  is  discussed  in  the  Appendix. 

What  about  the  opposite  situation,  when  the  unstable 
region  is  smaller  than  Profiles  of  the  form  (3.1 ) 


0 


LONGITUDE 

Fig.  1 1.  Structure  of  multicell  instability  in  weak  shear  7  =  0.5  X  10"5  s'1. 
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could  be  generalized  in  order  to  coniine  the  most  un¬ 
stable  flow: 

-cos(x:/a)  =  2sin2(x/2a)-  1  -*-2sin',(.\Y2a)  -  1. 

(3.3) 


interpretation  is  that  ii0  displaced  the  center  of  mode 
downstream  from  the  source  of  instability  near  ,v 
=  180°.  For  sufficiently  large  1 170|  the  mode  was  pushed 
into  the  sponge  region  and  could  not  grow.  Evidently 
global  modes  %vere  stabilized  by  the  sponge  region  also. 


Setting  n-  16.  for  example,  the  region  of  most  unstable 
flow  was  less  than  60°  wide.  Results  indicated  that,  for 
parameters  of  Fig.  4,  the  local  stationary  mode  was 
essentially  unchanged  by  reducing  the  size  of  the  un¬ 
stable  region,  although  growth  rate  was  reduced  slightly 
(by  about  25%  when  n  —  16).  (An  explicit  version  of 
the  time-dependent  code  was  used  to  obtain  this  result 
and  to  validate  other  results  of  the  semi-implicit 
model.)  This  supports  its  interpretation  as  a  local  mode 
that  depends,  not  on  global  average  stability,  but  on 
local  stability  within  the  most  unstable  part  of  the  do¬ 
main  ( Pierrehumbert  1984). 

To  be  sure,  the  concept  of  “local  stability”  is  im¬ 
precise  when  ( as  in  this  case)  instability  and  basic  state 
share  the  same  zonal  scale.  WKJB  analysis  like  that  of 
Pierrehumbert  ( 1984)  is  formally  valid  only  in  slowly 
varying  mean  flow.  Perhaps  a  more  accurate  measure 
of  stability  would  be  obtained  by  averaging,  say,  over 
a  half-wavelength  of  perturbation. 

h.  Sponge  region  in  x 

As  mentioned  in  section  2c,  the  advected  instability 
was  altered  by  a  sponge  region  designed  to  completely 
absorb  disturbances  crossing  the  domain  boundary 
near  a:  =  0.  Though  artificial,  the  sponge  may  simulate 
regions  of  stable  flow  (as  undoubtedly  exist  in  the  at¬ 
mosphere  )  or  the  tendency  of  instabilities  to  “break” 
and  dissipate  through  nonlinear  saturation  before 
making  a  complete  circuit  of  the  globe. 

Figure  12  shows  growth  rate  as  in  Fig.  5  but  with 
sponge  region  included.  Instabilities  were  unaffected 
at  small  |  «bl  but  were  stabilized  at  large  ]  z70l-  The 
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Fig.  12.  Growth  rate  as  a  function  of  fto  as  in  Fig.  5. 
but  with  sponge  region  included. 


4.  Discussion 

Results  of  section  3  documented  local  as  well  as 
global  modes  of  inertial  instability  in  nonparallel  flow. 
Local  modes  dominate  inhomogeneous  flow  and  are 
likely  more  important  in  observed  basic  states.  Local 
stationary  instabilities  had  the  largest  growth  rate,  de¬ 
pending  primarily  on  the  local  rather  than  global  sta¬ 
bility  of  the  basic  state:  they  occurred  when  (among 
other  things )  mean  flow  advection  was  small  in  the 
unstable  region.  Large  advection  caused  zonal  phase 
propagation,  here  referred  to  as  advected  instability. 
Like  stationary  modes,  advected  modes  had  zero  group 
velocity  but  the  center  of  the  wave  packet  was  displaced 
downstream  from  the  most  unstable  longitude,  and  the 
growth  rate  reduced.  If  ii0  was  sufficiently  large  and  the 
flow  stable  at  large  \x\,  advected  modes  were  stabilized. 

Even  for  a  simple  profile  like  (3.1)  the  parameter 
dependence  of  inertial  instability  is  complicated;  ex¬ 
amples  of  section  3  highlighted  behavior  along  a  few 
trajectories  in  parameter  space.  One  of  the  more  in¬ 
teresting  results  was  the  dependence  of  local  instability 
on  Wo,  illustrating  the  breakdown  of  Galilean  invariance 
in  nonparailel  flow.  Zonal  inhomogeneity  breaks  the 
symmetry  of  uniform  parallel  flow  and  introduces  a 
new  constant:  the  zonal  propagation  speed  of  basic- 
state  pattern.  This  was  zero  in  (3. 1 );  consequently,  in¬ 
stabilities  with  zero  group  velocity  remained  in  the 
same  location  relative  to  the  mean  flow  pattern  at  all 
times.  (Not  surprisingly,  their  growth  maximized  when 
the  mode  best  overlapped  the  most  unstable  region.) 
In  WKB  theory,  a  flow  is  said  to  be  absolutely  unstable 
when  such  instabilities  exist,  that  is, 


where 


Uj(ko)  >  0 

(4.1) 

da) 

!kk  ~  0 

(4.2) 

and  unstable  branches  originate  on  opposite  sides  of 
the  real-A;  axis  (e.g.,  Pierrehumbert  1984,  1986).  The 
importance  of  absolute  instability  derives  from  sym¬ 
metry  breaking,  for  example,  in  nonparailel /nonuni¬ 
form  flow,  or  when  there  is  a  “wave  maker”  at  some 
.r(Huerre  and  Monkewitz  1990). 

Pierrehumbert  ( 1 984)  analyzed  the  two-layer  quasi- 
geostrophic  model  for  baroclinic  instability  of  non¬ 
parallel  flow,  including  a  discussion  of  “absolute”  and 
“convective”  instability  in  terms  of  local  and  global 
modes,  supplemented  by  WKB  analysis.  It  is  unnec¬ 
essary  to  duplicate  this  work,  but  simple  to  quote  three 
of  his  main  conclusions  pertaining  to  local  modes: 
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Growth  rate  of  a  local  mode  is  determined  by  the  max¬ 
imum  baroclinicity  in  the  domain,  and  not  by  the  av¬ 
erage  baroclinicity. 

This  was  suggested  for  inertially  unstable  flow  ( sub¬ 
stituting  -‘inertial  instability”  for  “baroclinicity”)  by 
results  of  Figs.  5  and  12  demonstrating  the  insensitivity 
of  local  stationary  instability  to  a  sponge  region  near 
x  =  0,  and  by  discussion  at  the  end  of  section  3g  in 
which  growth  rate  was  only  weakly  dependent  on  the 
size  of  the  unstable  region.  As  in  Pierrehumbert  ( 1984), 
the  maximum  baroclinicity  or  maximum  inertial  in¬ 
stability  of  the  basic  state  must  take  into  account  factors 
such  as  boundary  conditions,  pressure  gradients,  and 
Ho  in  addition  to  the  parcel  growth  rate  (e.g.,  R).  On 
this  account  we  are  reminded  that  equatorial  inertial 
instability  is  not  a  pure  parcel  instability  (Dunkerton 
1983). 

Vanishing  absolute  growth  rate  at  infinity  is  not  nec¬ 
essary  for  localization. 

This  was  apparent  from  examples  based  on  the  pro¬ 
file  ( 3. 1 ),  which  was  unstable  at  all  x  but  displayed  a 
preference  for  local  instabilities  in  most  circumstances. 

The  contrast  between  maximum  and  minimum  baro¬ 
clinicity  in  the  flow  determines  the  extent  of  localiza¬ 
tion  of  the  eigenmodes.  with  high  contrast  favoring 
localization. 

In  a  similar  way  ( although  it  was  not  discussed  ex¬ 
tensively  above),  the  zonal  scale  of  local  instability  in¬ 
creased  as  $  was  decreased  below  1,  before  the  onset 
of  global  instability. 

Absolute  instability  does  not  depend  on  the  validity 
of  WKB  in  (4.1,  2)  since,  by  definition,  the  instability 
must  grow  indefinitely  at  a  fixed  point  in  x.  This  prop¬ 
erty  was  demonstrated  numerically,  but  not  analyti¬ 
cally,  for  local  stationary  instability.  WKB  analysis  of 
inertial  instability  on  an  equatorial  j3  plane  is  difficult, 
although  suggestive  of  absolute  instability  as  discussed 
in  the  Appendix.  Basic  states  like  Fig.  2  disallow  a  clean 
separation  of  wave  and  mean-flow  scales,  and  a  nu¬ 
merical  demonstration  of  instability  is  preferable  in 
this  case  (cf.  Jones  and  Thorpe  1992). 

5.  Conclusions 

A  simple  model  has  shown  that  inertial  instabilities 
in  zonally  nonuniform  cross-equatorial  shear  take  the 
form  of  local  and  global  modes.  Local  modes  may  be 
exactly  stationary  or  display  zonal  phase  propagation 
depending  on,  among  other  things,  advection  by  the 
basic  state.  Both  subtypes  of  local  instability  are  dis¬ 
tinguished  from  global  instability  by  their  zero  group 
velocity  and  concentration  of  amplitude  within,  or 
downstream  from,  the  region  of  most  unstable  flow. 
These  properties,  including  breakdown  of  Galilean  in¬ 
variance,  are  reminiscent  of  absolute  instability — al¬ 
though,  apart  from  the  Appendix,  none  of  our  results 
depended  on  a  WKB  approximation.  Exact  stationarity 


of  phase  over  a  finite  range  of  parameters,  and  mea¬ 
surable  dependence  of  growth  rate  on  the  size  of  un¬ 
stable  region,  are  possibly  attributable  to  a  non-WKB 
effect. 

It  was  suggested  that  for  equatorial  Rossby  wave 
critical  layers  in  the  mesosphere,  local  stationary  modes 
are  important  given  the  strength  of  background  shear 
7  and  weakness  of  in  situ  advecting  current  near  the 
zero-wind  line.  Such  a  simplification  is  desirable  in  this 
context  to  avoid  the  complexity  of  global  and  advected 
modes.  Existence  of  local  stationary  modes  locked  to 
PV  anomalies  is  consistent  with  middle  atmosphere 
simulations  of  inertial  instability  and  Rossby  wave 
breaking  by  O’Sullivan  and  Hitchman  (1992).  As  a 
follow-up,  it  will  be  worthwhile  to  show  the  evolution 
of  unstable  modes  in  an  authentic  Rossby  wave  critical 
layer,  and  their  effect  on  potential  vorticity  evolution. 
This  will  be  done  in  the  sequel. 

Local  and  global  modes  of  inertial  instability  share 
many  of  the  same  properties:  for  example,  dependence 
of  growth  rate  on  vertical  wavelength  and  downgradient 
momentum  flux.  Vertical  scale  selection  may  be 
achieved  with  diffusion  as  argued  previously,  or,  as  the 
new  analysis  suggests,  through  mean  meridional  ad¬ 
vection.  Nonmodal  instability  (Farrell  1982)  on  a 
transient  basic  state  will  allow  scale  selection  at  finite 
amplitude.  It  should  be  obvious  from  our  results  that 
when  the  time  scale  of  basic-state  variation  is  much 
less  than  y  1 ,  exponential  modes  of  inertial  instability 
are  irrelevant. 

Inertial  instability  can  be  expected  to  have  three 
consequences  in  the  tropical  middle  atmosphere:  hor¬ 
izontal  redistribution  of  angular  momentum  (and  po¬ 
tential  vorticity),  dissipation  of  laterally  propagating 
Rossby  waves,  and  vertical  mixing  of  constituents.  Be¬ 
cause  of  their  similarity,  local  and  global  modes  may 
contribute  alike  to  mean  flow  effects,  such  as  inertial 
adjustment  ( Dunkerton  1981)  and  mixing.  [The  effi¬ 
ciency  of  inertial  adjustment  is  open  to  question.  Our 
recent  numerical  results  in  the  middle  atmosphere 
demonstrate  that  an  adjustment  occurs  but  is  incom¬ 
plete.  See  also  Thorpe  and  Rotunno  ( 1989)  for  further 
discussion  in  the  context  of  frontal  CSI.]  Local  modes 
will  be  more  important  for  Rossby  wave  dissipation. 
According  to  Killworth  and  McIntyre  (1985),  baro- 
tropic  processes  are  insufficient  for  time-averaged  crit¬ 
ical  layer  absorption  as  t  -*  oo . 
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APPENDIX 

Dispersion  Relation  at  Complex  k 

Analysis  of  absolute  instability  in  multidimensional 
flow  is  difficult:  a  formal  proof  of  absolute  equatorial 
inertial  instability  using  WKB  theory  has  not  been 
given  and  will  not  be  attempted  here.  (I  doubt  that 
such  a  proof  is  possible  for  equatorial  instability  without 
approximation  to  the  governing  equations.)  It  is  inter¬ 
esting  nevertheless  to  calculate  eigenfrequencies  of 
(2.10)  at  complex  k.  Several  limitations  of  this  ap¬ 
proach  should  be  noted.  1 )  Analysis  begins  with  the 
complex  dispersion  relation  defined  numerically  by 
(2.10)  with  2)  =  0  rather  than  the  Laplace  transform 
of  the  time-dependent  system;  this  excludes  continuum 
modes  that  may  be  important  in  reality.  2)  The  exact 
dispersion  relation  is  determined  numerically  over  a 
finite  range  of  k,  rather  than  analytically  over  all  k. 
The  cubic  approximation  of  section  2c  (and  a  quadratic 
version  of  it)  gave  similar  results,  but  their  validity 
could  not  be  guaranteed  for  all  combinations  e,k.  3) 
The  infinite  0  plane  could  not  be  simulated;  results 
were  obtained  in  a  channel.  Channel-dependent  modes 
occur  outside  the  short-wave  cutoff  (defined  by  a  hy¬ 
perbola  in  the  complex-^  plane  extending  outward 
from  kmax  on  the  real  axis).  This  region  will  be  ignored. 
4)  Construction  of  WKB  solutions  in  some  cases  may 
require  matching  across  a  WKB  breakdown  point 
( Pierrehumbert  1984);  the  full  construction  will  not 
be  attempted  but  it  will  be  assumed  that  matching  is 
possible. 

To  simplify  matters,  the  shear  was  assumed  inde¬ 
pendent  of  latitude: 


U  -  y(y  —  yc) 


1  —  yp  cosx/a 
l  +  ip 


=  y(y  -  yc)%( x) 


and  V  m  0,  corresponding  as  closely  as  possible  to  the 
theoretical  results  of  section  2c.  Local  stationary  insta¬ 
bility  exists  near  yc  =  ys;  discussion  will  center  around 
this  choice  of  yc.  The  eigenproblem  (2.10)  was  solved 
at  x/ a  =  ir,  the  location  of  maximum  shear  (where  'fy 
=  1 ),  allowing  A:  to  be  complex.  Eigenfrequencies  were 
obtained  by  a  shooting  method  for  several  combina¬ 
tions  y,  e.  Results  shown  in  Fig.  13  were  representative 
when  e  ~  €neut  (in  this  case,  y  -  1.5  X  10-5s“‘,  e 
=  0.1657  m-2  s2).  Growth  rate  increased  in  both  di¬ 
rections  away  from  the  real  axis,  and  there  was  a  saddle 
point  near  s  =  ka  =  (5.55,  -1.50)— not  far  from  the 
most  unstable  wavenumber  on  the  real  axis.  Existence 
and  location  of  the  saddle  are  a  function  of  parameters 
( y,  e,  'Iff,  yc,  etc.)  so  there  is  nothing  special  about  this 
example.  In  the  inversion  of  Laplace  transform,  the 
saddle  is  considered  unavoidable  when  integrating  over 
kr  if  unstable  branches  originate  on  opposite  sides  of 
the  real-A;  axis  ( Pierrehumbert  1 986 ) .  This  was  appar¬ 
ently  the  case  based  on  numerical  results,  but  an  an¬ 
alytical  proof  is  beyond  the  scope  of  this  paper. 


Not  all  saddles  are  unavoidable  in  this  problem,  par¬ 
ticularly  those  on  the  imaginary  axis  resulting  from 
merger  of  unstable  branches  originating  in  the  upper 
half-plane.  This  sort  of  behavior  was  observed  when  u 

=  yy- 

Knowledge  of  £(A)  is  sufficient  to  determine  new 
eigenfrequencies  w(k)  from  a  complex  Doppler  shift 

u(k)  =  u>(k)  +  ku0 

where  £  is  given  by  the  complex  dispersion  relation  in 
the  control  case  yc  =  ys  (e.g.,  Fig.  13).  Other  effects 
could  be  added,  such  as  x  diffusion.  In  WKB  theory, 
breakdown  of  Galilean  invariance  is  evidently  due  to 
the  imaginary  part  of  k,  as  can  be  seen  by  translating 
a  wave  packet  of  exponential  shape  to  the  left  or  right. 
For  constant  shear,  variation  of  yc  introduces  a  constant 
mean  flow  change  u0  =  -y&yc%.  The  location  of  sad¬ 
dle  point  and  its  associated  eigenfrequency  change  as 
a  function  of  yc ;  this  is  illustrated  in  Fig.  14  where 
comparison  is  made  to  time-dependent  model  results 
obtained  with  t p  =  0.25.  The  agreement  is  rather  good. 
Growth  maximized  when  the  mean  flow  was  zero  near 
the  center  of  the  unstable  region.  Due  to  noticeable 
asymmetry  of  £(k)  between  upper  and  lower  half¬ 
planes,  results  were  not  exactly  symmetric  about  vc/ys 
=  1.0. 

Growth  rates  were  slightly  smaller  in  the  time-de¬ 
pendent  model,  by  a  few  percent — possibly  due  to  y 
diffusion,  included  in  the  time-dependent  model  but 
not  in  the  eigenproblem  (2.10).  (Horizontal  diffusivity 
in  the  time-dependent  model  was  approximately  2. 1 


co,  / 7 


kr 

Fig.  13.  Inviscid  growth  rate  as  a  function  of  complex  k  obtained 
from  the  eigenproblem  (2.10)  with  2)  =  0.  u  -  y(  v  -  yc)y  y  =  1.5 
X  10~5s“\  e  =  0.1657  m~2  s2,  and  yc  ~  ys.  (For  clarity,  not  all 
contours  are  shown. ) 
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X  10 5  nr  s~! .  Saddle-point  values  in  Fig.  14  were  cor¬ 
rected  for  .v  diffusion  only.)  There  was  also  a  non-WKB 
effect  such  that  growth  rates  in  the  time-dependent 
model  diminished  as  \p  increased,  that  is,  as  the  size  of 
the  unstable  region  was  reduced  in  .v.  Simulations  with 

^  0.25  underestimated  the  true  growth  rate  some¬ 
what. 

Calculated  and  observed  periods  were  in  reasonable 
agreement  for  advected  instabilities.  Note  that,  unlike 
the  situation  in  Fig.  3,  ur  was  nonzero  everywhere  ex¬ 
cept  near  yc/ys  —  1 .0;  this  was  due  to  the  smaller  value 
of  7  and  relatively  smaller  contribution  from  the  zon- 
ally  symmetric  component  of  instability. 

Displacement  of  the  saddle  point  off  the  real  axis 
implies  growth  of  the  wave  packet  between  x/ a  =  x 
and  its  center,  which  is  displaced,  in  physical  space, 
downstream  from  the  region  of  most  unstable  flow. 
This  expectation  was  verified  by  numerical  results,  in¬ 
cluding  a  more  subtle  prediction  that  zonal  wavelengths 
are  slightly  larger  (smaller)  for  packets  displaced  to  the 
left  (right)  of  center. 

Finally,  it  was  possible  to  vary  x/an  holding  other 
parameters  fixed,  and  determine  At(oj,  a:)  with  cu  set 
equal  to  its  saddle-point  value  as  in  Fig.  10  of  Pierre- 
humbert  ( 1984).  This  was  a  tedious  calculation  due 
to  rapid  variation  of  k;  only  two  examples  were  at¬ 
tempted.  one  at  eneut  and  another  above.  In  both  cases, 
it  was  striking  how  quickly  k  departed  from  the  real 
axis  in  both  directions  as  'Py  was  reduced  below  unity. 
This  may  help  explain  the  locally  symmetric  appear¬ 
ance  of  local  instability  having  few  cells  in  the  x  direc¬ 
tion.  Also,  it  may  explain  why  local  modes  can  be  found 
in  weakly  inhomogeneous  flow,  as  only  a  slight  reduc¬ 
tion  of  'P/  below  unity  immediately  moves  k  into  the 
opposite  half-plane — a  necessary  condition  for  the  ex¬ 
istence  of  local  modes  (Pierrehumben  1984). 


When  6  >  €neut,  the  saddle-point  trajectory  is  closer 
to  the  imaginary  axis,  consistent  with  numerical  results 
demonstrating  the  preference  for  quasi-stationary,  sin¬ 
gle-cell  instability  in  this  case. 
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ABSTRACT 

The  seasonal  ( wintertime )  development  of  middle  atmosphere  circulation  in  opposite  phases  of  the  equatorial 
quasi-biennial  oscillation  (QBO)  was  simulated  with  a  three-dimensional  nonlinear  numerical  model.  In  the 
stratosphere,  the  effect  of  equatorial  QBO  was  generally  consistent  with  the  extratropical  QBO  observed  by 
Holton  and  Tan,  namely,  a  stronger  midwinter  polar  vortex  in  the  westerly  phase,  and  vice  versa.  However,  the 
extratropical  response  to  the  QBO  was  sensitive  to  other  factors  such  as  mesospheric  gravity  wave  drag  and  the 
amplitude  of  Rossby  waves  specified  at  the  model’s  lower  boundary.  The  extratropical  QBO  was  realistic  only 
when  a  drag  parameterization  was  included  and  Rossby  wave  amplitudes  lay  in  an  intermediate  range  close  to 
the  observed.  At  somewhat  stronger  forcing,  the  model’s  response  was  largest  in  the  mesosphere  where  (in  this 
effect WeSterliCS  WerC  Str°nger  in  **  eaSterly  Phase  of  equatorial  QBO.  This  was  apparently  due  to  a  shielding 

The  theory  of  planetary  wave-mean  flow  interaction  suggests  that  the  sensitivity  to  equatorial  QBO  should 
be  greatest  for  wave  forcings  near  a  “bifurcation”  point.  Below  this  threshold  the  stratosphere  approaches 
radiative  equilibrium,  shutting  off  vertical  propagation  of  planetary  waves.  Supercritical  forcing  leads  to  a  major 
warming.  The  model’s  sensitivity  to  forcing,  while  consistent  with  this  idea,  was  most  apparent  in  perpetual 
solstice  runs  without  parameterized  wave  drag.  Seasonal  integrations  with  wave  drag  produced  a  more  realistic 
extratropical  QBO.  making  the  bifurcation  less  conspicuous. 


1.  Introduction 

It  is  currently  thought  that  the  equatorial  quasi-bien¬ 
nial  oscillation  ( QBO )  affects  the  wintertime  circulation 
of  the  extratropical  middle  atmosphere.  The  strato¬ 
spheric  polar  vortex  is  weaker  and  more  disturbed  by 
planetary  waves  in  the  easterly  phase  of  the  QBO  (Hol¬ 
ton  and  Tan  1980,  1982;  Wallace  and  Chang  1982;  La- 
bitzke  1982,  1987;  van  Loon  and  Labitzke  1987;  Dunk¬ 
erton  and  Baldwin  1991;  Baldwin  and  Dunkerton 
1991).  During  the  westerly  phase  there  is  a  colder, 
stronger  vortex  with  greater  potential  for  ozone  deple¬ 
tion  (Bojkov  1986;  Garcia  and  Solomon  1987;  Lait  et 
ai.  1989).  In  the  Northern  Hemisphere,  this  “extratrop¬ 
ical  QBO’  ’  is  strongest  in  early  to  midwinter  ( Dunkerton 
and  Baldwin  1991,  1992).  It  is  not  uniform  in  all  years 
but  ( for  reasons  not  well  understood)  seems  more  prom¬ 
inent  near  solar  minima  (Labitzke  and  van  Loon  1988) 
in  late  winter  (Dunkerton  and  Baldwin  1992). 

Planetary  waves  play  an  important  role  in  the  extra¬ 
tropical  QBO  as  measured  by  their  Eliassen-Palm  flux 
and  divergence  (Dunkerton  and  Baldwin  1991 — their 
Fig.  16).  Whether  they  are  the  main  cause  of  the  ex- 


Corresponding  author  address:  Dr.  Timothy  J.  Dunkerton.  North¬ 
west  Research  Associates.  P.O.  Box  3027.  Bellevue,  WA  98009- 
3027. 


©  1994  American  Meteorological  Society 


tratropical  QBO  is  not  readily  answered  from  obser¬ 
vations.  but  is  suggested  by  recent  numerical  investi¬ 
gations  ( O’Sullivan  and  Salby  1990;  Dameris  and  Ebel 
1990;  Holton  and  Austin  1991;  O’Sullivan  and  Young 
1992).  Starting  from  initial  conditions  differing  only 
in  the  phase  of  equatorial  QBO  (with  planetary  Rossby 
waves  forced  at  a  tropopause  lower  boundary ) ,  twin 
experiments  gradually  diverge  over  several  weeks  until 
the  extratropical  mean  flow  is  weakened  by  wave  trans¬ 
port — the  overall  deceleration  being  somewhat  larger 
in  the  easterly  phase  of  the  QBO.  There  can  be  a  sig¬ 
nificant  extratropical  difference  induced  by  the  equa¬ 
torial  QBO  that  is  most  realistic  for  intermediate  wave 
forcings  (Holton  and  Austin  1991;  O’Sullivan  and 
Young  1992).  Large  forcings  cause  the  vortex  to  break 
down  prematurely  ( as  in  a  major  warming )  irrespective 
of  the  QBO,  while  weak  forcings  produce  a  weak  or 
negligible  effect. 

According  to  numerical  results,  it  is  not  the  instan¬ 
taneous  linear  difference  of  planetary  wave  structure 
due  to  equatorial  QBO  that  is  important,  but  rather  the 
cumulative  effect  of  planetary  wave— mean  flow  inter¬ 
action  over  a  period  of  time.1  The  model’s  extratropical 


1  Linear  structural  changes  are  possible,  but  their  effect  on  EP  flux 
divergence  is  confined  mainly  near  subtropical  critical  latitudes 
(O’Sullivan  and  Young  1992). 
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QBO  is  therefore  a  quasi-linear  or  nonlinear  phenom¬ 
enon  rather  than  a  linear  one. 

By  the  same  token,  the  extratropical  QBO  evolves 
during  winter.  The  mean  zonal  wind  difference  be¬ 
tween  west  and  east  phases  is  maximum  in  midwinter 
( December- January ) ,  while  that  of  the  EP  flux  max¬ 
imizes  a  bit  earlier  (November— December)  (Dunker- 
ton  and  Baldwin  1991 — their  Figs.  10,  20).  To  model 
the  seasonal  development  realistically  requires  time- 
dependent  forcings  and  diabatic  cooling  of  the  ob¬ 
served  magnitude.  O’Sullivan  and  Young  (1992)  ig¬ 
nored  the  seasonal  development  and  relaxed  the  flow 
to  an  equilibrium  zonal  wind  similar  to  the  observed 
(50%  stronger  than  their  initial  condition  in  midlati¬ 
tudes),  rather  than  radiative  equilibrium.  Consequently 
their  radiative  forcing  was  too  weak,  as  was  the  induced 
meridional  circulation.  Holton  and  Austin  ( 1991 )  used 
a  sophisticated  radiative  algorithm,  but  began  their  ex¬ 
periments  in  mid- January,  which  is  after  the  maximum 
of  observed  extratropical  QBO. 

One  purpose  of  this  study  was  to  model  the  seasonal 
development  of  extratropical  QBO  with  realistic  ther¬ 
mal  drive;  another  was  to  understand  its  sensitivity  to 
various  parameters,  such  as  the  amplitude  of  Rossby 
waves  entering  the  stratosphere.  Theory  suggests  that 
weak  waves  allow  the  mean  state  to  approach  radiative 
equilibrium,  impeding  further  wave  propagation  Con¬ 
versely,  strong  waves  decelerate  the  mean  flow,  en- 
nancing  propagation  until  the  vortex  experiences  a  ma- 
jor  warming.  Simple  one-dimensional  models  predict  a 
bifurcation”  between  extremes  that  is  a  function  of 
wave  forcing  amplitude  and  other  parameters  ( Holton 
and  Mass  1976;  Holton  and  Dunkerton  1978'  Yoden 
1987a.b,c,  1990). 

10  Ijght  of  the  ID  theory,  it  is  interesting  to  consider 
,ct  ot  some  external  influence,  like  the  equatorial 
QBO,  in  conjunction  with  planetary  waves  forced  at 
the  extratropical  tropopause.  If  the  wave  forcing  is  such 
that  for  one  phase  of  the  QBO  the  forcing  is  subcritical. 
while  for  the  opposite  phase  it  is  supercritical,  we  ex¬ 
pect  maximum  sensitivity  to  the  equatorial  QBO.  If  the 
forcing  is  above  or  below  this  threshold  in  both  phases 
alike,  we  expect  minimum  sensitivity.  This  behavior  is 
consistent  with  previous  numerical  studies  of  the  ex¬ 
tratropical  QBO  but  deserves  further  investigation  for 
example,  to  determine  whether  the  bifurcation  occurs 
in  a  nonlinear,  three-dimensional  model  of  the  middle 
atmosphere,-  and  what  effect  the  seasonal  cycle  has 
Other  factors  may  also  affect  the  extratropical  response 
to  the  QBO.  The  effect  of  gravity  wave  drag  on  plan- 
etaiyt  wave-mean  flow  interaction  will  be  emphasized 
in  the  discussion  to  follow. 


Section  2  briefly  describes  the  numerical  model  and 
the  parameterization  of  thermal  and  mechanical  forc¬ 
ings.  Section  3  describes  the  seasonal  mean-flow  evo¬ 
lution  with  and  without  parameterized  wave  drag,  and 
the  variation  of  extratropical  QBO  as  a  function  of 
wave  forcing  in  a  baseline  series  of  experiments.  These 
results  are  discussed  further  in  section  4  along  with 
possible  causes  of  extratropical  QBO  variability. 

2.  Numerical  model 

The  numerical  model  used  in  this  study  was  the  same 
as  that  of  O’Sullivan  and  Young  (1992):  a  global, 
primitive  equation  model  extending  from  10  to  70  km 
with  geopotential  specified  at  the  lower  boundary.  Sim¬ 
ulations  were  performed  at  T42/2 1  horizontal  resolu¬ 
tion  ( spherical  harmonics  with  total  wavenumber  up  to 
42  and  zonal  wavenumber  up  to  21  are  included)  with 
20  equally  spaced  levels  in  log  p(Az  =  3  km).  Sixth- 
order  horizontal  diffusion  was  included  to  suppress 
grid-scale  motions  {v  =  2.5  x  102s  m6  s'1).  The  time 
step  was  15  minutes.  Time  integrations  began  on  “11 
September”  and  extended  thereafter  for  200  days.  The 
mean  thermal  forcing  and  perturbation  lower-boundary 
condition  varied  through  the  seasonal  cycle  as  de- 
scribed  below. 

The  main  differences  relative  to  O’Sullivan  and 
Young  (1992)  were  the  parameterization  of  thermal 
and  mechanical  drive,  modified  to  produce  a  more  re¬ 
alistic  wintertime  circulation,  and  integration  through  a 
seasonal  cycle.  Radiative  heating  and  unresolved  body 
forces  were  parameterized  as  linear  relaxation.  This  ap¬ 
proach  allowed  many  integrations  at  relatively  low 
cost. 

a.  Thermal  and  mechanical  forcings 

Zonally  averaged  diabatic  heating  was  of  the  form 

Q(y,  z.  t)  —  ocT(z)[fe(y ,  -,  t)  -  T(v,  t )],  (2.1 ) 

where  fE  is  the  “radiatively  determined”  profile  of 
Shine  ( 1987) .  This  is  not  the  true  radiative  equilibrium 
under  perpetual  conditions  but  a  radiatively  determined 
temperature  in  the  seasonal  cycle.  When  used  in  con¬ 
junction  with  a  temporally  fixed  Newtonian  coolin'7 
rate,  * 


ar(z) 


=  0.09  +  0.04  tanh  (  — -35j  J 


day"1,  (2.2) 


-Holton  and  Wehrbein  (1981)  demonstrated  that  vacillation  was 
possible  in  a  3D  quasi-linear  model,  and  nonlinear  barotropic  inte¬ 
grations  showed  evidence  of  bifurcation  (Salby  et  al  1990-  Yoden 
and  Ishioka  1993).  n 


realistic  solstitial  heating  and  cooling  are  obtained  This 
was  verified  by  comparison  to  K.  Shine’s  MIDRAD 
algorithm  in  a  2D  model.  Newtonian  cooling  rates  im¬ 
plied  by  that  algorithm  vary  in  latitude  and  time  an 
effect  which  is  ignored  in  (2.2)  but  is  implicit  in  our 
use  of  radiatively  determined,  rather  than  radiative 
equilibrium,  temperature.  Identical  damping  was  ap¬ 
plied  to  the  deviation  from  zonal-mean  temperature. 
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The  temporal  dependence  of  radiative  forcing  was  fur¬ 
ther  simplified  by  adopting  only  the  solstitial  profile  of 
Te,  holding  its  symmetric  part  about  the  equator  fixed 
while  letting  the  antisymmetric  part  vary  sinusoidally 
through  the  annual  cycle. 

As  noted  by  Shine  ( 1987)  and  many  others,  the  ra- 
diatively  determined  state  in  the  middle  atmosphere 
yields  an  enormous  zonal  wind  at  winter  solstice.  To 
overcome  this  problem  requires  body  forces  that  retard 
the  mean  flow.  Two  kinds  of  body  force  can  play  a 
role,  (i)  Planetary  Rossby  waves  forced  at  the  model’s 
lower  boundary  (section  2b)  decelerate  the  polar  vor¬ 
tex.  These  waves  can  have  a  large  effect,  but  (in  the 
absence  of  other  body  forces)  their  ability  to  retard  the 
vortex  turns  out  to  be  quite  sensitive  to  forcing  ampli¬ 
tude  ( section  3a ) .  ( ii )  Unresolved  disturbances  such  as 
gravity  waves  exert  a  substantial  drag  on  the  mean 
flow,  particularly  in  the  mesosphere. 

An  investigation  of  gravity  waves  in  the  middle  at¬ 
mosphere’s  three-dimensional  circulation  would  be  in¬ 
teresting  (Dunkerton  and  Butchart  1984)  but  outside 
the  scope  of  this  paper.  It  will  be  assumed  instead  that 
the  cumulative  effect  of  unresolved  waves  is  to  retard 
the  zonal-mean  westerlies  to  something  like  their  cli¬ 
matological  value.  This  was  achieved  in  the  numerical 
model  by  inserting  in  the  zonally  averaged  zonal  mo¬ 
mentum  equation  a  one-way  friction  ( Dunkerton  1991) 
designed  to  relax  overly  strong  westerlies  to  the  cli¬ 
matological  average  profile  of  mean  zonal  wind  uM 
(Barnett  and  Comey  1985)  if  and  only  if  u>  uM.  Pro¬ 
files  of  monthly  mean  uM  were  derived  from  Barnett 
and  Comey  for  each  of  the  12  calendar  months,  and 
interpolated  to  the  current  time  step.  The  mechanical 
damping  coefficient  was  zero  for  winds  below  25 
ms"1,  increasing  to  0.5  per  day  for  winds  in  excess  of 
50  ms'1.  One-way  friction  had  little  effect  on  the  polar 
lower  stratosphere  and  no  direct  effect  in  the  Tropics 
or  summer  hemisphere,  but  severely  restricted  the  me¬ 
sospheric  jet  and  thereby  maintained  a  strong  diabatic 
circulation  at  upper  levels.  Hereafter  the  friction  will 
be  referred  to  as  "parameterized  wave  drag.” 

Although  this  device  was  used  mainly  for  simplicity, 
gravity  wave  momentum  flux  divergence  can  be  written 
as  linear  relaxation  when  momentum  flux  is  propor¬ 
tional  to  the  negative  of  mean  flow  (as  suggested  by 
some  observations),  and  the  mean  flow  varies  slowly 
over  a  density  scale  height.  Similar  results  might  be 
achieved  with  a  more  sophisticated  gravity  wave  pa¬ 
rameterization  (Fritts  and  Lu  1993)  tuned  to  produced 
realistic  winds. 

b.  Rossby  wave  forcing 

Geopotential  was  specified  at  the  lower  boundary 
( 10  km) .  The  forcing  consisted  solely  of  stationary  zo¬ 
nal  wavenumber  1  with  maximum  amplitude  at  60°N 
as  in  Eq.  (2)  of  O’Sullivan  and  Young  ( 1992).  It  was 
turned  on  gradually  over  the  first  ten  days  to  an  initial 
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value  150  m  at  r  =  r„  ( 3  September  21 ).  Thereafter, 
a  slow  seasonal  amplification  was  assumed 


d>i 


60N 


1502  +  b1  sin 


ttjt  -  fo) 
1 80  days 


(2.3) 


as  suggested  by  Randel’s  (1992)  climatology.  Here 
d>60N  is  the  maximum  value  of  wavenumber  one  ampli¬ 
tude  on  the  lower  boundary.  The  parameter  b  was  var¬ 
ied  from  one  experiment  to  the  next  to  investigate  the 
effect  of  forcing  amplitude  on  the  extratropical  QBO. 
Hereafter  we  denote  this  variation  in  terms  of  a  param¬ 
eter  <bm,  the  maximum  value  of  d»«,N  attained  at  winter 
solstice.  Values  of  <f>m  ranged  from  150-350  m.  like 
those  observed  ( ~  200 -280  m). 

Slowly  varying  forcing  is  an  idealization  of  the  real 
atmosphere  with  its  transient  waves.  Pulsations  of  wave 
activity  cause  episodes  of  enhanced  mixing  and  decel¬ 
eration  (O’Sullivan  and  Salby  1990),  but  so  long  as 
quasi-stationary  waves  dominate  the  spectrum,  the°ex- 
tratropical  QBO  should  be  much  the  same. 


c.  Equatorial  QBO 

The  model  was  initialized  with  an  equinoctial  zonal 
wind  profile  typical  of  September,  together  with  a  con¬ 
tribution  from  the  equatorial  QBO  like  that  of  O’Sul¬ 
livan  and  Young  (1992) — a  Gaussian  profile  in  lati¬ 
tude  and  height,  centered  on  the  equator  between  about 
20  and  40  km — but  (in  our  case)  having  an  e-folding 
width  of  10°  ( 15°)  latitude  and  maximum  amplitude  15 
ms"1  (-20  ms"1)  in  westerly  (easterly)  phase.  Nu¬ 
merical  experiments  showed  that  the  extratropical  QBO 
was  unrealistically  large  for  wider  equatorial  QBOs. 
Interestingly,  most  of  the  extratropical  QBO  was  due 
to  the  westerly  equatorial  phase;  results  in  easterly 
phase  did  not  differ  much  from  control  experiments 
with  no  QBO. 

West-phase  experiments  included  a  small  restoring 
force  to  prevent  erosion  of  equatorial  winds.  This  de~ 
vice  might  represent  the  positive  acceleration  associ¬ 
ated  with  equatorial  Kelvin  waves  not  included  in  the 
model. 


3.  Results 

a.  Experiments  without  wave  drag  parameterization 

O’Sullivan  and  Young  (1992)  relaxed  mean  tem¬ 
perature  to  an  equilibrium  state  similar  to  the  observed. 
Consequently  their  mean  zonal  wind  was  realistic  but 
diabatic  cooling  and  mean  meridional  circulation  were 
not.  Here  mean  temperature  was  relaxed  to  a  radiativelv 
determined  state  so  that,  without  dynamical  heat  fluxes 
or  opposing  body  forces,  the  polar  vortex  becomes  ex¬ 
tremely  strong.  This  “cold  bias”  can  be  alleviated  by 
transport  due  to  planetary  Rossby  waves,  gravitv 
waves,  or  both.  It  is  interesting  (and  perhaps  amusing  ) 
to  consider  first  the  situation  without  gravity  waves. 
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where  vertically  propagating  planetary  waves  alone  are 
responsible  for  mean  flow  deceleration.  Our  largest  ex- 
tratropicai  QBO  was  obtained  in  perpetual  solstice  runs 
(temporally  constant  TE,  <bm  after  the  initial  i 0-day 
ramp)  without  parameterized  wave  drag.  In  some  of 
these,  the  time-averaged  difference  of  mean  zonal  wind 
between  west  and  east  QBO  phases  exceeded  160 
m  s  m  the  mesosphere,  the  pattern  being  similar  in 
other  respects  to  that  of  O’Sullivan  and  Youn<*  ( l  qqo- 
their  Fig.  6a).  s  v  “’ 

Time  series  illustrating  the  QBO’s  effect  are  shown 
in  Fig.  1.  The  mean  zonal  wind  in  both  cases  acceler¬ 
ated  under  radiative  control  until  day  50,  after  which 
die  east,  but  not  west,  case  decelerated — gradually  at 
first,  and  then  more  rapidly.  The  magnitude  of  this  “ex- 
tratropical  QBO”  at  days  120-150  (ten  times  the  ob¬ 
served)  revealed,  in  essence,  the  difference  between  a 
radiauvely  determined  state  (west  phase)  and  a  major 
warming  (east  phase).  The  time-averaged  difference 
was  greatest  when  the  averaging  interval  was  located 
after  the  first  major  warming  in  the  east-phase  run,  and 
wave  forcing  was  optimally  chosen  to  lie  near  a  “bi¬ 
furcation  point  so  there  was  no  major  warming  in  the 
west-phase  run.  The  bifurcation  was  anticipated  by 
simple,  one-dimensional  models  (Holton  and  Mass 
1976;  Holton  and  Dunkerton  1978)  and  its  effect  in 
this  case  was  quite  dramatic.  In  east-phase  experiments 
(similar  to  those  with  no  equatorial  QBO),  wave  forc¬ 
ings  st  300  m  led  to  a  major  warming.  Smaller 


forcings  (e.g..  d>m  275  m)  could  not  prevent  the  ap- 
proach  to  radiative  equilibrium.  The  effect  of  the  QBO 
was  to  shift  the  bifurcation  point  toward  greater  <bm  in 
the  west  phase,  so  that  twin  experiments  with  identical 
wave  forcing  near  the  original  bifurcation  produced  a 
large  west-minus-east  (W-E)  difference.  Remarkably 
the  bifurcation  was  just  as  sharp  in  the  3D  model  as  in 
simple  ID  calculations  of  Holton  and  Dunkerton 
( 1978;  see  also  Yoden  1987a,b). 

Comparable,  though  less  dramatic,  results  were  ob¬ 
tained  when  the  seasonal  cycle  was  included  (without 
parameterized  wave  drag)  as  in  sections  2a.b.  Figure  2 
shows  the  time-averaged  December-February  (DJF) 
mean  zonal  flow  in  a  pair  of  runs  with  <j>„  =  375  m 
(Figs.  2a.b),  and  their  difference  (Fig.  2c).  The  aver¬ 
aging  interval  was  90  days,  beginning  on  day  80  of  the 
model  run  (December  1);  this  interval  will  apply 
throughout  the  paper  unless  otherwise  noted.  None  of 
the  experiments  revealed  any  significant  difference  be¬ 
tween  east  and  west  phase  in  early  winter  (prior  to 
December),  even  with  parameterized  wave  drag,  so 
there  was  little  point  in  including  this  time  period  in 
the  averaging  interval.  [The  model’s  early  winter  was 
defective,  since  observations  show  a  QBO  effect  in  No¬ 
vember  (Dunkerton  and  Baldwin  1991).]  As  in  per¬ 
petual  solstice  runs  near  the  bifurcation,  there  was  a 
major  warming  in  east  phase  (Fig.  2b),  but  not  west 
phase  (Fig.  2a),  resulting  in  a  large  DJF  W-E  differ¬ 
ence  (Fig.  2c).  The  difference  was  exaggerated  by  the 
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DJF  u  a  W  QBO 


DJF  G  C  W-E  QBO 


Fig.  2.  DJF  mean  zonal  wind  in  seasonal  integrations  with¬ 
out  parameterized  wave  drag;  forcing  amplitude  4>m  =*  375  m. 
(a)  West  phase;  (b)  east  phase;  (c)  difference  between  west 
and  east  phase. 


approach  to  radiatively  determined  conditions  in  the 
west  phase,  resulting  in  a  mean  flow  (Fig.  2a)  almost 
the  same  as  if  there  were  no  planetary  wave  forcing 
at  all. 

Figures  3a,b  display  the  seasonal  evolution  of  mean 
zonal  wind  at  75°N,  0.5  mb  in  west  and  east  phase, 
respectively,  in  several  runs  without  parameterized 
wave  drag.  These  time  series  can  be  compared  to  Yod- 
en’s  (1990,  his  Fig.  5)  and  are  qualitatively  similar. 
West-phase  experiments  were  radiatively  determined, 
except  those  with  forcings  in  excess  of  400  m.  Mean 
zonal  wind  near  the  pole  increased  slightly  with  in¬ 
creasing  forcing,  due  to  asymmetric  displacement  of 
the  vortex  off  the  pole  (which  should  not  obscure  the 
fact  that  integrated  angular  momentum  is  diminished 
by  the  vertical  EP  flux  entering  the  domain  at  the  lower 
boundary).  East-phase  experiments  were  perturbed  by 
forcings  exceeding  300  m;  some  of  these  contained 
strong  vacillations  after  the  initial  breakdown.  Al¬ 
though  sharp  bifurcations  were  not  apparent,  there  was 


a  maximum  sensitivity  of  DJF  mean  flow  to  forcing 
amplitude,  located  between  400  and  450  m  in  the  west 
phase  and  325  and  375  m  in  the  east  phase.  Forcings 
between  350  and  400  m  therefore  maximized  the  effect 
of  the  QBO  (cf.  Fig.  2c). 

The  QBO’s  influence  could  be  visualized  as  a  shift 
of  the  bifurcation  point,  which  in  this  case  is  no  longer 
a  “point”  (as  in  perpetual  solstice  runs)  but  rather  a 
range  of  <j>m  values  over  which  the  time-averaged  flow 
is  most  sensitive  to  changes  in  <f>m.  The  shift  is  toward 
greater  <f>m  in  the  west  phase,  compared  to  east-phase 
integrations  and  control  runs  with  no  QBO.  Alterna¬ 
tively,  the  QBO’s  effect  is  analogous  to  a  change  in 
<£m,  a  flip  from  west  to  east  phase  being  roughly  equiv¬ 
alent  to  a  75-100  m  increase  of  forcing. 

In  seasonal  integrations,  the  response  of  time-aver¬ 
aged  flow  to  Rossby  wave  forcing  was  therefore  a 
smooth  function  of  cj>m ,  unlike  perpetual  solstice  runs 
with  a  near-discontinuity  at  the  bifurcation  point.  The 
radiative  drive  and  lower  boundary  forcing  varied  as  a 
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function  of  time.  It  was  not  merely  a  question  of  warming  would  occur  and  (if  so)  how  large  it  would 
whether  a  major  warming  would  occur  for  a  particular  be.  In  seasonal  integrations  there  are  more  degrees  of 
<pm  (as  in  perpetual  solstice  runs),  but  when  such  a  freedom,  as  it  were,  that  determine  the  response  to  d>„. 
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Despite  the  added  complexity,  it  is  useful  to  imasjine  a 
"fuzzy  bifurcation”  in  the  response,  demarcating  the 
range  of  <bm  values  over  which  the  stratospheric  flow 
changes  from  an  essentially  radiative  state  ( at  small 
d>m)  to  a  highly  disturbed  state  ( at  large  6m). 

The  DJF  time  period  was  obviously  important  to  the 
seasonal  development  of  extratropical  QBO  in  the 
model  and  is  emphasized  for  the  following  reasons.  As 
already  noted,  differences  were  insignificant  in  early 
winter,  prior  to  December.  Nor  was  there  a  large  QBO 
effect  in  late  winter,  after  February.  The  polar  vortex 
breaks  down  as  spring  approaches:  differences  due  to 
the  QBO  are  not  expected  once  this  has  happened.  Late 
winter  warmings  are  a  consequence  of  diminished  ra¬ 
diative  forcing  and  increasing  susceptibility  to  wave 
transport  as  the  vortex  erodes  away;  these  breakdowns 
are  expected  regardless  of  the  phase  of  the  QBO  ( Hol¬ 
ton  and  Austin  1991 ).  Large  W—E  differences  are  ob¬ 
tained  when  a  midwinter  warming  occurs  in  one  phase 
but  not  in  the  other,  and  the  averaging  interval  is  chosen 
accordingly,  as  in  Fig.  2.  The  90-day  DJF  averaging 
interval  (days  80-170)  captured  the  W-E  difference 
well,  as  in  observations  (Dunkerton  and  Baldwin  1991, 
1992).  The  modeled  seasonal  development  was  fairly 
realistic  but.  as  already  noted,  its  early  winter  was  too 
quiescent  and  devoid  of  any  QBO  effect. 

Needless  to  say,  the  model’s  extratropical  QBO  was 
well  defined  in  midwinter.  In  the  real  atmosphere  the 
situation  is  not  so  clear  cut.  While  there  is  a  tendency 
for  major  midwinter  warmings  to  occur  preferentially 
in  the  east  phase,  this  is  hardly  significant,  since  warm¬ 
ings  also  occur  in  the  west  phase  and.  in  any  case,  the 
WMO  definition  of  a  “major”  warming  is  arbitrarv. 
Observations  indicate  a  significant,  but  much  smaller, 
composite  difference  than  shown  in  Fig.  2c. 

Although  the  observed  W— E  difference  may  be  due, 
in  part,  to  interannual  variations  unrelated  to  the  QBO 
(e.g.,  tropospheric  forcing),  this  explanation  is  inade¬ 
quate.  The  simulated  difference  (Fig.  2c)  came  from 
an  unrealistic  radiatively  determined  state  in  the  west 
phase  (Fig.  2a),  which  is  not  observed  in  either  hemi¬ 
sphere.  Obviously  important  are  transport  mechanisms, 
other  than  vertically  propagating  planetary  waves,  that 
restrict  the  amplitude  of  the  mesospheric  jet.  The  meso¬ 
sphere  contains  barotropic  instabilities  (Hartmann 
1983)  that  redistribute  angular  momentum  horizontally 
and  homogenize  the  high-latitude  flow  ( Dunkerton  and 
Delisi  1985a)  and  large-amplitude  gravity  waves 
(Fritts  1984)  that  transport  momentum  vertically  and 
are  thought  to  decelerate  the  flow  on  the  order  of  100 
ms"'  day-1  (Lindzen  1981). 

While  interesting  in  their  own  right,  these  transport 
mechanisms  can  also  affect  planetary  wave  propaga¬ 
tion.  Diminished  mean  flow  in  high  latitudes  alters  the 
vertical  propagation  of  planetary  waves  and  affects  the 
bifurcation  property  of  the  vortex.  A  weakened  vortex 
enables  vertical  Rossby  wave  propagation  (Chamey 


and  Drazin  1961  ).  increasing  the  probability  of  wave  - 
mean  flow  interaction  and  vortex  breakdown. 

b.  Experiments  with  wave  drag  parameterization 

Numerical  integrations  with  parameterized  wave 
drag  produced  a  more  realistic  polar  night  jet  and,  as  a 
result,  a  more  realistic  extratropical  QBO.  Figures  4a.b 
show  latitude -time  plots  of  mean  zonal  wind  ( west  and 
east  phase,  respectively)  at  8  mb.  with  planetary  wave 
forcing  <bm  =  275  m.  Both  experiments  were  “subcrit- 
ical’  ’  in  that  neither  led  to  a  major  midwinter  warmina. 
There  was  some  contraction  and  deceleration  of  the  jet 
after  the  midpoint  of  integration,  the  details  differing 
slightly  between  the  two  experiments.  At  this  level, 
simulations  did  not  diverge  significantly  until  after  1 60 
days,  when  a  breakdown  occurred  in  the  east  phase. 
Similar  evolution  was  observed  at  <bm  =  250  m  ( not 
shown). 

With  slightly  larger  forcing  ( d>m  =  300  m)  there  was 
a  major  midwinter  wanning  in  the  east  ( but  not  west ) 
phase  (Figs.  5a,b).  As  shown  later,  this  pair  of  inte¬ 
grations  gave  the  maximum  W-E  difference  in  the 
stratosphere,  the  sign,  shape,  and  magnitude  of  ±u  be¬ 
ing  consistent  with  observations.  There  was  an  overall 
reduction  of  DJF  mean  flow  north  of  about  45°N  and 
slight  acceleration  to  the  south  (note  formation  of  a 
weak  jet  near  30°N  between  days  120  and  140  in  Fig. 
5b).  Similar  evolution  was  observed  at  <j>m  =  325  m 
( not  shown )  except  that  a  minor  midwinter  breakdown 
also  occurred  in  the  west  phase. 

At  still  higher  forcing  (</>„  =  350  m),  midwinter 
warmings  occurred  in  both  QBO  phases  (Figs.  6a.b), 
although  the  patterns  were  not  identical.  The  east  phase 
event  qualified  as  a  “major”  sudden  warming  by  day 
100,  followed  by  some  strong  vacillations,  while  the 
west  phase  evolution  (after  an  initial  deceleration  at 
day  85)  resembled  the  gradual  deceleration  and  break¬ 
down  seen  in  previous  examples.  Apparently  the  west 
phase  transition  or  “fuzzy  bifurcation”  lay  between 
300  and  350  m.  whereas  the  east  phase  transition  was 
between  275  and  300  m.  We  emphasize,  however,  that 
these  transitions  (as  a  function  of  4>„ ),  like  Fig.  3,  were 
nowhere  as  abrupt  as  those  observed  in  perpetual  sol¬ 
stice  runs  without  parameterized  wave  drag. 

DJF  W-E  mean  zonal  wind  differences  for  several 
experiments  ( <b„  =  275,  300.  325,  350  m)  are  sum¬ 
marized  in  Figs.  7a-d.  The  extratropical  QBO  began 
at  low  forcing  with  a  positive  anomaly  in  the  polar 
mesosphere  (Fig.  7a)  that  descended,  with  increasing 
forcing,  to  the  stratosphere  (Fig.  7b).  Simultaneously 
a  weak  negative  anomaly  appeared  in  the  midlatitude 
stratosphere,  arching  into  the  mesosphere  on  top  of  the 
positive  anomaly  (Fig.  7c).  At  large  forcing  both  pos¬ 
itive  and  negative  components  weakened  to  insignifi¬ 
cant  values  (Fig.  7d). 

This  sequence  can  be  contrasted  to  O’Sullivan  and 
Young  ( 1992),  where  the  basic  QBO  dipole— a  pos- 
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TIME  (days) 

Fig.  4.  Latitude -time  cross  sections  of  mean  zonal  wind  at  8  mb  in  seasonal  integrations  with  parameterized  wave 
drag;  forcing  amplitude  <b„  =  275  m.  (a)  West  phase:  (b>  east  phase.  Model  output  smoothed  with  3-day  running  mean. 


itive  (negative)  anomaly  in  polar  (middle)  latitudes — 
appeared  in  most  of  their  experiments  with  little  ver¬ 
tical  structure,  maximizing  in  the  mesosphere.  In  Fig. 
7,  the  extratropical  anomaly  (when  significant)  maxi¬ 
mized  in  the  middle  or  upper  stratosphere  and  was 
overlain  by  a  mesospheric  anomaly  of  opposite  sign 
contiguous  to  the  midlatitude  half  of  the  stratospheric 
dipole.  The  mesospheric  response  may  be  a  shielding 
effect,  whereby  the  planetary  waves,  having  caused  the 
stratospheric  vortex  to  weaken  or  break  down,  no 
longer  penetrate  the  mesosphere  well  enough  to  over¬ 
come  the  restoring  force  of  radiative  relaxation.  This 
interpretation  is  supported  by  Figs.  8a.b,  showing  the 
latitude-time  structure  of  mean  zonal  wind  at  0.2  mb 
with  <bm  =  325  m.  The  corresponding  plots  of  wave- 
number  1  Rossbv  wave  amplitude  at  1 .0  mb  are  shown 
in  Figs.  9a.b;  Fig.  9c  shows  the  DJF-average  difference 
of  wave  amplitude.  At  upper  levels  it  was  the  west 
phase  that  was  more  exposed  to  Rossbv  waves,  causing 
a  more  sustained  disruption  of  the  flow  than  in  the  east 
phase.  Relaxation  of  upper-level  flow  is  vaguely  anal¬ 


ogous  to  the  shielding  effect  observed  in  the  equatorial 
QBO  (Plumb  1977). 

As  in  the  previous  subsection,  the  effect  of  the  equa¬ 
torial  QBO  on  the  extratropical  stratosphere  could  be 
likened  to  a  shift  of  the  point  of  maximum  sensitivity 
to  <bm ,  or  “fuzzy  bifurcation.”  Figure  10  illustrates  the 
effect  schematically  for  both  series,  showing  the  time- 
averaged  mean  zonal  flow  at  70°N.  8  mb.  as  a  function 
of  <bm.  The  largest  W-E  difference  was  obtained  in  the 
middle  range  where  the  curves  are  farthest  apart. 

Due  to  parameterized  wave  drag,  the  mesospheric 
response  was  more  complicated.  Figure  1 1  shows  the 
time-averaged  profiles  of  mean  zonal  wind  at  70°N 
from  the  experiments  of  Fig.  8.  Mean  zonal  winds  were 
decelerated  to  a  greater  degree  in  the  west  phase,  re¬ 
versing  the  sign  of  W-E  difference  compared  to  the 
stratospheric  extratropical  QBO.  Figure  1 1  also  shows 
analogous  profiles  obtained  from  Fig.  2a.b  without 
wave  drag.  There  was  no  reversal  in  this  case. 

Two  distinct  effects  of  the  wave  drag  parameteriza¬ 
tion  are  therefore  apparent  when  compared  to  results 
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of  the  previous  subsection.  First,  drag  reduced  the  am¬ 
plitude  of  the  mesospheric  jet  to  begin  with,  so  any 
deceleration  caused  by  Rossby  waves  was  necessarily 
less  than  that  of  Fig.  2c.  Second,  by  reducing  the  am¬ 
plitude  of  the  jet,  the  specified  drag  allowed  planetary 
Rossby  waves  to  propagate  to  the  mesosphere  regard¬ 
less  of  <!>m  and  enabled  smaller  values  of  <f>m  to  cause 
vortex  breakdown.  With  parameterized  wave  drag  there 
is  no  sharp  bifurcation  between  extremes,  because  the 
radiatively  determined  state  is  no  longer  admissible  as 
a  solution.  This  is  especially  relevant  in  the  meso¬ 
sphere,  where  planetary  Rossby  waves  cannot  penetrate 
the  radiatively  determined  state  (Chamey  and  Drazin 
1961;  cf.  Fig.  2a)  but  easily  propagate  into  a  jet  decel¬ 
erated  by  unresolved  body  forces3  (Fig.  8a)  if  they  are 
not  first  absorbed  in  the  stratosphere. 


3  Gravity  wave  drag  presumably  damps  planetary  Rossby  waves 
in  the  mesosphere  (Miyahara  1985).  This  effect  was  ignored,  but  if 
included  would  not  alter  our  conclusion  that  gravity  wave  drag,  first 


Clear  evidence  of  a  reversal  of  the  Holton -Tan  os¬ 
cillation  in  the  mesosphere  has  not  been  obtained  for 
lack  of  observations,  although  Dunkerton  and  Baldwin 
(1991)  found  a  maximum  extratropical  zonal-wind 
QBO  at  5  mb  ( with  diminished  amplitude  at  2  and  l 
mb).  Above  5  mb  the  polar  temperature  QBO  reversed 
sign  as  required  by  thermal  wind  balance.  Within  in¬ 
dividual  years,  observations  of  mesospheric  cooling 
above  stratospheric  warmings  are  fairly  common  ( Qui¬ 
roz  1969;  Labitzke  1972;  Hirota  and  Barnett  1977; 
Matsuno  and  Nakamura  1979;  Holton  1983;  Delisi  and 
Dunkerton  1988). 

The  sequence  of  extratropical  QBOs  in  Fig.  7  was 
not  unique  to  these  experiments.  A  better-resolved  se¬ 
ries  with  forcing  increments  as  low  as  12.5  m  (using  a 
different  frictional  drag  coefficient )  also  gave  a  contin¬ 
uous  variation  of  Aw  as  a  function  of  <f>„.  It  is  note¬ 
worthy  that  the  extratropical  QBO — when  defined  over 


of  all,  enables  planetary  wave  propagation  by  preventing  a  radiatively 
determined  state. 
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TIME  (days) 

Fig.  6.  As  in  Fig.  4  but  with  <t>m  =  350  m.  (a)  West  phase;  (b)  east  phase. 


a  suitably  long  time  interval,  that  is,  90  days— depends 
smoothly  on  <j>m  when,  in  fact,  the  day-to-day  evolution 
is  chaotic  after  the  initial  breakdown.  This  is  probably 
attributable  to  the  averaging  interval  being  long  com¬ 
pared  to  the  radiative  and  mechanical  damping  times. 

Little  interannual  variability  is  expected  in  a  strongly 
damped,  nonlinear  dynamical  system  with  fixed  exter¬ 
nal  forcing  (Yoden  1990) — a  conjecture  that  remains 
to  be  verified  in  our  model.  The  point  is  academic  for 
the  terrestrial  atmosphere  since  tropospheric  forcings 
vary  from  year  to  year.  Composites  of  the  extratropical 
QBO  blend  together  these  interannual  variations  un¬ 
related  to  the  QBO  and  are  therefore  not  expected  to 
match  precisely  the  results  of  any  individual  experi¬ 
ment  as  shown  in  Fig.  7,  for  example.  To  account  for 
interannual  (non-QBO)  variability  we  constructed  a 
“model  composite”  of  the  extratropical  QBO  averag¬ 
ing  together  results  of  the  four  experiments  shown  in 
that  figure,  which  probably  represent  a  typical  inter¬ 
annual  range  of  tropospheric  forcings.  Composites  of 
W-E  differences  in  zonal-mean  wind,  temperature, 
and  wave  amplitude  are  shown  in  Figs.  12a-c.  The 


wind  and  temperature  differences  agreed  reasonably 
well  with  observations,  both  in  magnitude  and  (for 
temperature)  the  altitude  of  sign  reversal. 

4.  Discussion 

The  sensitivity  of  extratropical  QBO  to  the  ampli¬ 
tude  of  forced  planetary  Rossby  waves  suggests  that 
the  Holton-Tan  oscillation  will  vary  from  year  to  year, 
as  tropospheric  forcing  varies.  Observations  show  that 
the  oscillation  is  not  constant,  but  varies  on  a  decadal 
timescale  (Dunkerton  and  Baldwin  1992). 

Generally  speaking,  factors  that  influence  the  extra¬ 
tropical  QBO  can  be  categorized  in  one  of  four  ways: 
effects  from  above  (mesosphere),  from  low  latitudes 
(equatorial  stratosphere),  in  situ  anomalies  (polar 
stratosphere),  and  forcing  from  below  (troposphere). 

a.  Mesosphere 

The  mesosphere  can  affect  planetary  wave  propa¬ 
gation  and  alter  the  seasonal  development  of  strato- 
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Fig.  7.  DJF  mean  zonal  wind  difference  (west  minus  east  phase)  for  various  forcing  amplitudes  (f>m; 
(a)  275  m;  (b)  300  m;  (c)  325  m:  (d)  350  m. 


spheric  circulation.  Gravity  wave  drag  illustrates  how 
unresolved  body  forces  in  the  mesosphere  indirectly 
affect  the  stratosphere.  By  retarding  the  mean  flow, 
gravity  waves  establish  a  broad  waveguide  for  plane¬ 
tary  Rossby  waves  that  would  otherwise  be  absent.  To 
be  sure,  planetary  waves  decelerate  the  mean  flow  and 
thereby  maintain  their  own  waveguide.  This  is  true 
even  in  experiments  without  parameterized  wave  drag. 
However,  the  extreme  variation  between  radiatively  de¬ 
termined  and  disturbed  states  (section  3a)  suggests  that 
gravity  wave  drag,  or  some  other  retarding  force,  is 
necessary  for  a  realistic  seasonal  cycle  and  extratropi- 
cal  QBO. 

It  is  conceivable  that  other  dynamical  processes 
(e.g.,  barotropic  instability)  and  thermodynamic 
ones  (e.g.,  anomalous  heating  due  to  the  solar  cycle) 
could  have  some  downward  influence  on  the  strato¬ 


sphere.  Interest  in  the  solar  cycle  arises  from  an  ap¬ 
parent  correlation  with  the  amplitude  of  extratropi- 
cal  QBO  ( Labitzke  and  van  Loon  1988:  Dunkerton 
and  Baldwin  1992;  Kodera  1993).  Observations 
suggest  that  the  December  mesospheric  jet  is  anom¬ 
alously  strong  near  solar  maximum  (Kodera  and 
Yamazaki  1990).  (Whether  this  association  will 
survive  the  test  of  time  remains  to  be  seen.  Centuries 
will  pass  before  solar  cycle  effects  can  be  estab¬ 
lished  with  statistical  confidence.)  Kodera  et  al. 
( 1991 )  attempted  to  force  the  observed  wind  vari¬ 
ation  in  a  numerical  model  by  altering  the  solar 
heating  due  to  ozone.  The  magnitude  of  their  im¬ 
posed  change  was  significantly  larger  than  current 
estimates  (by  a  factor  of  10—20).  The  pattern  of 
December  wind  anomaly  in  Kodera  and  Yamazaki 
( 1 990 ) — fairly  deep  but  restricted  to  a  narrow  band 
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of  latitudes — suggests  an  anomalous  body  force 
(not  anomalous  heating)  as  the  source.4 

b .  Equatorial  stratosphere 

The  timing  of  the  QBO  will  vary  from  year  to  year 
since  it  is  not  synchronized  with  the  seasonal  cycle. 
Over  several  winters  the  high-latitude  flow  will  “see” 
different  profiles  of  tropical  mean  wind.  This  will  mod¬ 
ulate  the  extratropical  QBO  on  a  timescale  from  5  to 
11  years,  corresponding  to  the  length  of  time  required 
for  QBO  onsets  to  advance  12  months  through  the  cal¬ 
endar  year  (Gray  and  Dunkerton  1990).  Further  inves- 


4  Our  results,  to  be  reported  elsewhere,  demonstrate  that  it  is  pos¬ 
sible  to  modify  the  extratropical  QBO  if  a  wind  variation  like  that  of 
Kodera  and  Yamazaki  is  imposed  on  the  model.  Thermal  forcing 
alone  is  insufficient  since  the  magnitude  of  variation  is  unreasonably 
large,  and  thermal  forcing  is  accompanied  by  nonlocal  effects  in  lat¬ 
itude  such  as  a  stronger  diabatic  circulation  and  easterlies  (hence, 
relative  cooling)  near  the  tropical  stratopause  (Dunkerton  1991). 


tigation  of  modulation  will  require  multiyear  exten¬ 
sions  of  the  seasonal  integrations  performed  here. 

Anomalous  heating  due  to  volcanic  aerosol  in  the 
tropical  lower  stratosphere  was  shown  by  Pitari  ( 1993 ) 
to  enhance  planetary  wave  propagation  in  a  numerical 
model. 

c.  Polar  lower  stratosphere 

Anomalous  radiative  effects  in  the  polar  lower 
stratosphere  (e.g.,  due  to  ozone  depletion)  can  have  a 
positive  feedback  on  planetary  wave  propagation 
(Houben  1989,  unpublished  abstract).  Unlike  meso¬ 
spheric  anomalies,  perturbations  in  the  lower  strato¬ 
sphere  alter  the  mean  wind  at  the  point  of  entry  for 
planetary  waves  and  can  therefore  have  a  significant 
effect  even  when  the  imposed  changes  are  small. 

d.  Troposphere 

Tropospheric  variations  affect  the  stratosphere,  not 
only  via  forced  waves  but  as  the  whole  “lower  bound- 
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Fig.  9.  As  in  figs.  8a,b  but  for  geopotential  height  ampli¬ 
tude  of  zonal  wavenumber  1  at  1.0  mb.  (c)  DJF  W-E  differ¬ 
ence  of  wave  amplitude. 


ary  condition”  (including  the  mean  zonal  wind)  is 
modified  by  low-frequency  variability.  Certain  patterns 
of  low-frequency  variability  in  the  troposphere  are  fa¬ 
vorably  associated  with  planetary  wave  excitation  and 
vortex  breakdown  (Baldwin  et  ai.  1994). 


If  results  of  Holton  and  Dunkerton  (1978)  are 
valid  in  a  three-dimensional  context,  bifurcation  of 
the  polar  vortex  depends  also  on  the  specified  mean 
zonal  wind  at  the  tropopause.  An  investigation  of  the 
stratospheric  response  to  tropospheric  low-frequency 
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Fig.  10.  Values  of  DJF  mean  zonal  wind  at  70°N,  8  mb  in  several 
experiments  with  parameterized  wave  drag,  illustrating  the  shift  of 
“fuzzy  bifurcation”  from  east  to  west  phase. 


variability  is  in  progress,  and  results  will  be  reported 
elsewhere. 

5.  Conclusions 

A  three-dimensional  nonlinear  numerical  model  of 
the  middle  atmosphere  was  used  to  investigate  the  sea- 


OJF  MEAN  ZONAL  WIND 


Fig.  11.  Vertical  profiles  of  DJF  mean  zonal  wind  in  experiments 
with  parameterized  wave  drag  {<f>m  =  325  m),  and  without  parame¬ 
terized  wave  drag  (<f>m  =  375  m.) 


DJF  G  a  W-E  QBO 


DJF  T  b  W-E  QBO 


DJF  wove  1  C  W-E  QBO 


Fig.  12.  Model  composite  of  DJF  W-E  difference  combining  re¬ 
sults  from  the  four  experiments  of  Fig.  7.  (a)  Mean  zonal  wind:  (b) 
temperature:  (c)  geopotential  height  amplitude  of  zonal  wavenumber 
i. 
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sonal  development  of  the  extratropieai  quasi-biennial 
oscillation  ( QBO  > .  Several  twin  experiments  were  per¬ 
formed.  differing  initially  only  in  the  phase  of  equa¬ 
torial  QBO.  with  identical  planetary  Rossby  waves 
forced  at  a  tropopause  lower  boundary.  Integrations 
evolved  during  winter  in  a  manner  roughly  consistent 
with  observations,  but  were  sensitive  to  the  lower 
boundary  forcing  and  the  parameterization  of  radiative 
cooling  and  unresolved  body  force.  Stronger  radiative 
forcing  (compared  to  O’Sullivan  and  Young  1992) 
heightened  the  sensitivity  of  the  polar  vortex  to  tropo¬ 
spheric  forcing,  resulting  in  a  large  extratropieai  QBO 
when  the  wave  forcing  6,„  was  near  a  '‘bifurcation 
point.'*  As  in  Holton  and  Dunkerton  ( 1978),  the  bi¬ 
furcation  is  the  value  of  <f>,„  for  which  the  time-averaged 
flow  is  most  sensitive  to  changes  in  <f>m.  The  effect  was 
most  dramatic  in  perpetual  solstice  runs  where,  for  cer¬ 
tain  values  of  <bm.  the  equatorial  QBO  caused  the  model 
to  cross  from  a  radiativeiy  determined  state  to  a  major 
warming.  The  sensitivity  was  reduced  but  not  elimi¬ 
nated  in  experiments  with  a  seasonal  cycle,  so  it  is  ap¬ 
propriate  to  speak  of  a  "fuzzy  bifurcation”  in  this  case. 

The  unrealistically  large  amplitude  of  extratropieai 
QBO  was  reduced  substantially  by  parameterized  wave 
drag,  which  had  two  distinct  effects.  First,  maximum 
wind  speeds  in  the  polar  vortex  were  reduced  so  that 
any  composite  difference  would  be  considerably  less 
than  the  difference  between  radiativeiy  determined  and 
disturbed  states  alluded  to  in  the  previous  paragraph. 
Second,  propagation  of  Rossby  waves  was  enhanced 
(Chamey  and  Drazin  1961 ),  making  it  easier  for  the 
polar  vortex  to  break  down.  Equivalently,  less  tropo¬ 
spheric  forcing  was  required  to  achieve  a  breakdown. 
This  is  important  in  the  mesosphere  since  the  radia¬ 
tiveiy  determined  state,  unhindered  by  wave  drag,  pre¬ 
cludes  vertical  propagation  of  planetary  waves. 

Realistic  extratropieai  QBOs  were  obtained  for  ob¬ 
served  values  of  <bm  when  excessive  westerlies  were 
relaxed  to  climatological  values  ( Barnett  and  Comey 
1985)  using  a  simple  parameterization  scheme.  Unlike 
previous  model  results,  but  consistent  with  observa¬ 
tions  (Dunkerton  and  Baldwin  1991 ),  the  extratropieai 
QBO  maximized  in  the  upper  stratosphere  when  real¬ 
istic  <j>m  were  specified.  The  anomaly  reversed  sign  in 
the  mesosphere,  apparently  due  to  a  shielding  effect  as 
explained  in  section  3b.  Further  observations  are 
needed  to  determine  whether  the  phase  reversal  is  sig¬ 
nificant. 

Numerical  experiments  demonstrate  that  the  extra- 
topical  QBO  is  sensitive  to  the  amplitude  of  planetary 
Rossby  waves  entering  the  stratosphere,  and  requires 
additional  wave  drag  in  the  mesosphere  to  be  realistic. 
Other  factors  may  also  affect  the  oscillation,  such  as 
anomalous  heating,  timing  of  equatorial  QBO,  ozone 
depletion,  and  tropospheric  low-frequency  variability. 
Further  investigation  of  some  of  these  mechanisms  may 
be  necessary  if  we  are  to  understand  the  observed  vari¬ 
ation  of  extratropieai  QBO. 


Finally,  our  experiments  suggest  a  paradox  in  trop¬ 
ical -extratropieai  interaction:  how  the  equatorial 
QBO,  on  the  one  hand,  could  be  shielded  from  the  in¬ 
fluence  of  midlatitude  planetary  waves  while,  on  the 
other  hand,  having  a  significant  influence  on  these 
waves  in  extratropieai  latitudes.  (The  simulated  west 
phase  of  equatorial  QBO  would  be  eliminated  by  plan¬ 
etary  wave  transport  in  many  cases,  apart  from  a  west¬ 
erly  restoring  force.)  Would  a  multiyear  numerical  in¬ 
tegration,  containing  Kelvin  waves,  Rossby -gravity 
waves,  and  midlatitude  planetary  waves  equilibrate  to 
a  realistic  equatorial  QBO,  generate  an  extratropieai 
signal,  but  remain  unaffected  at  low  latitudes  by  mid¬ 
latitude  waves?  We  speculate  that  this  may  require  for¬ 
mation  of  a  region  of  weak  climatological  easterlies  in 
the  subtropics  ( Dunkerton  and  Delisi  1985b).  Experi¬ 
ments  with  more  realistic  QBO  structure  and  evolution 
would  shed  light  on  the  momentum  balance  of  the  QBO 
(as  affected  by  lateral  propagation  of  planetary  waves, 
particularly  in  the  subtropics),  the  nature  of  the  sub¬ 
tropical  “dynamical  barrier,”  and  its  effect  on  wave 
propagation  and  transport.  They  would  also  demon¬ 
strate  whether  the  extratropieai  effects  described  in  this 
paper  can  occur  under  more  realistic  conditions. 

We  anticipate  that  it  may  soon  be  feasible  to  perform 
a  few  such  integrations  on  supercomputers. 
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Abstract 

The  effect  of  conditional  heating  on  unstable  Kelvin  waves  generated  by  CISK  and  evaporation- 
wind  feedback  was  investigated  with  multi-level  models  on  the  equatorial  beta-plane.  In  agreement 
with  previous  analytic  results,  the  scale-selection  catastrophe  of  linear  heating  is  modified  but  not 
eliminated  by  conditional  heating.  There  is  only  one  region  of  upward  motion,  but  it  contracts  to 
the  smallest  zonal  scale  allowed  by  diffusion  or  model  resolution.  The  pathology  is  most  severe  in 
the  CISK-unstable  case.  With  evaporation-wind  feedback,  Kelvin-wave  instability  is  possible  apart 
from  CISK;  growth  rate  in  this  case  increases  slowly  with  decreasing  diffusion,  but  contraction  is  still 
observed.  We  note  incidentally  that  eastward  moving  evaporation-wind  modes  can  be  generated  on 
the  beta-plane  with  no  mean  wind,  contrary  to  previous  assumptions  requiring  an  easterly  background 
flow. 

While  evaporation-wind  feedback  is  an  improvement  over  CISK,  as  far  as  scale  selection  is  concerned, 
zonal  phase  speeds  remain  too  large  compared  to  the  observed  propagation  of  tropical  intraseasonal 
oscillations  and  superciusters,  unless  (1)  the  parameter  q  controlling  the  strength  of  moisture  conver¬ 
gence  feedback  lies  close  to  the  CISK-neutral  point,  (2)  the  parameter  A  controlling  the  strength  of 
evaporation-wind  feedback  is  small,  and  (3)  the  heating  profile  is  assumed  proportional  to  vertical 
velocity.  Small  positive  phase  speeds  are  obtained  under  these  conditions,  for  example,  in  a  two-level 
model. 


1.  Introduction 

Theoretical  models  of  the  tropical  intraseasonal 
oscillation  (TIO:  Madden  and  Julian,  1971,  1972: 
Knutson  and  Weickmann,  1987;  Rui  and  Wang, 
1990)  have  often  relied  on  CISK  (Lau  and  Peng, 
1987;  Chang  and  Lim,  1988)  or  evaporation- wind 
feedback  (Emanuel,  1987;  Neelin  et  a/.,  1987)  to 
generate  unstable  waves.  Similar  mechanisms  were 
shown  to  be  important  in  GCMs  containing  in¬ 
traseasonal  oscillations  (Hayashi  and  Sumi,  1986: 
Swinbank  et  a/.,  1988;  Lau  et  aL  1988).  Never¬ 
theless,  these  ‘explanations5  of  TIO  are  incomplete, 
as  they  fail  to  reproduce  the  slow  eastward  phase 
speed,  horizontal  structure,  and  observed  hierar¬ 
chy  of  convection  within  intraseasonal  oscillations. 
Furthermore,  the  exponentially-growing  modes  are 
pathological,  preferring  the  smallest  zonal  scale. 
This  is  true  whether  linear  (positive  and  negative) 
or  conditional  (positive-only)  heating  is  assumed. 

We  recently  documented  the  scale-selection  catas¬ 
trophe  analytically  with  conditional  heating  (Crum 
and  Dunkerton.  1992;  hereafter  CD)  and  now  ex¬ 
tend  the  demonstration  to  multi-level  models  on 
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an  equatorial  beta-plane  (Section  2b).  Our  conclu¬ 
sion  concerning  CISK  modes  is  in  agreement  with 
a  statement  by  Yoshizaki  (1991)  who,  apart  from 
this,  did  not  demonstrate  the  catastrophe  nor  inves¬ 
tigate  the  role  of  evaporation-wind  feedback.  Be¬ 
cause  linear  theory  predicts  a  weaker  catastrophe 
with  evaporation-wind  feedback,  we  were  motivated 
to  explore  the  effect  of  conditional  heating  and  beta- 
plane  geometry  (Section  2c).  In  this  context  we 
demonstrate  that  it  is  possible  to  avoid  the  lineariza¬ 
tion  of  evaporation  and  still  obtain  an  eastward- 
propagating  instability.  Variation  of  zonal  phase 
speed  as  a  function  of  moisture  convergence  and 
evaporation-wind  feedback  is  briefly  discussed  in 
Section  2d. 

2.  Model  results 

a.  Numerical  model 

We  consider  linear,  hydrostatic  perturbations 
from  a  time- independent  basic  state  on  an  equatorial 
beta-plane.  A  three-dimensional  primitive-equation 
model  was  formulated  in  pressure  coordinates  with 
two  or  five  vertical  levels  (Charney-Phillips  grid), 
finite  differences  in  pressure  and  latitude,  and  a 
Fourier  spectral  representation  in  longitude.  Solu- 
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tions  were  periodic  in  x,  dependent  variables  were 
set  to  zero  at  the  sidewalls  (±60°  latitude),  and  ver¬ 
tical  boundary  conditions  were  uj  =  0  at  p  =  0 
and  p  =  poo  =  1000  mb.  Static  stability  was  derived 
from  Houghton  (1986)  at  10°N.  A  small  second- 
order  horizontal  diffusion  was  inserted  on  the  rhs 
of  horizontal  momentum  and  temperature  equations 
(which  nevertheless  proved  to  be  essential  for  well- 
resolved  simulations).  Precipitation  was  parameter¬ 
ized  in  terms  of  total  moisture  entering  a  column, 
including  that  from  surface  evaporation: 

rPoo 

I  V  ‘{Vq)d'pfg  +  E0o,  if  >  0: 
o 

0,  otherwise. 


(2.1) 


where  P  is  precipitation  in  units  of  kgm~2  s"1,  q  is 
water  vapor  mixing  ratio,  V  is  horizontal  velocity, 
and  £oo  is  surface  evaporation  is  units  of  P.  Spa¬ 
tially  constant  q  was  assumed,  non-zero  only  in  the 
lowest  model  layer,  and  heating  was  distributed  ver¬ 
tically  by  a  specified  profile  (CD).  In  the  five-level 
model,  our  heating  profile  rj  was  defined  as  (moo, 
moo,  Veoo,  Vsoo)  (1/2,  1/2,  2,  1)  where  subscripts 
refer  to  pressure  level.  This  profile  generates  the 
‘slow  CISK’  mode  of  Chang  and  Lim  (1988). 

Readers  familiar  with  CISK  will  recall  that  the 
specified  heating  profile  is  the  mathematical  device 
responsible  for  vertical  mode-coupling;  in  unstable 
CISK,  two  stable  modes  have  coalesced  to  form  an 
unstable  mode  (CD).  It  should  be  noted  that  by 
specifying  the  heating  profile  a  priori ,  instead  of 
making  it  proportional  to  vertical  velocity,  zonal 
propagation  of  CISK  modes  can  occur  (as  explained 
in  Appendix  A  of  CD).  The  second  alternative  im¬ 
plies  nonpropagating  convection  unless  additional 
sources  of  heat  are  present,  such  as  might  be  due 
to  evaporation-wind  feedback. 

We  limit  our  application  of  the  five-level  model 
primarily  to  CISK  modes  obtained  with  a  speci¬ 
fied  heating  profile.  To  better  approximate  the  con¬ 
cept  of  reduced  moist  static  stability  in  evaporation- 
wind  feedback  (Yano  and  Emanuel,  1991)  a  two-level 
model  was  used  following  Neelin  et  al .  (1987).  In 
this  case  there  is  no  heating  ‘profile'  as  such;  instead, 
77  cc  1  at  a  single  internal  half-level,  tantamount  to 
a  reduced  moist  static  stability  (or  convective  insta¬ 
bility  in  unstable  CISK).  From  this  discussion  it  can 
be  seen  that  there  is  an  important  conceptual  and 
practical  distinction  to  be  made  between  (1)  the  use 
of  a  specified  heating  profile  proportional  to  verti¬ 
cal  velocity  at  cloud  base  and  (2)  the  assumption  of 
heating  proportional  to  vertical  velocity  at  all  levels. 
We  return  to  this  in  Section  2d  when  discussing  the 
zonal  phase  speed  of  unstable  modes. 

In  what  follows,  ‘conditional'  heating  refers  to  the 
nonlinearity  introduced  by  Eq.  (2.1).  ‘Linear'  heat¬ 
ing  employs  the  upper  formula  of  Eq.  (2.1)  irrespec¬ 
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tive  of  its  sign.  In  both  cases,  the  problem  is  dynam¬ 
ically  linear,  independent  of  solution  amplitude,  and 
exponentially-growing  normal  modes  are  possible. 

6.  Wave-CISK 

In  this  case,  evaporation  is  ignored  in  Eq.  (2.1). 
With  our  specified  heating  profile,  the  third  and 
fourth  internal  modes  coalesce  near  $=4gkg~l  to 
form  an  unstable  propagating  mode  with  phase 
speed  15-20  ms-1.  We  found  similar  behavior  in 
models  with  many  more  vertical  levels.  Due  to  coa¬ 
lescence,  vertical  structure  approximates  a  'modified 
first-internal  mode'  with  no  internal  zero-crossings 
of  vertical  velocity  u.  This  property  is  desirable  ob- 
servationally.  although  the  'slow-CISK'  heating  pro¬ 
file  generates  outflow  too  low  (near  400  mb),  while 
phase  speed  remains  too  large.  Taller  heating  pro¬ 
files  cause  faster  propagation.  Our  immediate  intent 
is  not  to  dwell  on  these  defects  but  to  demonstrate 
the  importance  of  spatial  filtering  in  zonal  scale  se¬ 
lection.  With  linear  heating,  analytic  Kelvin  wave 
solutions  are  nondispersive  so  that  growth  rate  is 
proportional  to  zonal  wavenumber  times  the  imagi¬ 
nary  part  of  the  eigenvalue  (complex  phase  speed) — 
as  is  generally  true  in  the  long-wave  approximation. 
With  conditional  heating,  the  scale-selection  catas¬ 
trophe  is  modified  but  not  eliminated.  Only  one 
region  of  upwelling  is  preferred,  but  fastest  growth 
occurs  at  the  smallest  zonal  scale.  CD  demonstrated 
this  pathology  analytically  using  simple  two-  and 
three-level  models  in  the  equatorial  x-z  plane. 

The  same  is  evidently  true  on  an  equatorial  beta- 
plane,  when  the  number  of  vertical  levels  is  in¬ 
creased.  Figure  1  shows  the  growth  rate  of  stable 
CISK  modes  with  conditional  heating,  q  =  8  gkg~l, 
obtained  numerically  with  the  five-level  beta-plane 
model  described  in  Section  2a.  The  initial  condi¬ 
tion  was  a  Gaussian  temperature  profile  in  latitude 
and  longitude,  maximum  in  mid-troposphere.  Var¬ 
ious  curves  refer  to  different  model  resolutions  (M, 
N)  where  M  is  the  number  of  zonal  harmonics  and 
N  is  the  number  of  latitudinal  grid  points  in  one 
hemisphere.  The  dotted  line  shows  the  width  of  the 
region  of  upward  motion  for  (Af,  N)  =  (128, 44).  Re¬ 
sults  of  Fig.  1  are  exactly  reminiscent  of  CD’s  Fig.  6, 
showing  that  only  a  few  of  the  simulations  are  ‘well- 
resolved’  (independent  of  resolution).  Coarse  res¬ 
olution  severely  underestimates  growth  rate  at  low 
diffusion.  In  well-resolved  simulations,  wet-region 
width  decreases  catastrophically  with  diffusion  as  in 
the  analytic  solutions  of  CD. 

Horizontal  structure  of  these  solutions  is  essen¬ 
tially  a  Kelvin  wave  with  negligible  meridional  ve¬ 
locity.  The  vertical  structure  is  a  modified  first  in¬ 
ternal  mode  similar  to  Lau  and  Peng  (1987). 

As  an  aside,  we  note  that  gravity  waves  do  not 
dominate  the  simulated  modes.  In  linear  CISK, 
gravity  waves  grow  slightly  faster  than  Kelvin  waves 
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Fig.  L  CISK  growth  rates  in  the  five-level 
model  as  a  function  of  diffusion  for  q  =  8 
gkg"*1.  Solid  circles,  open  circles,  and 
squares  denote  simulation  with  (M,  N)  = 

(32, 12),  (64,22),  and  (128,44)  respectively. 

The  e-folding  width  of  ascending  motion 
is  shown  by  the  dashed  line  for  (Af,  N)  = 
(128,44). 

(Crum  and  Stevens,  1983).  A  series  of  tests  with¬ 
out  the  ^/-component  of  diffusion  demonstrated  that 
it  was  the  conditional  heating  and  x-component  of 
diffusion,  and  not  the  damping  of  oscillatory  struc¬ 
ture  in  latitude,  that  suppressed  gravity  waves.  The 
viscous  beta-plane  result  is  interesting  because,  even 
with  small  diffusion,  conditional  heating  leads  to  a 
narrow,  propagating  convergence  zone  perpendicu¬ 
lar  to  the  equator  rather  than  a  circularly-symmetric 
singularity  as  assumed  in  the  /-plane  solutions  of 
Charney  and  Eliassen  (1964). 

c.  Evaporation-wind  feedback 
In  this  case  evaporation  in  Eq.  (2.1)  is  parameter¬ 
ized  according  to  the  standard  bulk  formula 

Eq0  =  pscD  (q*  -  q)  \V3\  (2.2) 

where  ps  is  surface  air  density,  cq  a  dimensionless 
coefficient,  q*  is  saturation  mixing  ratio  at  the  sur¬ 
face,  and  V9  is  horizontal  surface  velocity.  Following 
Neelin  et  al  (1987),  if  velocity  is  linearized  about  a 
mean  zonal  flow,  the  evaporation  perturbation  may 
be  written 

£00  =  A?4sgn(us)  (2.3) 

where  A=psCD(q*  —  q)  is  assumed  constant,  u3  is 
the  zonal  velocity  component,  and  the  overbar  and 


Fig.  2.  Evaporation-wind  feedback  growth 
rates  in  the  two- level  model  as  a  func¬ 
tion  of  diffusion  for  q  =  14gkg"*1  (CISK- 
stable  case)  and  A  =  10.0  x  10“6  kgm”2 
s~l  (m/s)*1.  Solid  and  open  circles  denote 
simulations  with  no  mean  wind  and  east¬ 
erly  mean  wind,  respectively.  Open  and 
solid  squares  show  the  e-folding  width  of 
ascending  motion.  Results  are  shown  for 
(iVf,iV)  =  (128,44). 

prime  denote  zonal-mean  and  perturbation,  respec¬ 
tively.  According  to  the  linearized  expression  Eq. 
(2.3),  an  easterly  (westerly)  perturbation  super¬ 
posed  on  an  easterly  (westerly)  mean  flow  acts  to 
enhance  heating  and  vice  versa.  In  the  numerical 
model,  horizontal  surface  velocities  were  approxi¬ 
mated  using  velocity  at  the  lowest  level. 

It  was  shown  by  CD  that  CISK-neutral  instabili¬ 
ties  with  evaporation-wind  feedback  and  conditional 
heating  exhibit  a  weaker  scale-selection  catastrophe 
that  CISK-unstable  modes  (weaker  in  the  sense  that 
although  wet-region  width  contracts  to  zero,  growth 
rate  asymptotes  to  a  constant  as  i/~+0).  Similar 
results  were  obtained  with  the  two-level  beta-plane 
model  as  shown  in  Fig.  2.  For  these  simulations  we 
chose  14  g  kg*1  (below  the  CISK-neutral  point 
without  evaporation-wind  feedback)  and  A  =  10.0x 
10*6kgm*2s*1  (m/s)*1.  This  value  of  A  corre¬ 
sponds  roughly  to  an  evaporation  rate  ~4mm/day 
induced  by  a  surface  perturbation  of  5  ms*1.  Ac¬ 
cording  to  Fig.  2,  growth  rate  increases  slowly  as 
diffusion  decreases,  and  the  width  of  wet  region  de¬ 
creases,  as  in  CD.  (The  e-folding  width  of  vertical 
velocity  was  plotted  in  Figs.  1,  2  instead  of  the  width 
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Fig.  3.  Structure  of  unstable  evaporation-wind  mode  with  conditional  heating  and  easterly  mean  flow  in 
the  twolevel  model  using  q  =  14gkg~\  A  =  10  x  10~6  kgm"3,  and  (Af,  AT)  =  (128,44).  Geopotential 
(thin  solid  and  dashed  lines)  and  horizontal  velocity  vectors  are  shown  at  250  mb  (a)  and  750  mb  (b). 
The  region  of  positive  heating  greater  than  10  %  of  its  maximum  value  lies  inside  the  thick  solid  line. 
The  horizontal  diffusion  in  this  case  was  1  x  106  m2  s~l. 


of  positive  lj  because,  in  the  case  of  evaporation- 
wind  feedback,  there  was  a  broad,  zonally  asym¬ 
metric  'tail’  of  small  lj  extending  far  to  the  west 
of  maximum  convergence.)  We  would  like  to  explore 
the  small  diffusivity  limit  more  rigorously.  However, 
analytic  progress  is  very  difficult  even  in  two  dimen¬ 
sions  (CD)  and  well  behaved  numerical  solutions  re¬ 
quire  some  finite  (and  not  necessarily  small)  amount 
of  diffusivity. 

As  shown  in  Fig.  2,  the  linearized  form  Eq.  (2.3)  is 
unnecessary  for  evaporation-wind  feedback.  Nearly 
identical  eastward-propagating  modes  are  obtained 
using  the  exact  formula  Eq.  (2.2)  with  mean  flow 
neglected.  (See  also  Xie,  et  a/.,  1993).  The  fol¬ 
lowing  discussion  may  help  explain  this  result.  A 
key  feature  of  the  mechanism — using  Eq.  (2.3)  with 
u  <  0 — is  that  warming  occurs  in  the  region  of  per¬ 
turbation  surface  easterlies,  east  of  the  temperature 
maximum.  This  causes  eastward  propagation.  It  is 
irrelevant  that  perturbation  westerlies  cause  cooling 
to  the  west  of  the  convergence  maximum  when  linear 
heating  is  assumed.  If  we  now  return  to  the  exact 
formula  Eq.  (2.2)  and  set  u  -  0,  the  situation  just 
described  would  instead  cause  warming  in  regions 
of  surface  easterlies  and  westerlies.  The  direction 
of  propagation  may  appear  ambiguous,  were  it  not 
for  an  asymmetric  response  on  the  equatorial  beta- 
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plane  (Gill,  1980).  The  response  is  to  heating  that  is 
propagating  and  due  primarily  to  a  superposition  of 
convergences  from  propagating  Kelvin  and  Rossby 
waves.  The  convergence  due  to  the  Kelvin  waves 
is  stronger  and,  consequently,  the  strongest  heating 
moves  eastward.  If  we  now  think  in  terms  of  the  cir¬ 
culation  induced  by  this  eastward  moving  heating, 
the  winds  associated  with  the  Kelvin  response  to  the 
east  are  stronger  than  the  winds  associated  with  the 
Rossby  response  to  the  west.  The  WISHE  mecha¬ 
nism  therefore  favors  an  eastward  moving  unstable 
mode.  There  are  two  crucial  points:  the  relative 
strength  of  the  eastward  moving  heating  and  the 
selective  enhancement  of  such  heating  by  WISHE. 
Thus,  eastward  propagation  is  possible  on  the  beta- 
plane  when  u  =  0.  In  fact,  solutions  are  nearly 
identical  to  those  obtained  with  Eq.  (2.3)  since  the 
Kelvin  response  dominates  that  case  as  well. 

Horizontal  structure,  shown  in  Fig.  3,  is  essen¬ 
tially  an  eastward-propagating  Kelvin  wave  as  in 
unstable  CISK.  A  slight  difference  is  that  heating 
extends  eastward  into  the  region  of  surface  easter¬ 
lies  due  to  the  contribution  from  Au's.  The  solution 
with  no  mean  wind  is  similar,  except  for  a  west¬ 
ward  extension  of  heating  into  the  region  of  sur¬ 
face  westerlies  to  the  west  (Fig.  4).  Using  a  westerly 
mean  flow  with  conditional  heating  produces  a  hy- 
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Fig.  4.  Horizontal  structure  as  in  Fig.  3,  but  with  no  mean  flow. 


Fig.  5.  Horizontal  structure  as  in  Fig.  3,  but  with  westerly  mean  flow. 


brid  Kelvin/Rossby  structure  propagating  westward 
(Fig.  5).  The  hybrid  structure  owes  its  existence  to 
the  fact  that  perturbation  surface  westerlies  enhance 
evaporation  to  the  west  of  the  heating,  dragging  the 


entire  disturbance  westward.  Not  surprisingly,  the 
Rossby  component  of  the  response  is  considerably 
larger  than  in  the  previous  two  examples. 

Solutions  were  also  obtained  with  evaporation- 
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Table  1.  Intrinsic  phase  speeds,  c  (ms'1),  and  e-folding  times,  r  (days),  of  evaporation-wind  feedback 
modes  with  an  easterly  basic  state  and  conditional  heating  as  a  function  of  .4  (kgm'2  s'1  (m/s)'1) 
and  q  (g/kg).  Actual  phase  speeds  relative  to  the  ground  can  be  obtained  by  adding  an  easterly  flow 
to  the  values  in  the  table.  Values  sure  displayed  as  an  ordered  pair  (c,  r).  Results  were  obtained  from 
the  numerical  model  using  (M,N)  =  (128,44)  and  u  =  2.5  x  105  m2  s'1.  The  CISK  neutral  value  in 
the  absence  of  evaporation- wind  feedback  is  q  =  19.1  g/kg.  A  dash  (— )  entry  means  that  no  unstable 
mode  could  be  identified  for  that  set  of  parameters. 

<T=  8  7=  14  q  =  19.1  q  ==  19.5 

A  =  15  x  105  45.3,  2.6  31.8,  2.4  12.8,  1.2  5.8,  0.6 

A  =  10  x  105  45.2,  4.1  31.6,  3.9  12.7,  2.3  4.5,  0.7 

A  =  5  x  105  45.2,  10.7  31.5.  19.2  —  2.6.  0.8 


wind  feedback  in  the  five-level  model,  with  heating 
proportional  to  the  sum  of  low-level  evaporation  and 
moisture  convergence  and  using  a  specified  vertical 
profile.  Results  were  similar  to  those  of  the  two-level 
model  (not  shown),  with  certain  differences  as  eluci¬ 
dated  in  the  next  subsection.  Unlike  unstable  CISK 
modes,  evaporation-wind  feedback  in  a  multi-level 
model  destabilizes  individual  vertical  modes  with¬ 
out  coalescence.  Consequently,  zonal  phase  speeds 
are  similar  to  those  of  free  modes,  and  are  gener¬ 
ally  larger  than  the  propagation  speed  of  CISK  in¬ 
stabilities  (which  are  already  greater  than  observed 
TIO  speeds) .  Examples  of  this  behavior  were  shown 
in  Appendix  B  of  CD.  However,  we  doubt  that 
this  generalization  of  evaporation-wind  feedback  to 
a  multi-level  model — with  specified  vertical  heat¬ 
ing  profile  as  in  wave-CISK — adequately  represents 
the  moist  entropy  budget  as  described  by  Yano  and 
Emanuel  (1991).  In  any  case,  the  phase  speeds  of 
unstable  Kelvin  waves  are  undesirably  large,  so  we 
do  not  discuss  these  modes  further. 

d.  Effective  moist  static  stability 

As  noted  in  the  Introduction,  theories  of  the  TIO 
that  rely  on  CISK  or  evaporation-wind  feedback  fail 
to  reproduce  the  slow  eastward  propagation,  hori¬ 
zontal  structure  (including  subtropical  gyres),  and 
hierarchy  of  convection  within  intraseasonal  oscil¬ 
lations.  The  last  two  observations  suggest  the  im¬ 
portance  of  latitudinal  shear  (Dunkerton,  1990)  and 
group  velocity  of  synoptic-scale  waves  (Liebmann 
and  Hendon,  1990),  and  show  that  a  purely  Kelvin- 
wave  interpretation  of  the  TIO  is  inadequate. 

Concerning  the  first  point,  a  general  shortcom¬ 
ing  of  CISK-unstable  Kelvin  waves  in  a  multi-level 
model  with  specified  vertical  heating  profile  is  that 
moist  stable  modes  (contiguous  to  dry  free  modes) 
coalesce  to  form  an  instability  at  unrealistically  high 
phase  speeds  similar  to  those  of  dry  modes.  When 
evaporation- wind  feedback  is  added,  CISK-stable  in¬ 
stabilities  remain  fast  (CD). 

An  exception  to  this  statement  occurs  in  the  two- 
level  model  destabilized  by  evaporation-wind  feed¬ 
back.  In  this  case  phase  speeds  can  be  made  ar¬ 


bitrarily  small,  according  to  linear  theory,  by  ap¬ 
proaching  the  CISK-neutral  point  from  below  and 
simultaneously  reducing  A.  Reduction  of  phase 
speed  (and  of  latitudinal  scale)  is  due  to  an  ef¬ 
fective  ‘moist  static  stability'  since  heating  in  this 
model  is  proportional  to  vertical  velocity.  At  the 
CISK-neutral  point,  latent  heat  release  is  exactly 
balanced  by  adiabatic  cooling.  This  state  of  af¬ 
fairs  is  sometimes  referred  to  as  4  free- ride’  balance 
(Praedrich  and  McBride,  1989)  and  represents  ap¬ 
proximate  convective  neutrality  of  the  tropical  at¬ 
mosphere.  It  is  not  limited  to  a  two-level  model,  as 
latent  heat  release  can  be  equal  to,  or  some  fixed 
fraction  of,  adiabatic  cooling  at  several  vertical  lev¬ 
els. 

We  now  inquire  whether  a  similar  reduction  of 
phase  speed  occurs  in  a  beta-plane  model  with  con¬ 
ditional  heating.  The  question  is  important  because, 
in  the  real  atmosphere,  effective  ‘moist  stability5  is 
a  function  of  the  prevailing  vertical  motion  and  la¬ 
tent  heat  release.  We  could  either  assume  some  fixed 
zonal  variation  of  stability  representing  the  Walker 
circulation  and  underlying  SST  gradient,  or  make 
the  stability  solution-dependent,  z.e.,  with  condi¬ 
tional  heating  in  a  uniform  basic  state.  The  sec¬ 
ond  approach  was  adopted  here  although  the  first  is 
equally  relevant  to  the  atmosphere. 

For  simplicity,  solutions  were  sought  with  lin¬ 
earized  evaporation-wind  feedback  and  easterly 
mean  flow  in  a  two-layer  model.  Table  1  displays 
the  growth  rate  and  intrinsic  phase  speed  of  unstable 
Kelvin  waves  near  the  CISK-neutral  point  obtained 
from  numerical  integrations  with  (M,  N)  =  (128, 44) 
and  v  =  2.5  •  10s  m2  s~x.  It  is  apparent  that  phase 
speeds  decrease  quite  dramatically  approaching  the 
neutral  point,  and  continue  to  decrease  entering  the 
CISK-unstable  region.  (The  advective  influence  of 
low-level  u  was  included  in  these  simulations,  but 
not  in  Table  1.  For  an  easterly  mean  flow,  ground- 
based  phase  speeds  c  =  c  4-  u  are  less  than  the  tab¬ 
ulated  values  of  c.)  Similar  results  were  recently 
obtained  by  Kirtman  and  Vemekar  (1993)  using  a 
Kelvin  wave  model  with  unconditional  heating,  and 
by  Xie  et  al.  (1993).  Horizontal  structure  in  most 
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cases  was  that  of  a  Kelvin  wave,  but  with  nonzero 
meridional  velocity  as  c  — *  0.  Tabie  1  suggests  that 
the  phase  speed  of  the  whole  disturbance,  includ¬ 
ing  dry  region,  is  controlled  by  an  effective  ‘moist 
stability'  of  the  wet  region.  Phase  speeds  were  also 
observed  to  decrease  slowly  with  decreasing  .4  as  ex¬ 
pected  from  linear  theory.  The  CISK-stable  region 
near  .4  =  0  was  problematic  since  nonzero  dilfusion 
was  always  required  for  computational  purposes  and 
growth  rates  are  expected  to  be  very  weak.  We  were 
therefore  unable  to  verify  whether,  as  expected,  the 
Kelvin  modes  described  in  Table  1  are  contiguous  to 
the  westward-propagating  hybrid  structure  obtained 
for  u  >  0. 

3.  Conclusions 

The  effect  of  conditional  heating  on  unstable 
Kelvin  waves  generated  by  CISK  and  evaporation- 
wind  feedback  was  investigated  with  multi-level 
models  on  the  equatorial  beta-plane.  In  agreement 
with  previous  analytic  and  numerical  results  (Crum 
and  Dunkerton,  1992),  the  scale-selection  catastro¬ 
phe  of  linear  heating  is  modified  but  not  eliminated 
by  conditional  heating.  There  is  only  one  region 
of  upward  motion,  but  it  contracts  to  the  smallest 
zonal  scale  allowed  by  diffusion  or  model  resolution. 
The  pathology  is  most  severe  in  the  CISK-unstable 
case.  With  evaporation-wind  feedback,  Kelvin-wave 
instability  is  possible  apart  from  CISK;  growth  rate 
in  this  case  increases  slowly  with  decreasing  diffu¬ 
sion,  but  contraction  is  still  observed.  We  note  in¬ 
cidentally  that  eastward  moving  evaporation-wind 
modes  can  be  generated  on  the  beta-plane  with  no 
mean  wind,  contrary  to  previous  assumptions  re¬ 
quiring  an  easterly  background  flow.  In  a  westerly 
mean  flow,  a  westward-propagating  hybrid  structure 
is  obtained. 

While  evaporation-wind  feedback  is  an  improve¬ 
ment  over  CISK,  as  far  as  scale  selection  is  con¬ 
cerned,  the  suggestion  that  unstable  Kelvin  waves 
explain  the  tropical  intraseasonal  oscillation  is  in¬ 
adequate.  The  first  and  most  obvious  problem  is 
that  zonal  phase  speeds  are  too  fast  compared  to 
the  observed  propagation  of  intraseasonal  oscilla¬ 
tions  and  superclusters.  In  multi-level  models  with 
specified  vertical  heating  profile  this  is  due  to  the 
coalescence  or  destabilization  of  moist  stable  modes 
having  phase  speeds  similar  to  those  of  dry  Kelvin 
modes.  In  a  two-level  model,  on  the  other  hand, 
destabilization  appears  possible  at  an  arbitrarily 
small  phase  speed  when  the  effective  moist  stabil¬ 
ity  is  reduced  to  zero  and  the  parameter  A  con¬ 
trolling  the  strength  of  evaporation-wind  feedback  is 
small.  This  state  of  affairs  represents  approximate 
convective  neutrality  of  the  tropical  atmosphere.  In 
the  limit  of  neutral  stability,  unstable  waves  on  the 
beta-plane  depart  from  a  pure  Kelvin  structure  as 
required  by  linear  theory.  While  slow  eastward  prop¬ 
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agation  appears  possible  in  this  case,  we  are  both¬ 
ered  by  the  sensitivity  to  parameters  (q.  A)  and  won¬ 
der  how  atmospheric  conditions  could  ever  conspire 
to  produce  a  robust  TIO  with  well-defined  spectral 
peak. 

Two  additional  problems  are  the  failure  to  gener¬ 
ate  any  significant  subtropical  structure  or  hierarchy 
of  convection  within  active  phases  of  the  TIO.  These 
observations  suggest,  respectively,  the  importance  of 
latitudinal  shear  (Dunkerton,  1990)  and  group  ve¬ 
locity  of  synoptic-scale  waves  (Liebmann  and  Hen¬ 
don,  1990).  The  fact  that  moisture  convergence  and 
evaporation-wind  feedback  actually  take  place  in  a 
hierarchy  of  squall  lines,  mesoscale  convective  com¬ 
plexes,  and  synoptic  waves  also  reveals  the  inade¬ 
quacy  of  a  cumulus  parameterization  that  depends 
only  on  large-scale  variables. 

Acknow  led  gement  s 

This  research  was  supported  by  the  National 
Science  Foundation.  Grants  ATM-8819582,  ATM- 
9123797,  ATM-9213437,  and  by  the  Air  Force  Of¬ 
fice  of  Scientific  Research,  Contracts  F49620-89-C- 
0051  and  F49620-92-C-0033.  Numerical  computa¬ 
tions  were  performed  at  the  National  Center  for  At¬ 
mospheric  Research,  which  is  sponsored  by  the  Na¬ 
tional  Science  Foundation. 

References 

Chang,  C.-P.  and  H.  Lim.  1988:  Kelvin  wave-CISK:  a 
possible  mechanism  for  the  30-60  day  oscillations. 
J.  Atmos .  Set,  45,  1709-1720. 

Charney,  J.G.  and  A.  Eliassen,  1964:  On  the  growth  of 
the  hurricane  depression.  J .  Atmos.  Set,  21,  68-75. 
Crum,  F.X.  and  D.E.  Stevens,  1983:  A  comparison  of 
two  cumulus  parameterization  schemes  in  a  linear 
model  of  wave-CISK.  J.  Atmos.  Set,  40,  2671-2688. 
Crum,  F.X.  and  T.J.  Dunkerton,  1992:  Analytic  and  nu¬ 
merical  models  of  wave-CISK  with  conditional  heat¬ 
ing.  J.  Atmos.  Set,  49,  1693-1708. 

Dunkerton,  T.J.,  1990:  Eigenvalues  and  horizontal 
structure  of  divergent  barotropic  instability  originat¬ 
ing  in  tropical  latitudes.  J.  Atmos.  Sci .,  47,  1288- 
1301. 

Emanuel,  K.A.,  1987:  An  air-sea  interaction  model  of 
intraseasonal  oscillations  in  the  tropics.  J.  Atmos. 
Set ,  44,  2324-2340. 

Fraedrich,  K.  and  J.L.  McBride,  1989:  The  physical 
mechanism  of  CISK  and  the  free-ride  balance.  J. 
Atmos.  Set ,  46.  2642-264 8. 

Gill,  A.E.,  1980:  Some  simple  solutions  for  heat-induced 
tropical  circulations.  Q.  J .  Roy.  Meteor.  Soc 106, 
447-462. 

Hayashi,  Y.-Y.  and  A.  Sumi.  1986:  The  30-40  day  oscil¬ 
lations  simulated  in  an  ‘aquaplanet’  model.  J.  Me¬ 
teor.  Soc.  Japan,  64,  451-467. 

Houghton.  J.T.,  1986:  The  Physics  of  Atmospheres.  Sec¬ 
ond  edition.  Cambridge  University  Press.  271  pp. 


18 


Journal  of  the  Meteorological  Society  of  Japan 


Vol.  72.  No.  1 


Kirtman.  B.  and  A.  Vernekar,  1993:  On  wave-CISK 
and  the  evaporation-wind  feedback  for  the  Madden- 
Julian  oscillation.  J.  Atmos.  Sci.,  50,  2811-2814. 

Knutson.  T.R.  and  K.M.  Weickmann,  1987:  30-60  day 
atmospheric  oscillation:  composite  life  cycles  of  con¬ 
vection  and  circulation  anomalies.  Mon.  Wea.  Rev., 
115.  1407-1436. 

Lau,  K.-M.  and  L.  Peng,  1987:  Origin  of  low  fre¬ 
quency  (intraseasonai)  oscillations  in  the  tropical  at¬ 
mosphere.  Part  I:  The  basic  theory.  J.  Atmos.  Sci., 
44,  950-972. 

Lau,  N.-C.,  I.M.  Held  and  J.D.  Neelin,  1988:  The 
Madden- Julian  oscillation  in  an  idealized  general  cir¬ 
culation  model.  J.  Atmos.  Sci.,  45.  3810-3832. 

Liebmann,  B.  and  H.H.  Hendon,  1990:  Synoptic-scale 
disturbances  near  the  equator.  J.  Atmos.  Sci..  47. 
1463-1479. 

Madden.  R.A.  and  P.R.  Julian,  1971:  Detection  of  a  40— 
50  day  oscillation  in  the  zonal  wind  in  the  tropical 
Pacific.  J.  Atmos.  Sci.,  28,  702-708. 

Madden,  R.A.  and  P.R.  Julian,  1972:  Description  of 
global  scale  circulation  ceils  in  the  tropics  with  a  40- 


50  day  period.  J.  Atmos.  Set 29.  1109-1123. 

Neelin,  J.D.,  I.M.  Held  and  K.H.  Cook,  1987: 
Evaporation-wind  feedback  and  low  frequency  vari- 
abilty  in  the  tropical  atmosphere.  J.  Atmos.  Set .. 
44,  2341-2348. 

Rui,  H.  and  B.  Wang,  1990:  Development  characteris¬ 
tics  and  dynamic  structure  of  tropical  intraseasonai 
convection  anomalies.  J.  Atmos.  Set.,  47,  357-379. 

Swinbank,  R.,  T.N.  Palmer  and  M.K.  Davey,  1987:  Nu¬ 
merical  simulations  of  the  Madden  and  Julian  oscil¬ 
lation.  J.  Atmos .  Sci.,  45,  774-788. 

Xie,  S.-P.,  A.  Kubokawa  and  K.  Hanawa,  1993:  Nonlin¬ 
ear  evaporation-wind  feedback  mode.  Preprint  Vol¬ 
ume,  Ninth  Conference  on  Atmospheric  and  Oceanic 
Waves  and  Stability,  San  Antonio.  Texas.  385-386. 

Yano,  J.-I.  and  K.  Emanuel,  1991:  An  improved  model 
of  the  equatorial  troposphere  and  its  coupling  with 
the  stratosphere.  J.  Atmos.  Sci.,  48.  377-389. 

Yoshizaki,  M.,  1991:  Selective  amplification  of  the 
eastward-propagating  mode  in  a  positive-only  wave- 
CISK  model  on  an  equatorial  beta  plane.  J.  Meteor. 
Soc.  Japan,  69,  353-373. 


Francis  X.  Crum  and  Timothy  J.  Dunkerton 

(  Northwest  Research  Associates  ) 


CISK  >  M^-lRy  <  -  ?  C £ 

ii£  b<x><r> — ?xh  9.  %<T)7k¥-X'r  — 

J: -oXfct  CISK  iz X  S^&Z%®&tz  i>  >  b^MX'b 

So  ZVt§&,  tfflg mu 

<  %6lZ-otliQ>6J<?frlz3:§.  <  & 9 /Js?  <  a.±x', 

% t  StiX '.'fzMMt -  Y't-vtT'&Z*:-  Ytf,  i  b  hbtiZX) 

ffynStitZo 

z  * -*&1RIZM  LX CISK  X  9  Uk  Lfr  L>  (1)  7m^My^- 

q  CISK  (2)  M&-®.y  <  -  K'<y 


44 


Horizontal  Buoyancy  Flux  of  Internal  Gravity  Waves  in  Vertical  Shear 


Timothy  J.  Dunkerton 


45 


46 


Journal  of  the  Mciroroiofftcal  Society  of  Japan .  Vol.  73,  Xo.  3.  pp.  747-755.  1995 


Horizontal  Buoyancy  Flux  of  Internal  Gravity  Waves  in  Vertical  Shear 


By  Timothy  J.  Dunkerton 

Northwest  Research  Associates.  P.0.  Box  3021,  Bellevue ,  WA  98009,  (J.S.A. 

(Manuscript  received  31  October  1994 >  revised  form  18  April  1995) 

Abstract 

The  equations  of  motion  describing  two-dimensional  internal  gravity  waves  were  analyzed  to  derive  an 
expression  for  the  horizontal  buoyancy  flux — i.e.,  the  correlation  between  buoyancy  (or  temperature)  and 
horizontal  velocity— in  various  cases  involving  vertical  shear.  It  is  shown  that,  although  this  correlation 
vanishes  at  zeroth  order  for  a  single  wave,  its  first  order  contribution  is  nonzero  due  to  shear,  even  for  steady, 
conservative,  incompressible  waves.  Departures  from  steady,  conservative  motion  or  quasi-compressibility 
also  cause  a  nonzero  correlation.  Several  cases  were  analyzed  and  some  numerical  results  obtained  for 
waves  approaching  a  critical  layer  of  reduced  intrinsic  phase  speed.  With  weak  shear,  the  buoyancy  flux 
is  small  relative  to  vertical  momentum  flux,  as  expected  from  the  perturbation  theory.  Strong  vertical 
shear  enhances  the  buoyancy  flux  within  the  shear  zone  and  causes  partial  reflection  beneath,  producing  a 
nonzero  correlation  (at  zeroth  order)  in  this  region.  These  effects  may  explain  recent  observations  of  zonal 
wind  and  temperature  cospectra  in  the  equatorial  lower  stratosphere. 


1.  Introduction 

Gravity  waves  are  important  for  the  transport  of 
momentum,  heat  and  constituents  in  the  middle  at¬ 
mosphere.  Their  vertical,  rather  than  horizontal, 
transport  is  emphasized  in  the  literature  since  the 
waves  grow  with  height  and  are  able,  for  example, 
to  transport  a  relatively  small  amount  of  momen¬ 
tum  from  the  lower  to  upper  atmosphere,  causing 
significant  mean  flow  acceleration  in  the  mesosphere 
(Lindzen,  1981;  Fritts,  1984).  Little  attention  has 
been  paid  to  the  horizontal  transport  of  momentum 
by  gravity  waves,  and  even  less  to  the  horizontal 
transport  of  heat  which  is  formally  smaller  (by  an 
order  of  magnitude)  in  WKB  theory. 

Whether  observations  assign  a  significant  role  to 
the  horizontal  transport  is  uncertain,  but  the  in¬ 
terpretation  of  observations  may  require  an  under¬ 
standing  of  horizontal  fluxes.  Recently  Sato  et  ai 
(1994)  discovered  that  zonal  wind  (uf)  and  tempera¬ 
ture  ( T" )  fluctuations  in  the  equatorial  lower  strato¬ 
sphere,  with  period  ~l-3  days,  display  a  remark¬ 
able  correlation  that  changes  sign  according  to  the 
phase  of  the  quasi-biennial  oscillation  (QBO).  Equa¬ 
torial  Kelvin  waves  were  previously  identified  from 
the  quadrature  spectrum  Q{uf which  has  a 
significant  peak  at  ~5-20  days  in  the  equatorial 
lower  stratosphere  (Wallace  and  Kousky,  1968:  Hi- 
rota,  1978;  Salbv  et  al .  1984;  Maruyama,  1991).  Re- 
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cent  observations  suggest  an  important  contribution 
to  the  momentum  flux  from  disturbances  with  pe¬ 
riods  below  5  days  (Pfister  et  ai ,  1993;  Maruyama, 
1994;  Bergman  and  Salbv,  1994).  Noting  the  pecu¬ 
liar  behavior  of  the  cospectrum  in  their  data.  Sato 
et  al.  (1994)  question  whether  the  observed  phase 
relation  between  u!  and  V  is  consistent  with  Kelvin 
or  internal  gravity  waves. 

Sato  et  al.  were  mainly  concerned  with  the  inter¬ 
pretation  of  short-period  disturbances,  but  it  is  obvi¬ 
ous  from  their  analysis  that  a  similar  problem  arises 
with  respect  to  ~5~20  day  disturbances  traditionally 
identified  as  Kelvin  waves.  At  Singapore,  the  inte¬ 
grated  ~o-20  day  quadrature  Q(u',Tf)  and  cospec¬ 
trum  C(u'.T')  power  are  greatly  enhanced  during 
the  onset  of  the  QBO  westerly  phase  at  each  level.1 
Their  ratio  is  typically  1.5:1.  The  spikes  of  ^5-20 
day  quadrature  power  are  consistent  with  those  of 
eastward-propagating  waves  encountering  a  critical 
layer  in  descending  QBO  westerlies.  Meridional  ve¬ 
locity  fluctuations  (z/)  are  uncorrelated  (or  negligi¬ 
ble),  supporting  a  Kelvin- wave  interpretation.  But 
the  cospectrum  of  these  waves,  and  that  of  short- 
period  disturbances,  requires  explanation. 

The  purpose  of  this  paper  is  to  discuss  the  rela¬ 
tion  between  horizontal  wind  and  temperature  (or. 
equivalently,  buoyancy)  for  two-dimensional  internal 

1  As  shown  in  Fig.  10  of  Sato  et  al.  (1994).  Spectral  analy¬ 
sis  of  Singapore  rawinsonde  data  by  the  author  supports 
their  conclusions. 
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gravity  waves  in  vertical  shear.  It  is  shown  theo¬ 
retically  and  numerically  that  the  correlation  uT 
is  nonzero  and  proportional  to  the  shear— even  for 
steady,  conservative,  incompressible  waves.  For  ver¬ 
tical  shears  typical  of  the  QBO,  the  effect  of  com¬ 
pressibility  is  small,  other  than  causing  the  waves  to 
grow  with  height.  Realistic  damping  likewise  does 
not  significantly  alter  the  effect  of  shear.  Numerical 
solutions  demonstrate  that  superposition  of  upward 
and  downward  propagating  waves  creates  a  signif¬ 
icant  horizontal  buoyancy  flux  beneath  the  shear 
zone,  consistent  with  zeroth  order  wave  solutions. 
Finally,  we  note  that  although  the  horizontal  flux 
is  formally  nondivergent  in  the  situation  examined 
here,  real  wavepackets  are  of  finite  extent  such  that 
the  horizontal  transport  of  buoyancy  is  nonzero. 

2.  Theoretical  results 

The  equations  of  motion  for  two-dimensional, 
nonrotating,  hydrostatic,  quasi-compressible,  linear 
perturbations  on  a  zonal-mean  flow  are 


Dtu'  4-  w'uz  +  p'x  =  X'  (2.1a) 

Dtd’.  +  w'N2  =  Q'  (2.1b) 

ux  d-  -“( Pow')z  =  0  (2.1c) 


in  the  notation  of  Holton  (1975),  where  u,  w  are 
zonal  and  vertical  velocity,  o  is  geopotential,  Oq  — 
psexp(—z/H)  is  basic  state  density,  and  X.  Q  rep¬ 
resent  sources  or  sinks  of  momentum  and  heat. 
Primed  quantities  are  deviations  from  the  zonal 
mean,  denoted  by  an  overbar;  subscripts  x,  z  in¬ 
dicate  differentiation,  and 


Dt  = 


dt  '  ^  dx 


(2.2) 


a.  Wave  action  law 

Following  Andrews  and  McIntyre  (1976),  parcel 
displacements  and  auxiliary  variables  may  be  de¬ 
fined  as 


D£  =  ue 

(2.3a) 

AC'  =  w' 

(2.3b) 

to 

II 

<o 

(2.3c) 

where 

ue  =s  v!  -r  C'uz- 

(2.4) 

Assuming  that  zonally  averaged  quantities  are  tem¬ 
porally  constant  to  leading  order  in  wave  amplitude, 

(2.1a-c)  may  be  written 

Dtue  +  <j>'x  =  X' 

(2.5a) 

4>'z  +  C'N2  =  q' 

Ci  +  —{PoC,')z  =  o. 

Po 

(2.5b) 

(2.5c) 

(£V)t  +  (&o'x  +  Cxo'z)  =  £i-Y'  +  c;  q‘  (2.6) 

noting  that 

)t  - 

=  -  a«  =  (&u*)t 

and  CiC  =  0.  Substitution  of  (2.5c)  in  (2.6)  implies 
that 

+  =  ^  +  (2.7) 

This  equation  describes  conservation  of  wave  action, 
a  quadratic  (second-order)  quantity  in  wave  am¬ 
plitude.  A  three-dimensional  version  of  this  equa¬ 
tion  (in  y.  z,  t)  was  derived  by  Andrews  and  Mcln- 
tyre(1976): 

dA  1 

»+^VB  =  ^  (2-8) 

where  A  is  wave  action  density,  B  is  wave  action 
flux,  and  T  represents  the  departure  from  conser¬ 
vative  motion.  For  steady,  conservative  waves,  B  is 
nondivergent.  The  famous  theorem  of  Eliassen  and 
Palm  (1961),  extended  by  Andrews  and  McIntyre, 
relates  the  mean  flow  tendency  ut  to  the  divergence 
of  B  and  other  departures  from  steady,  conservative 
motion. 

The  vertical  flux  in  (2.7)  is  akin  to  a  ‘radiation 
stress’  describing  the  impulse  to  an  overlying  layer 
of  fluid  by  pressure  forces  o'  averaged  over  a  wave- 
distorted  fluid  interface  z  +  ^{x,z,  t).  (Refer  to  An¬ 
drews  et  al.,  1987,  p.  325  for  a  diagram.)  This  stress 
is  nonzero,  since  the  quantities  Q'x,  o'  are  in  phase 
for  internal  gravity  waves,  but  it  does  not  vary  in 
height  if  the  waves  are  steady  and  conservative. 

The  radiation  stress  may  be  related  to  the  vertical 
flux  of  horizontal  momentum  as 

=  C{Dtue  -  X') 

=  (CV)t-^+i(C;2)tuI-C7F  (2.9) 

from  which  it  is  apparent  that  the  momentum 
flux  differs  from  the  radiation  stress  only  on  ac¬ 
count  of  departures  from  steady,  conservative  mo¬ 
tion.  In  general,  the  components  of  Eliassen-Palm 
flux  satisfy  this  type  of  relationship  which — together 
with  the  wave  action  law  (2.7)  and  mean  ten¬ 
dency  equations — leads  to  a  ‘nonacceleration’  the¬ 
orem  (Andrews  and  McIntyre,  1976). 

b.  Correlation  of  zonal  wind,  and  temperature 
For  hydrostatic  waves  the  temperature  perturba¬ 
tion  is  proportional  to  o',  (Holton.  1975).  The  cor¬ 
relation  of  zonal  wind  and  temperature  is  therefore 
proportional  to 


Multiplication  of  (2.5a, b)  by  £.,  Q'x  gives 


pou'O'.  =  p0u'{q'  -  C'-V2) 
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=  pou'q'  - 
=  pou/c?  - 


(2.10) 


for  each  harmonic  exp  i(kx  -  ut),  where  k  is  zonal 
wavenumber,  ^  is  frequency,  c  =  uj/k  is  phase  speed, 
and  c  =s  c  -  u  is  intrinsic  phase  speed.  (Use  of  the 
harmonic  form  is  convenient,  but  unnecessary  to  the 
argument  as  shown  at  the  end  of  this  subsection.) 
The  continuity  equation  (2.1c)  was  inserted  on  the 
last  line  of  (2.10),  together  with  the  vertical  dis¬ 
placement  (2.3b),  z.e., 


(2.11) 


for  steady  waves. 

A  similar  expression  for  p§u‘ o' z  may  be  derived 
using  the  streamfunction 

-pou'  =  (poib')z  (2.12a) 

p0w'  =  (po^)x  (2.12b) 

From  (2.11),  -cQf  =  so  that 

Pou'ftz  =  Pou'itf  -  C'-V2) 

_  iV2 _ 

=  Pou'q' - -(poV')zPo'tp'  (2.13) 

Poc 

which  agrees  with  (2.10)  on  account  of  the  harmonic 
form  assumed  in  the  derivation. 

Apart  from  a  thermal  damping  term,  the  corre¬ 
lation  of  zonal  wind  and  temperature  evidently  de¬ 
pends  on  the  vertical  variation  of  wave  amplitude 
squared,  as  determined  by  p oifr'  or  pow': 

_  /V2  _ 

pou'<p’z  =  +  O(q')  (2.14) 

Consider  now  the  special  case  of  a  steady,  conser- 
vative  upward- propagating  wave  in  a  basic  state 
(po,  u,N)  independent  of  height,  for  which 


w  = 


(2.15) 


where  m  is  vertical  wavenumber.  Substituting  for  a 
factor  w'  in  (2.14), 


Pou'<j>'z 


N2 

2pok2c 


(Pou'tu'&c/iV)* 


(2.16) 


The  quantity  pou'w'  is  independent  of  height  by  the 
wave  action  law  (2.7)  and  (2.9),  as  are  po,k,c,  and 
N,  so  the  correlation  pQu'(p'z  vanishes  identically. 
This  is  obvious  from  the  perfect  correlation  (or  anti¬ 
correlation)  of  u*  and  0'  to  begin  with,  and  the  fact 
that  <pz  is  exactly  in  quadrature  with  0'. 

Another  special  case  involves  the  superposition  of 
identical  upward  and  downward  waves  with  unequal 


amplitude: 

ip'  =  Re[wr  exp  i(kx  4-  mz  -  u ;t) 

4-iiq  exp  i(kx  —  mz  —  ut)] 

=  Rew^  [exp  i(kx  4-  mz  —  ut) 

+Rexpi(kx  -  mz  —  ut  +  <5)]  (2.17a) 

ip'2  =  ^I^tI2  [l  +  R2  +  2Rcos(2mz  -  5)]  (2.17b) 

The  correlation  pou'oi  oscillates  in  the  vertical 
about  zero,  and  if  R  =  1,  the  momentum  flux  u'w' 
vanishes.  Perfect  reflection  is  analogous  to  a  zon- 
ally  propagating  circulation  cell  in  the  x  —  z  plane 
confined  between  two  rigid  plates,  such  that  verti¬ 
cal  velocity  vanishes  on  the  boundaries.  Suppose 
the  propagation  is  eastward,  and  the  motion  adia¬ 
batic.  Behind  the  upwelling  there  is  cold  air,  due 
to  adiabatic  cooling,  which  partially  overlaps  the 
eastward  (westward)  motion  at  lower  (upper)  lev¬ 
els  required  by  continuity.  Warm  air  is  ahead  of 
the  upwelling,  overlapping  the  westward  (eastward) 
motion  at  lower  (upper)  levels.  The  correlation  v/T' 
is  therefore  negative  (positive)  in  the  lower  (upper) 
layer,  and  zero  at  the  midpoint  and  boundaries.  By 
contrast,  equal  amounts  of  positive  and  negative  mo¬ 
mentum  flux  cancel  identically  in  each  layer  when 
the  zonal  average  is  taken. 

The  WKB  approximation  provides  a  third  exam¬ 
ple: 


ip'(x,z,t)  =s  ReA(Z)  expi$(x,z,  t)  (2.18) 


where  k  =  d$/dx ,  m  =  d$/dz,  u ;  =  - d$/dt ,  and 
Z  =  pz  is  a  'slow5  vertical  coordinate  such  that  A,  m 
vary  slowly  in  height,  i.e.,  as  functions  of  Z.  Then 


ip'xW 


/ 

X 


\pk  Im(A*A') 


(2.19) 


This  is  analogous  to  the  first  case,  but  using  a  local 
value  of  m,  together  with  O(p)  correction  that  varies 
slowly  in  height.  If  departures  from  steady,  conser¬ 
vative  motion  are  0(p?)  or  less,  the  momentum  flux 
is  invariant  to  this  order,  and  a  useful  result  is 


_  JV2  _ 

Pou'Qi  =  -—-(pou'w*)  *  (pokc/N)z  4-  O(pr) 

(2.20) 


In  (2.18-20)  we  have  also  assumed  that  the  basic 
state  density  variation  is  0(p2),  z.e.,  that  the  mean 
flow  varies  rapidly  over  a  density  scale  height,  as  ob¬ 
served  in  the  QBO.  Supposing  for  consistency  that 
variations  of  N  are  likewise  small,  (2.20)  reduces  to 


_  ,y  _ . 

PqU'0'z  as  —tPou'w'Uz  4-  0(p2).  (2.21) 

The  correlation  between  zonal  wind  and  tempera¬ 
ture  is  proportional  to  vertical  shear,  an  order  of 
magnitude  smaller  (in  p)  than  the  momentum  flux. 


49 


750 


Journal  of  the  Meteorological  Society  of  Japan 


Vol.  73.  No.  3 


For  upward  waves,  the  ratio  of  momentum  flux  to  in¬ 
trinsic  phase  speed  is  positive,  so  the  sign  of  pouVf 
in  this  case  is  determined  by  the  vertical  shear.  Up¬ 
ward  propagating  waves  with  phase  speeds  of  oppo¬ 
site  sign  will  contribute  in  the  same  way  to  the  cor¬ 
relation  of  zonal  wind  and  temperature,  but  their 
momentum  fluxes  will  tend  to  cancel.  For  down¬ 
ward  waves,  the  sign  of  the  correlation  is  unchanged 
since  the  sign  of  m  in  (2.15),  and  of  u'w'/c ,  are  both 
reversed. 

Under  the  assumptions  leading  to  (2.21)  the  same 
result  may  be  obtained  from  the  WKB  approxima¬ 
tion  noting  that  |^'|  cc  |m|“ 1/2,  substituting  the 
value  of  m  from  (2.15)  together  with  this  propor¬ 
tionality  into  (2.14).  If  the  shear  is  positive,  ip* 
decreases  (increases)  with  height  for  waves  with  pos- 
itive  (negative)  intrinsic  phase  speed;  the  correlation 
Pqu'T'  is  therefore  positive.  If  the  shear  is  negative, 
ip'2  increases  (decreases)  with  height  for  waves  with 
positive  (negative)  intrinsic  phase  speed;  the  corre¬ 
lation  is  then  negative.  It  is  assumed  that  the  sign 
of  c  is  invariant,  i.e that  waves  cannot  propagate 
through  critical  levels. 

If  the  unsteadiness,  damping  and  density  varia¬ 
tion  are  0{p),  a  more  complicated  expression  for  the 
single  upward  wave  may  be  obtained  from  (2.14)  by 
writing 


w'2  —  u'w'kc/N  -f-  0(/i) 


(2.22) 


so  that 

_  jy2  _ 

P^'Vz  =  (Pou'w'kc/N); 

+pou'q'  +  0{p?) 

N2 _ 

=  -2fa.u,w'(p0c/N)z 

N  _ 

-2^(Pou'w')z 

+P0U3?  +  0(n2)  (2.23) 

The  first  term  represents  variations  of  the  basic  state 
at  O(p);  the  second  term  contains  all  the  departures 
from  steady,  conservative  motion  at  O(p)  implied  by 
(2.7)  and  (2.9). 

Three  simple  examples  demonstrate  the  effects  of 
density  variation,  damping,  and  saturation  at  O(p): 
1)  If  the  waves  are  steady  and  conservative,  and  u, 
N  are  constant  in  height  (but  not  p0),  then  from 
(2.14,15) 


Pou'0'z  = 


_ fa 

N  _ 

— po-u'w'  +0(pr) 


The  incompressible  value  of  m  was  substituted  in 
this  equation.  2)  If  the  waves  are  steady  but  damped 
by  mechanical  and  thermal  relaxation,  with  equal 
coefficients  au  =  ctr  =  and  the  basic  state  is 


again  constant  in  height  (including  p0).  then  it  may 
shown  from  the  wave  action  law  (2.7)  and  zeroth 
order  wave  solutions  that 


B{z)  =  B0exp-P(z) 


P(z)  = 


D(z')dzf 


(2.25a) 

(2.25b) 


D(z)  =  aN/kc2 


where  B  =  pou'w',  In  this  case  B  decays  with 
height,  such  that  Bz  =  - BPZ  =  —  BD,  therefore 


Pou'd)'z  =  ~2 p^£PoBzlv  +  poU'q/  + 


2  kHD 


Pou'w'  +  Pou'q'  +  0(p2)  (2.26) 


where  Hq  =  kc2/Na  (Dunkerton,  1981).  Note  the 
similarity  to  the  previous  case,  with  H  replaced 
by  Hq.  This  formula  could  easily  be  extended  to 
include  0{p)  basic  state  variations.  3)  Following 
Lindzen  (1981),  suppose  that  the  waves  are  steady 
and  conservative,  growing  with  height  until  'saturat¬ 
ing5  at  some  level,  the  criterion  for  saturation  being 
that  the  waves  are  marginally  unstable  to  convec¬ 
tion.  Above  this  level, 

B  =  Bs  =  po  (2.27) 

so  that 

_  JV  2 

Pou'<p':  =  --—^(plk2c4/N2)z  +  pou'q' 


2pouzc2  -h  pou'q'. 


(2.28) 


On  the  last  line  it  was  assumed  for  illustration  that 
the  mean  flow  varies  rapidly  over  a  density  scale 
height.  Note  that  the  buoyancy  flux  at  the  breaking 
level,  apart  from  the  thermal  damping  term,  is  four 
times  that  implied  by  (2.21),  but  decreases  rapidly 
above  this  level.  More  complete  analysis  of  satura¬ 
tion  may  require  consideration  of  the  explicit  forms 
of  X ',  Q'  that  cause  the  waves  to  saturate  in  the 
first  place. 

The  contribution  from  pou'q'  may  be  evaluated 
for  steady  waves  noting  that  - cq'x  =  Q',  or 

k2cq'  =  —olt<Pzx  =  arN2^x  =  —arN2wf/c 


so  that 


Pou'q'  = 


oltN2 


k 2£2  POU'W'. 


The  sign  is  opposite  to  the  momentum  flux.  The 
effect  of  thermal  damping  is  to  shift  tempera¬ 
ture  anomalies  forward  in  the  direction  of  intrinsic 
phase  propagation,  producing  a  negative  correlation 
u'T'  <  0  if  u'w'  >  0. 
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Finally,  a  nondimensional  measure  of  the  horizon¬ 
tal  buoyancy  flux  may  be  obtained  from  the  ratio 


ku'oL 

u'o'zx 


kc  u'oL 
N2  u! w 1 


0{p?) 


Re(0V)  olt  2^ 
T'~  7  — :  — rr  -r  0(m  ) 
Im(u/*u>~)  «c 


(2.31) 


For  a  wave  of  harmonic  form,  this  is  equivalent  to  the 
ratio  of  cospectrum  and  quadrature  power.  In  the 
three  simple  examples  above,  the  respective  nondi¬ 
mensional  contributions  are  (2 mi?)"*1,  (2 m/W 
and  2 uz/N.  In  (2.21)  the  nondimensional  value  is 
uz/2N. 

From  the  variety  of  cases  analyzed  here,  it  is  obvi¬ 
ous  that  a  horizontal  buoyancy  flux  can  arise  at  0(1) 
from  wave  reflection,  or  at  0(p)  within  a  region  of 
basic  state  variation,  wave  transience,  saturation, 
or  damping.  The  0(/i)  contribution  should  not  be 
downplayed,  particularly  if  it  is  due  to  critical-level 
interaction  within  a  zone  of  strong  vertical  shear.  In 
this  case,  zonal-wind  perturbations  grow  as  m1^2  up 
to  the  breaking  level.  The  observed  horizontal  buoy¬ 
ancy  flux  might  actually  be  larger  within,  rather 
than  below,  the  shear  zone  on  account  of  the  fac¬ 
tor  of  c  appearing  in  the  denominator  of  (2.21),  or 
because  of  wave  saturation  (2.28).  Saturation,  how¬ 
ever,  would  rapidly  diminish  its  magnitude  above 
the  breaking  level. 

The  harmonic  form  has  been  used  for  convenience, 
but  is  unnecessary.  Starting  from  (2.4)  the  second 
term  on  the  rhs  of  (2.10)  may  be  expanded  as  fol¬ 
lows: 


-poCV  —  poC2uz  ~  PoCu£ 

=  p0^2-poCW  _ 

=  PoC*Uz  -  Po(¥C)t  + 

-  poc yz  -  p0(Fc)t  -  poWx 

=  PoC2Uz  -  Po(f'C')t 

+J-(/W,)(PoC%  (2.32) 

Po 

Using  (2.12a)  it  can  be  shown  that 

-2poCV  -  PoC2uz  -  Po(F^)t  +  -Tj-  (pfflC) 

Po  Oz 

(2.33) 

Quadratic  terms  on  the  rhs  are  evaluated  from  the 
lowest  order  wave  solutions.  The  last  term  con- 
tains  the  vertical  momentum  flux  (proportional  to 
PoNip for  steady,  conservative,  slowly  varying 
waves)  and  a  basic  state  variation  (po/N)  inside 
the  vertical  derivative.  Substitution  of  the  harmonic 
form  in  this  equation  recovers  (2.13). 

c.  Numerical  solutions 

The  linear  equations  (2.1a-c)  may  be  combined  to 
give 
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V*sz  4-  mV  =  0 


(2.34) 


for  a  single  harmonic,  where 


N 2  uzz 
(c  —  u)2  c  — fi 


(2.35) 


for  hydrostatic,  incompressible,  steady,  conserva¬ 
tive  waves.  The  nonhydrostatic  and  quasicompress- 
ible  versions  may  be  found  in  Booker  and  Brether- 
ton  (1967)  and  Holton  (1975),  respectively.  If  me¬ 
chanical  and  thermal  damping  coefficients  are  equal 
and  independent  of  height,  then  c  — ►  c  +  zco  where 
Co  =  a/k .  To  illustrate  the  effect  of  damping,  a 
small  imaginary  part  may  be  added  to  the  phase 
speed.  Eq.  (2.32)  was  solved  using  the  tridiagonal 
algorithm  (with  801  equally  spaced  gridpoints)  sub¬ 
ject  to  boundary  conditions 


il>(zB)  =  1  (2.36a) 

(zt)  =  im(zr)w(zT)  (2.36b) 

choosing  the  sign  of  m  in  (2.34b)  to  correspond  to  an 
upward  propagating  wave.  Mean  flow  profiles  were 
of  the  form 


u  =  uq  4*  (fir  -  «b)  • 


1  +  tanh 


(2.37) 


and  N  =  .02  s  l.  In  the  following,  zq  =  0,  z\f  = 
7.5  km,  and  zt  =  15  km.  For  case  1  (weak  shear): 


ub  —  0 
ut  =  10  ms'1 
D  =  2  ms-1  km'1 

c  =  15  ms-1.  (2.38) 

For  case  2  (strong  shear): 

Uq  =  —20  ms'1 
ut  =  20  ms'1 
D  =  2  ms^km'1 

c  —  25  ms'1.  (2.39) 

These  profiles  are  shown  in  Figs.  la.  lc.  Case 
2  is  representative  of  the  QBO  near  25  km,  with 
maximum  shear  10  ms'1  km"”1,  corresponding  to  a 
Richardson  number  Ri  =  4. 

The  horizontal  buoyancy  flux  is  shown  in  Figs,  lb. 
Id  for  cases  1,  2  respectively,  along  with  (vertically 
constant)  momentum  flux.  The  WKB  approxima¬ 
tion  (2.21)  is  also  shown  (dashed  line).  In  case  1, 
there  is  a  small  buoyancy  flux  in  the  shear  zone,  and 
none  outside.  Nor  is  there  any  evidence  of  reflection. 
The  WKB  solution  is,  for  practical  purposes,  exact. 
The  momentum  flux  is  .00133366,  almost  identical 
to  the  no-shear  value  (m  =  .0013).  In  case  2.  there 
is  a  larger  buoyancy  flux  in  and  beneath  the  shear 
zone,  with  wavy  structure  below  due  to  interference 
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of  upward  and  downward  waves  in  this  region.  The 
wavelength  of  the  pattern  is  approximately  half  that 
of  the  wave  itself  (c/.  2.17).  The  WKB  solution 
omits  the  small  reflected  component  but  is  other¬ 
wise  accurate.  The  momentum  flux  is  .00039424. 
slightly  less  than  the  no-shear  value  (m  =  .00044). 

Two  variations  of  case  2  were  examined.  In  the 
first,  a  small  damping  was  added,  Co  =  1  ms^1.  The 
momentum  flux  now  decays  with  height,  markedly 
so  as  c  diminishes  entering  the  shear  zone  (not 
shown).  The  buoyancy  flux  is  slightly  more  posi¬ 
tive  than  before;  the  only  major  difference  is  that  it 
is  relatively  larger  when  measured  against  the  mo¬ 
mentum  flux  in  the  shear  zone.  In  the  second  vari¬ 
ation  (shown  in  Figs,  le,  If),  Z?  =  0.5  ms~lkm~l; 
the  shear  is  four  times  that  of  case  2  (Ri  — *  1/4), 
with  other  parameters  the  same.  The  reflected 
wave  is  larger,  and  the  momentum  flux  smaller 
(.000148514),  so  that  the  buoyancy  flux  is  relatively 
much  larger  than  in  Fig.  Id.  The  ratio  of  horizontal 
buoyancy  flux  to  vertical  flux  of  horizontal  momen¬ 
tum  has  no  limit,  since  the  momentum  flux  vanishes 
with  perfect  reflection.  The  amount  of  reflection  ob¬ 
served  in  Fig.  If  is  close  to  that  predicted  by  Ache- 
son’s  (1976)  piecewise-continuous  model  (R  =  0.8). 

Other  cases  using  the  shear  of  Fig.  le  were  ana¬ 
lyzed  for  larger  phase  speeds  (|c|  >50  ms^1)  that 
might  be  associated  with  fast  Kelvin  or  gravity 
waves  and  tides.  Buoyancy  fluxes  are  sharply 
peaked  in  the  shear  zone,  with  less  evidence  of  re¬ 
flection.  When  mean  flows  are  reversed,  the  buoy¬ 
ancy  fluxes  reverse  sign,  but  are  similar  in  magni¬ 
tude.  When  phase  speeds  are  reversed,  the  buoy¬ 
ancy  fluxes  do  not  change  sign  and  are  again  similar 
in  magnitude. 

d.  Comments 

In  order  to  apply  these  examples  to  the  atmo¬ 
sphere,  wave  amplitudes  must  be  weighted  by  a  fac¬ 
tor  exp(z/2H)  and  corrections  due  to  com¬ 
pressibility  should  be  retained.  If  the  major  ef¬ 
fect  is  wave  growth,  the  cospectrum  and  quadra¬ 
ture  power  corresponding  to  that  shown  in  the  fig¬ 
ures  will  have  an  exponential  envelope,  growing  with 
height.  Moreover,  observable  fields  like  u\  V  will 
grow  relative  to  ty* ,  as  c  — *  0,  by  a  factor  c~l .  These 
two  effects  will  tend  to  emphasize  features  in  the 
shear  zone  over  those  below.  Apart  from  reflection, 
ulTr  will  maximize  between  the  base  of  the  shear 
zone  and  some  point  just  above  the  breaking  level  (if 
one  exists).  A  more  detailed  analysis  of  observations 
of  Kelvin  and  gravity  waves  in  relation  to  the  de¬ 
scending  QBO  shear  zones  is  warranted  before  draw¬ 
ing  specific  conclusions.  Concerning  Kelvin  waves, 
it  should  be  noted  that  their  latitudinal  scale  con¬ 
tracts  approaching  the  critical  layer,  causing  wave 
amplitude  to  pile  up  near  the  equator  more  dramat¬ 
ically  than  in  the  two-dimensional  problem.  This  is 


a  WKB  effect  (with  slowly  varying  approximation  in 
time  and  height,  but  not  latitude,  as  in  Boyd,  1987) 
and  does  not  require  vertical  reflection  from  slop¬ 
ing  surfaces  Qj  =  /,  although  such  reflection  may  be 
important  if  our  two-dimensional  examples  are  rep¬ 
resentative.  Mean  winds  are  sometimes  asymmet¬ 
ric  about  the  equator,  so  maximum  wave  amplitude 
may  be  shifted  slightly  off  the  equator  as  a  result  of 
reflection. 

Time  series  of  cospectrum  power  observed  by  Sato 
et  al  (1994)  reveal  a  qualitative  difference  between 
long  and  short-period  disturbances  at  the  equator. 
The  cospectrum  and  quadrature  power  of  ~5-20  day 
disturbances  are  positively  enhanced  at  the  onset 
of  QBO  westerlies,  but  are  much  smaller  at  other 
times.  By  contrast,  the  cospectrum  power  of  short- 
period  disturbances  oscillates  with  the  QBO  be¬ 
tween  small  positive  and  negative  values  of  nearly 
equal  magnitude.  Although  the  quadrature  power 
at  short  periods  is  markedly  less  than  that  of  long- 
period  waves,  the  relatively  large  magnitude  and  os¬ 
cillatory  behavior  of  the  cospectrum  C(u',T ')  sug¬ 
gests  that  it  may  represent  a  real  phenomenon.  If 
so,  two  possible  interpretations  are: 

1)  Fast  Kelvin  waves,  although  transparent  to 
QBO  westerlies  (Hirota,  1978;  Salby  et  a/.,  1984; 
Maruyama,  1991),  might  contribute  to  a  small  hor¬ 
izontal  buoyancy  flux  as  they  encounter  QBO  shear 
zones.  In  order  for  C(u\T')  to  be  enhanced  apart 
from  Q (ti',  T"),  the  vertical  shear  must  be  extremely 
strong  (e.g.,  Fig.  If). 

2)  The  gravity  wave  spectrum  contains  eastward  and 
westward  propagating  waves  (Pfister  et  a/.,  1993). 2 
If  waves  of  equal  but  opposite  phase  speed  are  in¬ 
coherent,  they  will  contribute  alike  to  the  horizon¬ 
tal  buoyancy  flux  (2.21),  but  momentum  fluxes  will 
tend  to  cancel.  The  vertical  shear  need  not  be 
strong. 

A  possible  problem  with  either  of  these  explana¬ 
tions  is  that,  unless  phase  speeds  are  large  compared 
to  the  QBO,  waves  will  be  selectively  amplified  (as 
measured  by  v!)  in  the  QBO  shear  zones  as  their 
intrinsic  phase  speed  is  reduced.  For  these  waves,  it 
is  difficult  to  see  how  C(u',r')  could  be  amplified 
without  a  similar  (though  possibly  smaller)  amplifi¬ 
cation  of  Q(u',Tf).  (Downward  propagating  waves 
would  be  attenuated  by  damping  or  saturation  at 
higher  levels,  unless  phase  speeds  are  large.)  Con¬ 
sistent  with  the  role  of  vertical  shear  is  the  observa¬ 
tion  in  Fig.  11  of  Sato  et  al  (1994)  that  short-period 
C{u\T)  is  amplified  at  slightly  different  altitudes 
depending  on  the  phase  of  the  QBO. 

Given  the  likelihood  that  many  gravity  waves  are 
present  at  the  same  time,  it  may  be  worthwhile 
to  construct  spectral  analogues  of  the  monochro- 

2  The  role  of  gravity  waves  in  the  QBO  was  emphasized  by 

Dunkerton  (1991)  and  Takahashi  and  Bovilie  (1992)  and 

is  part  of  the  original  Lindzen  and  Holton  (1968)  theory. 
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matic  formulae  for  iTT' .  One  might  assume,  for 
example,  that  individual  waves  saturate  just  be¬ 
low  critical  levels,  but  that  the  resulting  turbulence 
causes  mechanical  and  thermal  damping  of  other 
(non-breaking)  waves.  Presumably  the  waves  most 
susceptible  to  damping  in  the  westerly  acceleration 
phase  of  the  QBO  would  be  those  with  phase  speeds 
of  positive  sign  and  vice  versa.  Simultaneous  obser¬ 
vations  of  wind  and  temperature  in  the  equatorial 
lower  stratosphere,  with  better  temporal  resolution 
(e.g.,  radar,  lidar),  would  provide  valuable  informa¬ 
tion  with  which  to  evaluate  these  effects  in  the  con¬ 
text  of  a  gravity  wave  spectrum. 

3.  Conclusion 

The  equations  of  motion  describing  two-di¬ 
mensional  internal  gravity  waves  were  analyzed  to 
derive  an  expression  for  the  horizontal  buoyancy 
flux — that  is,  the  correlation  between  buoyancy  (or 
temperature)  and  horizontal  velocity — in  various 
cases  involving  vertical  shear.  It  was  shown  that,  al¬ 
though  this  correlation  vanishes  at  zeroth  order  (in 
fi)  for  a  single  wave,  its  first  order  contribution  is 
nonzero  due  to  shear,  even  for  steady,  conservative, 
incompressible  waves.  Departures  from  steady,  con¬ 
servative  motion  or  quasi-compressibility  also  cause 
a  nonzero  correlation.  Several  cases  were  analyzed 
and  some  numerical  results  obtained  for  waves  ap¬ 
proaching  a  critical  layer  of  reduced  intrinsic  phase 
speed.  With  weak  shear,  the  buoyancy  flux  is  small 
relative  to  momentum  flux,  as  expected  from  the 
perturbation  theory.  Strong  vertical  shear  enhances 
the  buoyancy  flux  within  the  shear  zone  and  causes 
partial  reflection  beneath,  producing  a  nonzero  0(1) 
correlation  in  this  region. 

These  effects  may  explain  recent  observations  of 
zoned  wind  and  temperature  cospectra  in  the  equa¬ 
torial  lower  stratosphere  (Sato  et  al. ,  1994).  Such 
observations,  at  least,  do  not  immediately  preclude 
their  interpretation  as  Kelvin  or  gravity  waves.  Tak¬ 
ing  the  analysis  of  this  paper  into  account,  the  be¬ 
havior  of  cospectrum  and  quadrature  power  relative 
to  the  QBO  seems  qualitatively  consistent  with  such 
an  interpretation  for  long  and  short-period  distur¬ 
bances,  respectively.  For  long-period  disturbances, 
equatorial  C(u',Tf)  and  Q{u,tf)  are  selectively  en¬ 
hanced  during  the  onset  of  QBO  westerly  phases, 
but  otherwise  small,  as  expected  for  Kelvin  waves. 
For  short-period  disturbances,  C{u*,Tf )  alternates 
sign  with  the  QBO,  as  expected  of  a  gravity  wave 
spectrum  containing  both  eastward  and  westward 
propagating  components  (Pfister  et  al. ,  1993). 

Although  the  theoretical  results  suggest  that  the 
cospectrum  will  be  measurable,  particularly  in  a 
critical  layer,  further  quantitative  analysis  of  the  ra¬ 
tio  C:Q  is  desirable  as  it  may  reveal  useful  informa¬ 
tion  about  the  anisotropy  of  the  spectrum  or  the  na¬ 
ture  of  wave  saturation  and  absorption.  Simultane¬ 


ous  radar  and  lidar  measurements  in  the  equatorial 
lower  stratosphere  would  be  particularly  valuable.  A 
three-dimensional  analysis  of  Kelvin  waves  in  verti¬ 
cal  shear  and  their  associated  horizontal  buoyancy 
flux  would  significantly  extend  these  results.  Such 
topics  are  beyond  the  scope  of  this  paper  but  fruit¬ 
ful  avenues  for  exploration.  A  final  point  is  that, 
although  the  horizontal  buoyancy  (or  temperature) 
flux  has  been  emphasized  for  its  interpretive  value — 
it  is,  after  all,  exactly  nondivergent  since  indepen¬ 
dent  of  x — real  wavepackets  are  of  finite  horizontal 
extent  and  may  therefore  contribute  to  a  net  hori¬ 
zontal  transport  under  conditions  more  general  than 
those  considered  here. 
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ABSTRACT 

The  excitation  and  propagation  of  inertia-gravity  waves  (IGWs)  generated  by  an  unstable  baroclinic  wave 
was  examined  with  a  high-resolution  3D  nonlinear  numerical  model.  IGWs  arose  spontaneously  as  the  tropo¬ 
spheric  jetstream  was  distorted  by  baroclinic  instability  and  strong  parcel  accelerations  took  place,  primarily  in 
the  jetstream  exit  region  of  the  upper  troposphere.  Subsequent  propagation  of  IGWs  occurred  in  regions  of  strong 
windspeed — in  the  tropospheric  and  stratospheric  jets,  and  in  a  cutoff  low  formed  during  the  baroclinic  lifecycle. 

IGWs  on  the  flanks  of  these  jets  were  rotated  inward  by  differential  advection  and  subsequently  absorbed  by  the 
model’s  hyperdiffusion.  Although  absorption  of  IGWs  at  the  sidewalls  of  the  jet  is  an  artifact  of  the  model,  1GW 
propagation  was  for  the  most  pan  confined  to  regions  with  an  intrinsic  period  shorter  than  the  local  inertial 
period.  Only  a  few  IGWs  were  able  to  penetrate  the  middle  stratosphere,  due  to  weak  winds  or  an  unfavorable 
alignment  of  wavevector  with  respect  to  the  mean  flow. 

IGWs  are  important  both  as  a  synoptic  signal  in  the  jetstream,  which  may  influence  subsequent  tropospheric 
developments,  and  as  a  source  of  isen tropic  or  cross-isentropic  mixing  in  the  lower  stratosphere.  The  authors’ 
results  demonstrated  for  the  first  time  numerically  a  significant  isentropic  displacement  of  potential  vorticity 
isopleths  due  to  IGWs  above  the  tropopause.  Since  conditions  for  IGW  propagation  are  favorable  within  a  jet, 
a  region  of  strong  isentropic  potential  vorticity  gradient,  it  is  likely  that  inertia -gravity  waves  affect  the  per¬ 
meability  of  the  lower  stratospheric  vortex  and  may  in  some  instances  lead  to  stratosphere- troposphere  exchange. 


1.  Introduction 

Inertia-gravity  waves  (IGWs)  are  ubiquitous 
within  the  atmosphere,  especially  the  upper  tropo¬ 
sphere  and  lower  stratosphere.  These  motions  are 
characterized  by  short  vertical  wavelengths  (  —  1-5 
km)  and  small  intrinsic  frequencies  Cj  near  the  Cor¬ 
iolis  frequency  /(Barat  1 983;  Yamanaka  and  Tanaka 
1984;  Thomas  et  al.  1992:  Sato  1994).  Their  hori¬ 
zontal  scale,  which  is  usually  estimated  rather  than 
measured,  is  thought  to  be  large  (-200-1000  km) 
compared  to  mesoscale  gravity  waves.  While  mo¬ 
mentum  transport  is  ascribed  to  mesoscale  waves 
(Fritts  1984),  it  is  the  role  of  IGWs  in  constituent 
mixing,  and  not  momentum  transport,  that  is  proba¬ 
bly  more  relevant.  IGW  mixing  may  be  isentropic 
(Pierce  and  Fairlie  1993;  Led  well  et  al.  1993),  like 
that  of  synoptic-scale  baroclinic  waves,  or  cross  is¬ 
entropic  due  to  instabilities  arising  within  the  wave 
field  (Dunkerton  1984;  Fritts  and  Rastogi  1985). 
Notwithstanding  the  importance  of  IGW  propaga¬ 
tion,  breaking,  and  absorption,  one  would  also  like 
to  know  something  about  the  excitation  of  these 
waves  and  whether  their  geographical  and  climato- 
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logical  variability  can  be  understood  in  terms  of  the 
postulated  forcing  mechanisms. 

The  origin  of  inertia-gravity  waves  is  ambiguous 
because  regions  of  strong  mean  flow  are  a  favored  locus 
of  IGW  excitation  and  propagation  (Dunkerton  1984). 
IGWs  observed  within  or  above  the  tropospheric  jet¬ 
stream  were  not  necessarily  excited  there:  due  to  their 
slow  propagation,  the  source  may  have  been  as  far 
away  as  the  Tropics  (Ushimaru  and  Tanaka  1990). 
Case  studies  often  suggest  a  connection  between  IGWs 
and  nearby  synoptic  weather  systems  (Sato  1 989, 1993, 
1 994;  Chan  et  al.  1 99 1 ;  Thomas  et  al.  1 992;  Eckermann 
and  Vincent  1993).  Therefore  it  is  likely  that  extra- 
tropical  sources  dominate.  That  gravity  wave  variance 
increases  above  oceanic  frontal  systems  and  jetstreams 
(Fritts  and  Nastrom  1992)  indicates  that  sources  other 
than  topography  are  important — an  essential  observa¬ 
tion  for  understanding  gravity  wave  transport  in  the 
Southern  Hemisphere  and  Tropics.  Part  of  the  increase 
is  undoubtedly  due  to  high-frequency  waves  launched 
by  moist  convection  or  shear  instabilities  within  frontal 
zones.  (It  was  primarily  mesoscale  variance  that  was 
observed  by  Fritts  and  Nastrom.)  However,  one  also 
expects  a  low-frequency  IGW  component  due  to  the 
unstable  baroclinic  wave  that  would  exist  in  the  ab¬ 
sence  of  convection,  arising  from  the  transiency  of  tro¬ 
pospheric  flow  as  it  relaxes  to  a  balanced  state.  In  a 
rotating  fluid  this  process  is  referred  to  as  geostrophic 
adjustment  (Blumen  1972;  Luo  and  Fritts  1993,  and 
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references  therein).  During  the  adjustment,  gravita¬ 
tional  oscillations  emanate  and  disperse  away  from  the 
source.  Geostrophic  adjustment  was  interpreted  by  Val- 
lis  ( 1992)  as  a  process  that  minimizes  the  energy  for  a 
fixed  distribution  of  potential  vorticity  (PV).  This  in¬ 
terpretation  assumes  that  the  adjustment  is  rapid  com¬ 
pared  to  the  evolution  of  PV.  General  methods  for  find¬ 
ing  energy  extrema  that  allow  rearrangements  of  PV  as 
a  materially  conserved  quantity  were  discussed  by  Val- 
lis  et  al.  (1989)  and  Shepherd  (1990). 

Many  studies  examined  geostrophic  adjustment  in  an 
idealized  way,  as  linear  relaxation  from  an  initial  im¬ 
balance  of  mass  or  momentum,  with  simplified  (Car¬ 
tesian  or  axisymmetric)  geometry.  In  real  life  the  ad¬ 
justment  involves  complicated,  three-dimensional  mo¬ 
tion  with  short  vertical  and  (ultimately)  horizontal 
scales.  The  need  for  consistent  vertical  resolution  in 
primitive  equation  models  was  stressed  by  Lindzen  and 
Fox-Rabinovitz  (1989),  and  their  comments  pertain 
equally  well  to  geostrophic  adjustment.  Only  within  the 
last  few  years  have  computational  resources  become 
available  to  make  a  realistic  modeling  study  of  cyclo- 
strophic  or  geostrophic  adjustment  practical  (e.g.,  Le- 
long  and  McWilliams  1993).  One  rather  fundamental 
point  that  has  not  been  extensively  addressed  is  the 
mechanism  by  which  an  initially  balanced  flow  evolves 
into  an  unbalanced  state,  followed  by  geostrophic  ad¬ 
justment  (Van  Tuyl  and  Young  1982;  Duffy  1990).  It 
is  expected  that  inertia-gravity  waves  arise  spontane¬ 
ously  in  primitive  equation  models  (Warn  and  Menard 
1986).  Since  IGWs  are  problematic  to  weather  predic¬ 
tion,  the  emphasis  has  been  mainly  on  how  to  eliminate 
them,  rather  than  explicit  modeling  of  IGW  propaga¬ 
tion  and  absorption.  This  omission  is  surprising  in  view 
of  the  large  body  of  literature  on  ageostrophic  circu¬ 
lations  within  synoptic  weather  systems  (e.g.,  Cammas 
and  Ramond  1989;  Keyser  et  al.  1989,  and  references 
therein)  and  their  role  in  frontogenesis,  tropopause 
folding,  and  stratosphere-troposphere  exchange  (Sha¬ 
piro  1 980;  Keyser  and  Shapiro  1986;  Shapiro  and  Key¬ 
ser  1990).  By  analogy  to  tidal  theory,  the  ageostrophic 
motion  contains  a  “forced”  evanescent  component 
(the  meridional  circulation)  and  a  “free”  radiating 
component  (the  inertia-gravity  waves). 

Recently  the  authors  performed  numerical  simula¬ 
tions  to  examine  the  evolution  of  IGWs  in  a  lifecycle 
of  midlatitude  baroclinic  instability.  Patterned  after  the 
classic  experiments  of  Simmons  and  Hoskins  (1978, 
1980),  our  simulations  described  the  nonlinear  lifecy¬ 
cle  of  an  unstable  baroclinic  wave  assuming  wavenum¬ 
ber  6  symmetry  on  the  sphere  but  incorporated  much 
higher  vertical  resolution  and  a  domain  extending  into 
the  lower  stratosphere  in  order  to  explicitly  model  the 
evolving  inertia-gravity  waves  and  their  radiation  into 
the  middle  atmosphere.  This  work  was  complementary 
to  several  other  studies  (beginning  with  Ley  and  Peltier 
1978;  see  also  Uccelini  and  Koch  1987;  Gall  et  al. 
1988;  Gamer  1989;  Benard  et  al.  1992a,b;  Blumen 


1992;  Snyder  et  al.  1993)  concerned  with  gravity  wave 
excitation  in  relation  to  frontogenesis.  We  attempted  to 
simulate  inertia-gravity  waves  excited  in  the  upper  tro¬ 
posphere  by  an  evolving  baroclinic  instability,  consid¬ 
ering  their  vertical  as  well  as  horizontal  propagation. 
Numerical  results  demonstrated  the  role  of  baroclinic 
instability  in  exciting  inertia-gravity  waves,  due  to 
rapid  deformation  of  upper-tropospheric  flow  and 
strong  parcel  accelerations. 

Although  our  integrations  were  designed  to  yield 
well-resolved  inertia-gravity  waves  independent  of 
model  resolution,  only  the  excitation  and  initial  prop¬ 
agation  of  IGWs  could  be  realistically  simulated.  Their 
asymptotic  propagation  and  breakdown  requires  verti¬ 
cal  and  horizontal  resolution  well  beyond  that  used  here 
(e.g.,  Dunkerton  and  Robins  1992).  Resolved  patterns 
of  horizontal  velocity  divergence  were  nonetheless 
quite  complex,  containing  an  assortment  of  gravity 
waves  and  IGWs  due  to  surface  frontogenesis  as  well 
as  parcel  accelerations  in  the  upper- level  jet.  This  paper 
emphasizes  geostrophic  adjustment  as  the  numerically 
best  resolved,  and  probably  most  fundamental,  mech¬ 
anism  of  IGW  excitation.  Section  2  describes  the  nu¬ 
merical  model.  Section  3  discusses  the  baroclinic  life- 
cycle  and  the  resulting  IGWs,  examining  IGW  wave 
parameters  and  Lagrangian  parcel  accelerations.  Some 
IGW-related  issues  are  discussed  briefly  in  section  4. 

2.  Model  description 

Simulations  were  performed  with  a  slightly  modified 
version  of  the  3D,  global,  hydrostatic  primitive  equa¬ 
tion  model  of  Young  and  Villere  (1985).  This  is  a 
spectral  transform  model  with  finite  differencing  in  the 
vertical.  There  is  no  topography  or  moisture.  Subgrid- 
scale  processes  are  represented  by  sixth-order  horizon¬ 
tal  hyperdiffusion  acting  on  the  vorticity,  diver¬ 
gence,  and  temperature  fields.  The  spectral  truncation 
in  runs  described  here,  unless  stated  otherwise,  was  tri¬ 
angular  at  total  horizontal  wavenumber  126  (denoted 
T126) — approximately  equivalent  to  1°  of  horizontal 
resolution.  Hyperdiffusion  is  needed  to  prevent  non¬ 
physical  behavior  as  sharp  horizontal  gradients  are  gen¬ 
erated  at  the  smallest  resolvable  scales.  In  practice,  vt 
corresponding  to  an  e-folding  time  (r)  of  1  h  or  less 
for  the  smallest  scale  is  sufficient  for  stable  integra¬ 
tions.  At  T126,  v6  =  1.36  x  1025  m6s“1  was  used, 
corresponding  to  a  damping  time  of  about  20  min  for 
the  smallest  resolved  scale. 

The  vertical  coordinate  was  uniformly  spaced  in  log 
pressure  with  A z  -  700  m,  unless  stated  otherwise,  and 
the  domain  extended  vertically  from  z  -  0-35  km  (51 
levels).  Boundary  conditions  assumed  zero  geometric 
vertical  velocity  at  the  lowest  level,  while  log  pressure 
vertical  velocity  was  zero  at  the  uppermost  level.  There 
was  no  surface  friction. 

As  discussed  in  section  3d,  comparison  simulations 
were  performed  at  various  combinations  of  horizontal 
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Fig.  1 .  Axisymmetric  initial  state  used  in  model  simulations,  showing  isopleths  of  potential  temperature  (solid 
contours  ranging  from  270  to  640  K)  and  mean  zonal  wind  (shading  and  alternating  solid/dotted  contours). 


and  vertical  resolution,  for  example,  T42,  T84,  T126, 
and  T156,  and  Az  =  350  or  700  m.  Our  custom  has 
been  to  maintain  approximately  the  same  hyperdiffu¬ 
sion  damping  time  for  the  smallest  resolved  scale,  re¬ 
gardless  of  horizontal  truncation.  Higher  horizontal  res¬ 
olutions  therefore  imply  smaller  hyperdiffusion  coef¬ 
ficients  and  vice  versa,  unless  stated  otherwise. 

Our  primary  interest  was  the  generation  of  IGWs  by 
geostrophic  adjustment  of  the  tropospheric  jetstream  as 
it  is  distorted  by  a  developing  baroclinic  wave.  For  this 
purpose,  the  evolution  should  be  realistic  but  not  overly 
complicated.  Following  Simmons  and  Hoskins  (1978), 
the  model  was  initialized  with  a  small-amplitude  wave- 
number  6  normal  mode  superposed  on  a  zonally  uni¬ 
form,  baroclinically  unstable  flow.  This  procedure 
yields  within  several  days  of  model  integration  a  well- 
resolved  and  fairly  realistic  finite  perturbation  to  the 
jetstream,  followed  by  a  barotropic  decay  phase  as  de¬ 
scribed  in  section  3a.  There  are  of  course  six  identical 
lifecycles  around  the  Northern  Hemisphere. 

For  time  integration  the  model  used  a  semi-implicit 
scheme  (Hoskins  and  Simmons  1975)  with  time  step 
of  600  s,  as  required  by  the  CFL  condition  for  wind- 
speeds  attained  in  the  simulation.  As  we  were  mainly 
interested  in  the  inertia-gravity  waves,  and  it  is  known 
that  the  semi-implicit  scheme  affects  (high  frequency) 
gravity  wave  phase  speeds,  two  T84  simulations  of  the 
same  lifecycle  were  compared  using  either  the  semi- 


implicit  scheme  or  a  fully  explicit,  centered-difference 
time  integration  scheme.  These  simulations  used  time 
steps  of  1200  and  150  s,  respectively  (demonstrating 
the  computational  advantage  of  the  semi-implicit 
scheme).  Comparing  simulations  during  the  mature 
stage  of  the  lifecycle  showed  identical  lifecycles  oc¬ 
curring.  Development  of  the  lifecycle’s  mature  stage 
was  delayed  by  about  1  day  in  the  semi-implicit  case, 
probably  reflecting  the  different  treatment  of  initial 
condition  by  the  two  schemes.  Once  this  stage  was 
reached,  subsequent  evolution  did  not  differ  noticeably. 
In  particular,  the  low-frequency  inertia-gravity  waves 
of  interest  were  insensitive  to  the  choice  of  time  inte¬ 
gration  scheme  and  appeared  nearly  identical  in  the  two 
cases. 

3.  Baroclinic  lifecycle  and  IGW  excitation 

a.  Baroclinic  lifecycle  simulation 

As  described  above,  a  baroclinic  lifecycle  analogous 
to  that  of  Simmons  and  Hoskins  ( 1978)  was  simulated 
with  resolution  adequate  to  model  inertia-gravity 
waves.  The  unstable  baroclinic  wave  developed  on  an 
initially  axisymmetric  jet  centered  at  45°N  as  shown  in 
Fig.  1.  The  initial  state  contained  several  realistic  fea¬ 
tures,  namely,  a  change  of  static  stability  across  the 
tropopause  (N2  =  1.5  x  10-4  s-2  in  the  troposphere. 
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N2  =  3.5  X  1 0 “4  s’*2  in  the  stratosphere),  a  sloping 
tropopause,  a  stratospheric  polar  night  jet,  and  wind- 
speed  minimum  between  stratospheric  and  tropo¬ 
spheric  jetstreams.  The  initial  state,  with  its  strong  jet 
displaced  north  of  the  climatological  Jetstream,  might 
represent  a  cross  section  through  midoceanic  (Pacific 
or  Atlantic)  storm  tracks  during  winter.  It  was  dis¬ 
turbed  by  a  small  wavenumber  6  perturbation,  which 
subsequently  grew  via  baroclinic  instability  to  finite 
amplitude.  After  7  days  of  linear  growth  the  baroclinic 
wave  appeared  as  undulations  on  the  westerly  current. 

Discussion  of  the  baroclinic  wave  and  accompany¬ 
ing  IGWs  begins  on  day  9  of  the  model  run,  by  which 
time  the  baroclinic  wave  at  upper  levels  attained  finite 
amplitude.  At  midtropospheric  levels  an  amplifying 
ridge,  and  deepening  trough  and  low  pressure  center, 
show  that  the  wave  is  still  in  its  growth  phase.  Days  10 
and  1 1  may  be  considered  the  wave’s  mature  stage  as 
growth  ceases  and  the  elongated  SW-trailing  trough 
sheds  a  cutoff  low  near  30°N  (Fig.  2a) .  After  formation 
of  the  cutoff  low,  the  lifecycle  enters  its  decay  stage 
with  the  main  westerly  jetstream  shifting  northward  on 
day  12  to  50°~60°N  and  becoming  more  zonally  sym¬ 
metric  and  intense  (Fig.  2b).  At  this  time  the  baroclin- 
icity  has  been  largely  removed  from  midlatitudes  where 
areas  of  weak  horizontal  temperature  gradient  predom¬ 
inate,  flanked  to  the  north  and  south  by  frontal  zones. 

At  low  levels  the  characteristic  poleward  heat  flux 
of  baroclinic  instability  is  seen  as  the  cold  air,  is  ad- 
vected  southward  and  the  warm  air  northward  (Figs. 
2c,d).  Cyclogenesis  and  frontogenesis  proceed  rapidly 
with  a  surface  low  and  warm/cold  fronts  well  formed 
by  day  9.  By  this  time,  the  surface  fronts  have  sharp¬ 
ened  to  the  point  that  their  width  is  limited  by  the  mod¬ 
el’s  small-scale  hyperdiffusion. 

The  low-level  cyclogenesis  and  frontogenesis 
strongly  resemble  the  observed  development  of  syn¬ 
optic  systems  as  discussed,  for  example,  by  Shapiro 
and  Keyser  (1990).  The  near-surface  temperature  on 
day  9  or  10  shows  the  so-called  T-bone  frontal  stage 
of  the  baroclinic  lifecycle,  a  sharp  front  at  60°-65°N 
extending  back  into  the  low  pressure  center,  while  the 
cold  front,  advancing  eastward  of  the  low,  forms  the 
vertical  stem  of  the  T.  A  pocket  of  warm  air  has  been 
secluded  behind  the  cold  front  near  55°N,  15°W  on  day 
9  and  1 0,  as  has  been  observed  over  oceans  when  cy¬ 
clones  reach  the  mature  stage  of  their  lifecycle  (appar¬ 
ently  because  of  the  lower  surface  friction  there). 

The  baroclinic  lifecycle  also  resembles  that  of  Sim¬ 
mons  and  Hoskins  (1980)  and  Thomcroft  and  Hoskins 
(1990),  even  showing  the  secondary  frontal  develop¬ 
ment  discussed  by  the  latter  authors.  This  occurs  on 
days  1 1  and  12  as  the  upper-level  circulation  associated 
with  the  cutoff  low  drifts  over  the  surface  cold  front, 
initiating  a  small-scale  frontal  wave  near  40°N,  20°W 
(Fig.  2d) — an  example  of  Peterssen’s  “type  B”  cy¬ 
clogenesis. 


In  the  tropopause  region  the  evolution  is  like  that 
described  for  the  midtroposphere.  The  jetstream  is  sim¬ 
ilarly  distorted  as  the  baroclinic  wave  radiates  upward. 
After  a  period  of  linear  growth,  the  wave  saturates  and 
a  cutoff  low  forms  by  day  1 1 .  This  is  a  period  of  rapid 
geostrophic  adjustment  throughout  the  troposphere 
but,  especially,  at  upper  levels  where  the  jetstream  is 
strongest. 

We  next  examine  evidence  of  gravity  wave  activity 
associated  with  surface  frontogenesis  and  geostrophic 
adjustment  in  the  upper  troposphere.  The  field  of  hor¬ 
izontal  velocity  divergence  is  useful  for  isolating  grav¬ 
ity  waves  from  the  large-scale  quasigeostrophic  flow 
because  of  the  divergent  nature  of  such  waves.  (By 
contrast,  the  signature  of  IGWs  in  geopotential  is 
hardly  visible.)  The  horizontal  divergence  field  in¬ 
cludes  forced  ageostrophic  circulations  accompanying 
the  varying  jetstream.  As  shown  next,  waves  located 
away  from  this  apparent  source  region  show  a  radiating 
pattern  of  divergence  distinct  from  that  associated  with 
the  forced  circulations. 

b.  IGW  generation  during  lifecycle 

The  divergent  component  of  tropospheric  flow 
evolves  quickly  and  in  a  complicated  manner  during 
the  baroclinic  lifecycle.  Analysis  of  horizontal  velocity 
divergence  (hereafter  V-  \H)  shows  gravity  wave  ac¬ 
tivity  associated  with  the  surface  fronts  and  cross-fron¬ 
tal  ageostrophic  circulations.  These  surface-related  fea¬ 
tures  weaken  with  increasing  height  because  fronto¬ 
genesis  is  strongest  at  ground  level,  and  the  radiating 
waves  are  dissipated  as  the  mean  flow  turns  with 
height.  At  the  tropopause  level  and  above  (10-25  km) , 
gravity  wave  activity  increases  dramatically  during  the 
lifecycle’s  mature  stage.  Figures  3a-d  show  V*  \H  at 
130  mb  (near  14-km  height)  on  days  9-12.  The  geo¬ 
potential  field  (Figs.  4a-d)  shows  the  location  of  the 
jetstream  at  this  level.  The  first  appearance  of  signifi¬ 
cant  gravity  wave  activity  in  the  upper  troposphere  oc¬ 
curs  on  day  9.  Over  the  next  day  this  activity  quickly 
increases  as  a  gravity  wavepacket  reaches  1 30  mb.  The 
wavepacket  is  strongest  on  day  11,  having  advanced 
eastward  with  a  secondary  wave  train  splitting  and 
propagating  southward  around  the  west  side  of  the  cut¬ 
off  low.  By  day  12,  the  main  packet  has  advanced  rap¬ 
idly  eastward  due  to  acceleration  of  the  jet.  A  third 
wave  train  appears  on  the  east  side  of  the  cutoff  low, 
unrelated  to  the  first  packet.  Over  subsequent  days, 
fragments  of  V  •  \H  become  more  widely  spread  and 
less  coherent  (not  shown). 

Figures  3  and  4  show  that  the  upper-level  gravity 
waves  tend  to  be  confined  to  the  vicinity  of  the  jet¬ 
stream,  propagating  with  wavevector  parallel  to  the  jet 
axis  at  midstream  but  at  an  angle  to  the  jet  along  its 
flanks.  This  is  also  observed  on  the  east  and  west  sides 
of  the  cutoff  low.  Interestingly,  at  late  times  the  rem¬ 
nants  of  IGWs  in  the  cutoff  low  are  carried  back  into 
the  jet  as  the  low  itself  is  absorbed  (not  shown). 
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Fig.  2.  (a,  b)  Geopotential  at  503  mb,  and 
on  days  10  (a,  c)  and  12  (b,  d)  of  the  TI26 
following  figures,  two  identical  60°  sectors  ai 


Fig.  3.  Horizontal  velocity  divergence  at  130  mb  on  consecutive  days  of 


Fig.  4.  As  in  Fig.  3  but  for  130-mb  geopotential.  Vectors  show  the 
horizontal  wind,  the  longest  vector  corresponding  to  62  m  s~*  in  (c). 
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While  model  gravity  waves  in  the  mid-  and  lower 
troposphere  can  be  related  to  frontal  activity,  those  in 
the  upper  troposphere  and  lower  stratospheric  generally 
cannot  and  appear  to  have  a  different  origin.  In  order 
to  provide  a  3D  view  of  these  waves,  and  to  clarify 
their  origin,  cross-sectional  views  of  the  wavepacket 
are  presented  in  Figs.  5  and  6.  Figures  5a,b  show  lati¬ 
tude-height  sections  of  V-  v#  on  days  10  and  11  at 
longitudes  that  intersect  the  gravity  wavepacket  on 
each  day.  In  these  cross  sections  the  dominant  feature 
is  the  enhanced  gravity  wave  activity  at  upper  levels, 
between  the  tropopause  and  ~20  km.  Frontal  features 
at  low  levels  are  also  evident,  largely  unconnected  to 
the  activity  aloft.  The  upper-level  waves  appear  like  a 
stack  of  pancakes,  with  phase  fronts  sloping  down  on 
either  side.  These  surfaces  also  tend  to  parallel  the  iso- 
pleths  of  constant  windspeed. 

Figures  6a-d  show  longitude-height  cross  sections 
ofV-v„  on  days  9-12  at  55°N,  a  latitude  intersecting 
the  main  wavepacket  on  those  days.  This  again  shows 
the  upper-level  wavepacket  to  be  the  dominant  feature. 
The  waves  have  phase  surfaces  sloping  uniformly  up¬ 
ward  to  the  west.  It  is  convenient  here  that  the  jetstream 
exit  region, 1  located  near  the  crest  of  the  geopotential 
ridge,  is  propagating  eastward  almost  parallel  to  a  lat¬ 
itude  circle  during  this  stage  of  the  lifecycle.  Thus  the 
longitude-height  cross  section  at  55°N  contains  the 
exit  region  on  these  days.  Figure  6  illustrates  an  im¬ 
portant  feature  of  the  upper-level  gravity  waves, 
namely  that  their  phase  and  group  speed  are  small  (or 
zero)  relative  to  the  jetstream  exit  region.  This  is  best 
seen  in  animations  of  the  cross-sectional  view,  but  it 
can  also  be  discerned  from  the  plots  for  day  9-11, 
where  the  first  three  distinct  wave  crests  (counting  up¬ 
ward  from  the  bottom  of  the  wavepacket)  are  advected 
eastward  with  the  jetstream  exit  region.  In  animation  it 
appears  that  the  jetstream  pattern  and  attendant  gravity 
waves  are  jointly  advected  eastward  as  a  coherent 
structure.  While  the  waves  are  approximately  station¬ 
ary  relative  to  the  jetstream  pattern,  they  are,  of  course, 
in  a  strong  westerly  flow,  so  their  intrinsic  phase  ve¬ 
locity  is  rapidly  westward.  The  intrinsic  group  velocity 
is  also  westward  as  discussed  below.  (Likewise,  the 
intrinsic  propagation  of  the  jetstream  pattern  is  west¬ 
ward  with  respect  to  the  fluid  in  the  upper  troposphere.) 
An  upward  group  velocity  can  be  discerned  for  these 
waves  over  the  3-day  period  of  Fig.  6,  as  the  wave- 
packet  expands  upward  with  time. 

The  impression  gained  from  Fig.  6  is  of  gravity  wave 
generation  in  the  upper  troposphere  near  the  level  of 
maximum  wind  in  the  vicinity  of  the  jetstream  exit  re¬ 
gion.  There  is  a  striking  vertical  asymmetry  with  most 


1  The  term  “exit  region”  is  used  to  describe  the  relatively  weak, 
tightly  curved  flow  downstream  of  the  jet  maximum.  Parcels  in  this 
area  do  not  decelerate  to  zero  but  are  deflected  southward  into  the 
west  side  of  the  developing  cutoff  low. 


of  the  IGWs  on  the  upper  side  of  the  jet,  where  static 
stability  increases  dramatically  with  height.  Over  the  3 
days  of  the  “mature  lifecycle”  stage  the  gravity  waves 
radiate  upward  into  the  lower  stratosphere  where,  ex¬ 
cept  for  one  or  two  small  packets,  most  of  them  ap¬ 
proach  a  critical  surface  and  are  absorbed.  (Note  that 
critical  layer  absorption  can  occur  when  the  mean  flow 
turns  with  height — the  windspeed  need  not  approach 
the  wave’s  phase  speed.)  There  seems  to  be  continual 
generation  of  gravity  wave  activity  near  the  level  of 
maximum  wind  during  this  time,  since  the  lower  edge 
of  the  wavepacket  does  not  move  upward.  Van  Tuyl 
and  Young  ( 1982)  noted  similar  stationarity  of  unbal¬ 
anced  motions  (IGW)  relative  to  the  jet  pattern  in  their 
two-level  mechanistic  experiments  of  geostrophic  ad¬ 
justment.  As  noted  by  those  authors,  IGWs  are  a  sub¬ 
synoptic  motion  that  can  have  significant  amplitude  in 
the  exit  region  of  strong  jetstreams.  In  our  high-reso¬ 
lution  simulations  IGWs  explain  a  substantial  part  of 
the  divergent  wind  component  and  at  maximum  am¬ 
plitude  (3-5  m  s-1)  constitute  about  10%  of  the  total 
wind  near  the  tropopause. 

To  test  the  conclusion  that  IGWs  are  due  to  geo¬ 
strophic  adjustment  of  the  upper  level  jet  and  are  not 
caused  by  gravity  wave  generation  near  the  surface 
frontal  regions,  the  simulation  was  repeated  with  strong 
damping  of  horizontal  velocity  divergence  in  the  lowest 
5  km.  Applying  additional  sixth-order  hyperdiffusion 
to  divergence  only — 500  times  stronger  at  the  surface 
than  the  standard  value  and  decreasing  linearly  to  zero 
at  5  km — gave  results  similar  to  the  standard  case  but 
with  much  smoother  horizontal  velocity  divergence  at 
low  levels,  effectively  eliminating  any  gravity  wave  ra¬ 
diation  from  the  surface  fronts.  At  upper  levels  the  flow 
was  very  similar  except  for  a  slight  reduction  of  jet  core 
speed  (~5%).  Most  importantly,  IGWs  were  again 
found,  unchanged  other  than  a  slight  reduction  of  am¬ 
plitude  (-10%). 

The  IGWs  are  well  resolved  in  this  T126  simulation, 
at  least  near  the  wavepacket  center  where  their  hori¬ 
zontal  wavelength  is  600- 1 000  km.  This  was  borne  out 
by  comparison  to  an  identical  T84  simulation  showing 
a  very  similar  IGW  packet  with  identical  phase  struc¬ 
ture  on  day  11  at  1 30  mb.  Differences  were  noticeable 
at  the  jet  flanks  where  the  IGW  horizontal  wavelength 
undergoes  contraction  as  discussed  in  sections  3c,d. 

c.  Description  of  IGWs  at  upper  levels 

The  gravity  waves  found  in  the  simulated  baroclinic 
lifecycle  are  of  inertia-gravity  type;  that  is,  N  >  |u>| 
^  1/ 1 ,  where  N  is  the  Brunt-Vaisala  frequency,  w  is 
the  intrinsic  wave  frequency,  and  /  is  the  Coriolis  pa¬ 
rameter.  Examining  the  center  of  the  main  IGW  packet 
between  days  10  and  12  indicates  that  the  IGWs  have 
a  period  with  respect  to  the  ground  of  about  12-24  h 
and  propagate  eastward  at  a  slower  rate  than  the  speed 
of  the  (westerly)  jet.  Following  Dunkerton  ( 1984),  we 
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Fig.  5.  Latitude -height  cross  section  of  horizontal  velocity  divergence  at  two  longitudes  following  the  main 
IGW  packet  on  consecutive  days:  (a)  day  10,  20°E;  (b)  day  1 1,  20°W.  Contour  interval  is  10“5  s‘L  Windspeed 
superposed  with  shading  increment  10  m  s_l  as  in  Fig.  1. 
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adopt  the  convention  that  the  wavevector  points  op¬ 
posite  the  mean  flow,  reflecting  the  fact  that  phase 
propagation  is  westward  with  respect  to  parcels  in  the 
jet.  Thus  the  intrinsic  frequency  is  positive  by  conven¬ 
tion,  and  the  frequency  relative  to  the  ground  is  mostly 
in  the  range  -2fl7day  *=  u>  *£  -7r/day.  For  a  zonally 
oriented  wavevector  with  zonal  wavenumber  k  and  a 
zonal  flow  u,  the  intrinsic  frequency  lu  -  w  -  ku .  Tak¬ 
ing  k  «  — 27t/(600  km),  w  =  — 27r/(  14  h),  and  u 
«  30  m  s~\  representative  of  the  zonal  wind  through 
the  wavepacket  at  1 30  mb  on  day  1 1 ,  yields  u>  «  ]  .9 
X  10-4  s-1.  Thus,  u>  > /(« 1.2  X  10"4  s-1  at  55°N). 
The  waves  propagate  westward  relative  to  the  wind 
with  an  intrinsic  period  of  about  9  h,  which  is  shorter 
than  the  inertial  period  of  14.6  h  at  this  latitude. 

This  quick  estimate  of  intrinsic  frequency  may  be 
compared  to  that  derived  from  the  dispersion  relation 
for  plane  hydrostatic  IGWs  of  the  form  =  <t>0 
exp/ ( k  •  x  -  ut);  namely, 


N2(k2  +  l2) 
U  ~f  mH-  l/(4//2) 


(3.1) 


with  /  and  m  being  the  meridional  and  vertical  wave- 
numbers,  respectively,  and  H  is  the  density  scale 
height.  (The  second  term  in  the  denominator  is  negli¬ 
gible  compared  to  the  first.)  On  day  1 1  the  wavepacket 
at  130  mb  is  centered  near  55°N,  and  the  horizontal 
mean  wind  and  wavevector  are  both  directed  zonally. 
The  horizontal  wavelength  is  approximately  600  km 
with  k2t>  /2,  and,  from  Fig.  6c,  the  vertical  wavelength 
is  about  4  km.  Substitution  into  the  dispersion  relation 
yields  w  «  1.7  x  10"4s_I,  which  agrees  reasonably 
well  with  the  other  estimate.  The  two  values  of  wave 
frequency  at  1 30  mb  indicate  that  gravity  wave  param¬ 
eters  are  consistent  with  an  IGW  interpretation  pro¬ 
vided  that  estimates  are  made  near  the  center  of  the 
wavepacket.  Exact  agreement  is  not  expected,  since  the 
waves  are  not  plane  and  the  dispersion  relation  ignores 
latitudinal  and  vertical  shear.  A  more  detailed  analysis 
of  to  is  presented  below. 

The  group  velocity  of  the  IGW  packet  may  be  de¬ 
rived  from  the  dispersion  relation: 


c,=jqM,<£±i!>i+5.  (3.2} 

m~io  [  ~m 


For  wave  parameters  near  the  130-mb  jet  axis  (where 

1*1  *>  Ul): 

Cg  =  (-7.1  x  102,  -0,4.7)  km  day"1  +  u.  (3.3) 

The  vertical  group  velocity  is  positive  since  m  is  neg¬ 
ative.  Upward  group  velocity  is  also  indicated  by  a 
wind  hodograph  taken  at  day  1 1  through  the  center  of 
the  wavepacket  at  53°N,  23°W  (Fig.  7).  The  hodograph 
of  perturbation  winds  (a  five-point  running  mean  hav¬ 
ing  been  removed)  shows  anticyclonic  (clockwise)  ro¬ 
tation  of  wind  with  increasing  height,  implying  that  the 


group  velocity  is  directed  upward.  The  major  axis  of 
the  hodograph  ellipse,  which  is  parallel  to  the  direction 
of  wave  propagation,  is  oriented  from  slightly  south  of 
east  to  north  of  west.  The  perturbation  horizontal  winds 
associated  with  the  IGWs  have  magnitudes  of  about  (5, 
3)  m  s”1  for  components  parallel  and  perpendicular  to 
the  direction  of  wave  propagation.  At  other  locations 
we  find  similarly  that  the  major  axis  of  the  hodograph 
ellipse  is  oriented  in  the  direction  of  wave  propagation 
suggested  by  Fig.  3  (not  shown). 

The  horizontal  group  velocity  is  directed  upstream 
relative  to  the  wind  and  has  a  magnitude  of  about  -8 
ms-1.  The  wavepacket  is  seen  to  propagate  eastward 
at  about  15-25  m  s“\  accelerating  to  the  higher  value 
after  day  1 1 .  Although  this  agrees  with  predictions,  the 
observed  distribution  of  the  wavepacket  does  not  nec¬ 
essarily  indicate  the  group  velocity  since  it  is  affected 
by  mechanisms  of  excitation  and  absorption,  as  well  as 
propagation. 

Upward  radiating  IGWs  such  as  these  will  encounter 
a  “Jones  critical  level”  when  |£|  |/|.  Using  wave 

parameters  typical  of  the  center  of  the  packet,  oriented 
in  the  direction  of  the  mean  flow,  we  estimate  uCl 
~  25  m  s-1  as  the  critical  mean  wind  speed.  In  this 
simulation  IGWs  are  seen  to  penetrate  slightly  higher 
than  the  Jones  critical  level,  to  levels  where  the  west¬ 
erly  windspeed  has  decreased  into  the  range  u  ~  10- 
25  ms-1.  This  penetration  above  the  Jones  critical 
level  is  likely  due  to  transience  in  the  forcing  of  the 
IGWs,  the  derived  uCL  being  appropriate  for  steady- 
state  conditions.  IGW  breakdown  is  thought  to  occur 
via  Kelvin -Helmholtz  (KH)  instability  due  to  vertical 
shear  of  the  wave’s  transverse  wind  component.  Our 
large-scale  model  cannot  resolve  the  KH  instabilities, 
of  course,  so  the  IGWs  are  absorbed  by  hyperdiffusion 
as  they  approach  the  critical  level.  This  could  happen 
either  through  contraction  of  the  vertical  wavelength  or 
rotation  of  the  horizontal  wavevector  away  from  the 
direction  of  the  mean  flow,  both  of  which  would  reduce 
the  intrinsic  frequency  and  vertical  group  velocity.  We 
believe  in  this  instance  (as  explained  further  in  section 
3d)  that  diffusion  of  model  IGWs  is  due  primarily  to 
contraction  of  horizontal  wavelength  and  rotation  of 
the  horizontal  wavevector  away  from  the  direction  of 
the  mean  flow.  Comparing  a  T84  version  of  this  sim¬ 
ulation  ( with  Az  =  700  m)  to  an  identical  T84  run  with 
double  vertical  resolution  (A z  =  350  m)  gave  strik¬ 
ingly  similar  IGW  and  vertical  wavelengths.  Propa¬ 
gation  of  IGWs  in  the  model  is  therefore  not  limited  by 
vertical  resolution .  Furthermore,  by  repeating  the  sim¬ 
ulation  at  higher  horizontal  resolution  while  keeping 
the  hyperdiffusion  coefficient  unchanged  gave  almost 
identical  IGWs.  Thus,  the  simulated  IGWs  are  well  re¬ 
solved  and  are  not  significantly  affected  by  inconsistent 
vertical  or  horizontal  resolution  as  discussed  by  Lind- 
zen  and  Fox-Rabinovitz  (1989). 

Along  the  flanks  of  the  jet  the  horizontal  wavevector 
is  rotated  inward  in  response  to  differential  advection 
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by  the  jet.  As  seen  in  Fig.  3,  the  horizontal  wavelength 
along  the  flanks  of  the  Jetstream  has  contracted  some¬ 
what  relative  to  that  along  the  jet  axis.  This  is  consistent 
with  waves  refracting  into  the  jetstream  as  discussed 
by  Dunkerton  (1984);  for  example, 

dl 

—  *  -kuy,  (3.4) 

etc.  On  the  southern  flank  of  the  jet  /  >  0  so  cgv  >  0, 
while  /  <  0  and  cgy  <  0  on  the  northern  flank. 

A  consequence  of  rotation  and  contraction  is  that 
wave  activity  along  the  jet  flanks  is  rapidly  dissipated 
by  the  model’s  hyperdiffusion.  This  mainly  accounts 
for  the  sensitivity  of  simulated  IGWs  to  the  model’s 
horizontal  resolution:  at  lower  resolution  (T84),  IGW 
activity  away  from  the  jet  axis  is  lost  even  more  readily 
to  diffusion.  We  expect  that  with  finer  horizontal  res¬ 
olution  well  beyond  T126,  the  trailing  IGWs  along  the 
flanks  of  the  jet  would  be  resolved  for  a  longer  time, 
evolving  to  even  shorter  horizontal  wavelengths.  Some 
evidence  in  support  of  this  claim  is  provided  by  results 
of  a  T156  run  as  discussed  in  section  3d. 

That  the  waves  are  being  absorbed  along  the  lateral 
boundaries  of  the  wavepacket  is  also  suggested  by  cal¬ 
culations  of  intrinsic  frequency,  determined  objectively 
from  the  spatial  distribution  and  temporal  behavior  of 
V-  \H.  Wavevector  orientations  and  magnitudes  were 
calculated  at  approximately  1°  intervals  within  the 
wavepacket  (at  130  mb)  by  applying  a  high-pass  spa¬ 
tial  filter  to  V  •  \H  ( to  eliminate  large-scale  ageostrophic 
circulations)  and  evaluating  the  horizontal  gradient  of 
this  quantity  at  the  zero  crossing.  Although  the  phase 
function  is  ill  behaved  in  a  complicated  wave  field  such 
as  this,  a  local  wavevector  direction  may  be  determined 
uniquely  as  parallel  ( or  antiparallel )  to  V(  V  •  vw)  at  the 
zero  crossing.  Half-wavelengths  were  estimated  as  the 
distance  between  adjacent  zero  crossings  in  the  direc¬ 
tion  of  the  wavevector.  To  get  reasonably  smooth  re¬ 
sults,  individual  estimates  within  a  4  x  4  matrix  of  grid 
points  were  binned  together  and  interpolated  back  to 
the  original  grid.  A  similar  binning  procedure  was  used 
for  frequency  w,  crudely  estimating  the  local  half-wave 
period  as  the  time  interval  between  adjacent  zero  cross¬ 
ings  in  a  bandpass-filtered  time  series  of  V  •  vw  at  each 
grid  point.  The  wavevector  direction  and  sign  of  u  were 
carefully  determined  by  examining  the  phase  propa¬ 
gation  around  the  time  of  interest.  Wavevectors  were 
aligned  in  the  direction  of  intrinsic  phase  propagation, 
that  is,  with  a  component  antiparallel  to  the  mean  flow. 
The  sign  of  w  was  negative  in  all  but  a  few  places;  this 
quantity  was  typically  smaller  than  the  advective  term 
k’U  by  a  factor  of  3—4.  To  evaluate  the  intrinsic  fre¬ 
quency  u>  =  u>  —  k*u,  a  small  amount  of  spatial 
smoothing  was  applied  to  the  horizontal  wind,  since 
IGWs  are  not  part  of  the  “mean  flow.’’ 

Figures  8a-c  show  V-  \H  ( with  high-pass  spatial  fil¬ 
ter)  at  130  mb  on  days  9.5,  1 1,  and  12.5,  illustrating 


the  early,  middle,  and  late  stages  of  IGW  evolution  at 
this  level,  respectively.  In  the  early  stage  (Fig.  8a)  there 
is  a  single  wavepacket  in  the  jet  stream  exit  region,  with 
most  wavevectors  pointing  backward  and  to  the  right 
with  respect  to  the  local  mean  flow  (for  an  observer 
facing  downstream).  By  day  1 1  (Fig.  8b),  the  packet 
has  begun  to  split,  part  of  it  going  south  into  the  cutoff 
low.  Wavevectors  in  both  halves  converge  into  the  jet, 
and  a  few  waves  manage  to  sneak  across  the  saddle  of 
the  cutoff  low  into  a  region  they  do  not  belong.  By  day 
12.5  (Fig.  8c),  the  main  part  has  advected  rapidly  east¬ 
ward  to  a  point  immediately  northwest  of  the  cutoff 
low  (reentering  from  the  left  side  of  the  diagram), 
while  the  wavepacket  trapped  within  the  cutoff  low  re¬ 
mains  largely  stationary.  Waves  crossing  the  saddle 
rapidly  disappear,  replaced  by  a  third  packet  ( unrelated 
to  the  others),  on  the  east  side  of  the  cutoff  low,  on  a 
collision  course  with  the  first. 

The  intrinsic  frequency  on  day  9.5  (Fig.  8d)  shows 
most  values  lying  in  the  range  f-4f,  with  a  flat  ledge 
of  smaller  values  (slightly  below /)  along  the  north¬ 
eastern  boundary  of  the  wavepacket.  The  variation  of 
w  is  due  mainly  to  the  angle  of  wavevector  with  respect 
to  the  mean  flow  and  not  the  wavelength  (which  is 
rather  uniform  on  this  day).  Although  penetration  be¬ 
yond  w  =  /  would  require  that  a  steady-state  wave  be 
evanescent,  in  this  case  the  phenomenon  does  not  per¬ 
sist  (as  noted  above,  in  connection  with  the  vertical 
propagation )  and  is  probably  due  to  transients  excited 
by  geostrophic  adjustment.  On  day  1 1  (Fig.  8e)  the 
region  of  subinertial  frequency  has  disappeared,  al¬ 
though  by  this  time  there  is  another  such  region  lying 
within  the  saddle  of  the  cutoff  low,  which  is  likewise 
temporary  and  quickly  damped.  Most  values  of  Cj  are 
again  in  the  range/- 3/.  By  day  12.5,  a  larger  per¬ 
centage  of  wavevectors  are  becoming  nearly  orthogo¬ 
nal  to  the  mean  flow,  and  values  of  w  are  somewhat 
smaller,  mostly /-2/(  Fig.  8f). 

It  must  be  noted  that  the  absorbing  action  of  the  side- 
walls  of  the  jet  is  an  artifact  of  the  model,  being  sen¬ 
sitive  to  the  model’s  horizontal  resolution  and  hyper¬ 
diffusion.  These  results  cannot  address  the  (still  un¬ 
answered  )  issue  of  whether  the  surface  u  =  /  acts  as 
a  critical  layer  to  absorb  IGWs  (as  in  vertical  shear)  or 
as  a  turning  point  to  reflect  them  (as  in  tidal  theory). 
This  matter  was  discussed  by  Kitchen  and  McIntyre 
( 1980)  among  others  but  deserves  further  study.  There 
is  little,  if  any,  evidence  for  reflection  of  model  IGWs, 
although  such  an  effect  might  be  masked  by  the  rapid 
evolution  of  the  flow.  It  was  suggested  by  Dunkerton 
(1984)  that,  due  to  horizontal  and  temporal  variations 
of  the  basic  state,  any  natural  waveguides  for  IGWs  are 
likely  to  “leak”  wave  activity  rather  than  contain  it.  In 
the  simulations  some  wave  activity  manages  to  escape 
to  the  stratosphere,  but  scale-dependent  damping  is  the 
dominant  effect. 

Radiation  to  the  middle  stratosphere  occurs  only  in 
the  trailing  part  of  the  main  wavepacket,  where  IGWs 
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Fig.  6.  Longitude- height  cross  section  of  horizontal  velocity  divergence  at  55°  latitude  on  consecutive  days: 
(a)  day  9;  (b)  day  10;  (c)  day  1 1;  (d)  day  12.  Contour  interval  is  I0~5  s"1.  Windspeed  superposed  as  in  Fig.  5 
with  contour  interval  10  m  s'1. 
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Fig.  6.  ( Continued) 
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Hodogroph  ot  23  W/53  N.  Time  =  11.0  days. 


Fig.  7.  Hodograph  of  the  vertical  profile  of  perturbation  horizontal 
wind  components  on  day  1 1  at  53°N,  23°W.  Perturbation  denotes 
deviation  from  winds  smoothed  with  a  five-point  running  mean  ver¬ 
tically.  Altitude  of  model  levels  indicated  in  km. 


tilt  northward  into  the  polar  night  jet  (not  shown). 
Propagation  at  this  location  is  possible  on  account  of 
increasing  mean  windspeed  despite  the  wavevector  be¬ 
coming  nearly  orthogonal  to  the  mean  flow.  Another 
packet  approaches  20  km  on  the  west  side  of  the  cutoff 
low  (not  shown).  These  results  demonstrate  that  ver¬ 
tical  propagation  of  model  IGWs  is  limited  to  the  re¬ 
gion  of  strong  mean  windspeed  where  tropospheric  and 
stratospheric  flows  have  approximately  the  same  direc¬ 
tion.  The  baroclinic  wave  distorts  the  tropospheric  flow 
from  a  zonal  orientation  but  is  itself  evanescent  in 
height  (i.e.,  does  not  similarly  distort  the  middle  strato¬ 
sphere).  Only  a  few  tropospheric  IGWs  can  therefore 
propagate  into  the  midstratospheric  circumpolar  vortex. 

d.  Effects  of  horizontal  and  vertical  resolution 

Model  IGWs  are  unaffected  by  vertical  resolution  if 
the  grid  spacing  is  sufficiently  fine  to  describe  the  ro¬ 
tation  of  perturbation  velocity  with  height,  as  shown  in 
Fig.  7.  Smaller  A z  did  not  affect  the  evolution  of  model 
IGWs,  indicating  that  the  IGWs  are  primarily  sensitive 
to  horizontal  hyperdiffusion. 

Integrations  were  performed  at  horizontal  resolution 
T42,  T84,  T126,  and  T156  with  hyperdiffusivity  co¬ 
efficients  z/6  =  1.37  x  1027,  1.02  x  1026,  1.36  x  1025, 
2.52  X  1 024  m6  s  ,  respectively,  such  that  the  damping 
time  at  the  smallest  resolved  scale  was  nearly  the  same 
in  each  run  (—20-30  min).  These  integrations  illus¬ 
trate  quite  dramatically  the  sensitivity  of  model  IGWs, 
as  shown  in  Figs.  9a~d.  Baroclinic  lifecycles  in  each 


run  were  at  virtually  the  same  stage  at  day  1 0,  but  major 
differences  are  seen  in  model  IGWs.  Contour  intervals 
are  the  same  in  each  panel,  indicating  a  large  increase 
of  IGW  amplitude  between  T42  (when  IGWs  were 
largely  absent)  and  T126.  Figure  9  also  demonstrates 
that  certain  parts  of  the  divergence  pattern  are  ade¬ 
quately  resolved  at  lower  resolution,  while  others  are 
not.  Increasing  horizontal  resolution  beyond  T84  did 
not  significantly  affect  the  horizontal  scale  of  IGWs  at 
the  rear  of  the  wavepacket  nor  the  structure  of  forced 
ageostrophic  circulations  evident  at  T42.  IGW  structure 
at  the  center  of  the  wavepacket  was  largely  unchanged 
going  from  T126  to  T156.  Differences  were  observed 
mainly  on  the  fringes  of  IGW  packets,  for  reasons  dis¬ 
cussed  in  the  previous  subsection. 

Interpretation  of  these  results  is  complicated  by  pos¬ 
sible  changes  of  the  baroclinic  wave  as  resolution  is 
increased,  for  example,  sharper  curvature  of  the  jet  and 
stronger  gradients  of  velocity.  An  example  of  wind- 
speed  at  day  10  for  T84  is  shown  in  Fig.  9e.  Wind- 
speeds  and  velocity  gradients  were  slightly  stronger  at 
higher  resolution,  by  a  few  percent.  The  basic  pattern 
of  horizontal  velocity  divergence  within  the  source  re¬ 
gion  (Fig.  9f)  was  similar,  except  for  substantially 
more  finestructure  at  higher  resolution  due  to  IGWs 
(not  shown).  Therefore,  changes  of  the  jetstream  pat¬ 
tern  were  minor  compared  to  the  observed  changes  of 
IGWs,  leading  us  to  believe  that  the  primary  effect  of 
increasing  resolution  in  this  range  (T42-T156)  is  to 
allow  the  region  of  imbalance  to  project  better  onto  the 
IGW  manifold,  irrespective  of  any  change  in  IGW 
sources.  We  cannot  rule  out  the  possibility  that  such 
sources  may  contract  further  in  horizontal  scale  or  be¬ 
come  more  intense  locally,  as  resolution  is  increased 
beyond  T156,  with  a  resulting  contraction  of  IGWs  and 
further  increase  of  amplitude.2  Indeed,  upper-level 
frontogenesis  is  expected  to  excite  mesoscale  gravity 
waves.  Model  experience  suggests,  however,  that  such 
waves  would  add  to  the  spectrum  of  IGWs  already  re¬ 
solved.  Important  to  this  argument  is  the  assumption 
that  model  IGWs  are  attributable  to  an  adequately  re¬ 
solved  region  of  imbalance  and  do  not  arise  entirely 
from  frontal  collapse  or  from  instabilities  lacking  a  pre¬ 
ferred  horizontal  scale. 

e.  Comparison  with  observations 

The  upper-level  IGWs  seen  in  this  model  study  are 
broadly  consistent  with  observations  of  IGWs  near 
midlatitude  jetstreams,  except  possibly  regarding  their 
direction  of  horizontal  energy  propagation.  Typically, 
the  lower  stratospheric  gravity  wave  field  is  dominated 
by  quasi-monochromatic  waves  with  horizontal  wave- 


:  Requirements  for  vertical  resolution  become  more  stringent  as 
horizontal  hyperdiffusion  is  decreased,  allowing  closer  approach  of 
IGWs  to  critical  levels. 
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lengths  of  several  hundred  kilometers  and  vertical 
wavelengths  of  1  -4  km  with  upward  energy  propaga¬ 
tion.  These  IGWs  appear  to  be  related  to  the  tropo¬ 
spheric  jet  stream,  though  it  has  proven  difficult  to  di¬ 
rectly  relate  them  to  a  specific  source  mechanism. 
Thomas  et  al.  (1992)  reported  a  case  study  based  on 
VHF  radar  observations  made  at  Aberystwyth,  Wales 
(52°N,  4°W).  They  found  southward  IGW  energy 
propagation  away  from  the  jet  stream,  which  was  pass¬ 
ing  to  the  north.  Prichard  and  Thomas  (1993)  exam¬ 
ined  stratospheric  IGWs  over  23  days  and  also  noted 
the  dominant  southward  IGW  propagation  direction.  In 
apparent  contrast  to  this  are  the  results  of  Sato  ( 1 994) , 
who  statistically  analyzed  IGWs  observed  by  radar  at 
Shigaraki,  Japan  (35°N,  136°E).  She  found  the  lower 
stratospheric  IGWs  to  be  most  energetic  in  wintertime 
and  to  propagate  southward  toward  the  strong  subtrop¬ 
ical  jet,  which  is  persistently  located  just  south  of  there 
during  winter.  Finally,  in  situ  observations  of  IGWs 
made  with  an  ER-2  aircraft  off  the  California  coast  in 
April  1986  found  large-amplitude  IGWs  in  a  jet  stream 
exit  region  with  propagation  roughly  parallel  to  the  jet 
axis  (X.  Pfister  19XX,  personal  communication).  At 
13-15-km  height  the  vertical  wavelength  of  these 
waves  was  1  -2  km,  and  the  horizontal  wavelength  was 
deduced  to  be  about  250  km. 

Thus,  while  lower  stratospheric  IGWs  primarily  tend 
to  propagate  meridionally,  there  is  some  uncertainty 
about  their  relation  to  the  jetstream.  There  may  be  large 
geographical  variability  in  the  relation  of  IGWs  to 
nearby  jetstreams.  Similarly,  it  is  likely  that  there  are 
several  different  types  of  jet  stream  developments  that 
result  in  significant  IGW  generation.  For  zonally  elon¬ 
gated  imbalances  to  the  jet  stream,  meridionally  radi¬ 
ating  (away  from  the  jet)  IGWs  should  result.  How¬ 
ever,  when  the  zonal  extent  of  the  imbalance  region  is 
reduced  relative  to  its  meridional  extent,  then  zonally 
propagating  IGWs  should  become  significant  (Luo  and 
Fritts  1993). 

Although  not  discussed  in  section  3c,  the  simulation 
at  a  later  time  shows  evidence  of  zonally  aligned  IGW 
phase  fronts  propagating  sideways  into  the  tropical 
lower  stratosphere.  The  waves  are  quite  weak  com¬ 
pared  to  those  near  the  jet,  but  their  (N-S)  wavevector 
orientation  is  consistent  with  IGW  observations  at  Are- 
cibo  Puerto  Rico,  described  by  Sato  and  Woodman 
( 1 982 )  and  later  authors.  This  would  suggest  that  not 
all  observed  IGWs  in  the  Tropics  are  due  to  underlying 
convection  but  that  some  of  them  may  originate  in  the 
midlatitude  jet. 

f.  Generation  of  IGWs  in  the  upper  troposphere 

Given  that  IGWs  are  difficult  to  observe  in  atmo¬ 
spheric  data  lacking  fine  vertical  and  temporal  resolu¬ 
tion,  it  would  be  desirable  to  relate  the  excitation  of 
IGWs  to  observable  parameters  in  the  large-scale  flow 
field.  The  following  thoughts  are  based  partly  on  earlier 


literature  and  partly  on  our  modeling  experience,  sug¬ 
gesting  how  this  might  be  done. 

As  discussed  in  section  1,  IGW  generation  is  ex¬ 
pected  during  periods  of  geostrophic  adjustment  of  the 
tropospheric  jet  stream  accompanying  synoptic-scale 
developments.  As  the  flow  evolves,  imbalances  de¬ 
velop  between  the  mass,  temperature,  and  velocity 
fields  causing  ageostrophic  motions  that  cannot  be  sup¬ 
ported  in  a  balanced  system  (quasi-geostrophic,  semi- 
geostrophic,  nonlinear  balance,  etc.).  In  this  event  the 
balance  breaks  down,  and  IGWs  are  radiated  away 
from  the  region  of  imbalance.  There  is  no  completely 
reliable  way  to  diagnose  the  conditions  required  for 
IGW  generation  however,  largely  for  lack  of  knowl¬ 
edge  of  the  balanced  state  to  which  the  flow  is  relaxing 
( Koch  and  Dorian  1988). 

The  horizontal  wind  divergence  field  near  the  level 
of  maximum  windspeed  ( -250  mb)  shows  the  signa¬ 
ture  of  slowly  varying  ageostrophic  circulations  that 
accompany  the  jet  stream  (Fig.  9f).  These  patterns  are 
similar  to  those  analyzed  by  Cammas  and  Raymond 
( 1 989 )  for  jet  streams  where  flow  curvature  was  strong, 
for  example,  their  Fig.  6.  In  such  cases  the  familiar 
quadrupole  pattern  of  upper-level  divergence  associ¬ 
ated  with  the  entrance  and  exit  regions  of  a  localized 
zonal  jet  stream  are  obscured  by  the  divergence  asso¬ 
ciated  with  the  flow  curvature.  The  appearance  of  a 
strong  divergence  center  in  the  jet  exit  region  at  250 
mb,  just  upstream  of  the  ridge,  implies  that  air  parcels 
rapidly  advected  through  this  region  will  experience  a 
large  local  rate  of  change  of  divergence,  suggestive  of 
flow  imbalance.  In  analogy  to  wave  generation  by  flow 
over  an  obstacle,  it  is  to  be  expected  that  the  dominant 
waves  generated  by  flow  through  this  jet  exit  region 
should  propagate  with  a  wavevector  approximately 
parallel  to  the  jet  axis,  as  observed  in  the  model  near 
the  center  of  the  jet.  Also  consistent  with  this  idea,  we 
find  that  the  IGW  pattern  is  approximately  phase 
locked  to  the  jet  exit  region  over  a  period  of  several 
days  (Fig.  6).  Similarly,  the  IGW  structure  in  the  jet 
exit  region  is  fairly  insensitive  to  improvements  in  res¬ 
olution. 

Typical  measures  of  imbalance  rely  on  estimating 
the  magnitude  of  those  terms  that  are  assumed  to  be 
negligible  in  any  particular  balanced-flow  state.  Thus, 
one  might  have  a  small  Rossby  number  (R„  =  U/fL 
<  1 )  for  geostrophic  balance,  a  small  Lagrangian 
Rossby  number  (RiL>  =  \d\/dt\/f  |V|  <  l)forsemi- 
geostrophic  balance  (since  it  implies  dY/dt  —  dYf/ 
dt),  and  a  small  material  derivative  of  horizontal  di¬ 
vergence  (dV-\H/dt  =€  10-9s“2)  for  the  nonlinear 
balance  equation. 

In  this  study  the  plots  of  dV  •  \H/dt  tended  to  mimic 
those  of  V  •  \H  ( e.g..  Fig.  3 ),  which  we  use  to  show  the 
IGWs.  This  quantity  is  therefore  an  accurate  but  re¬ 
dundant  estimate  of  imbalance.  Instead,  one  would  like 
a  measure  of  flow  imbalance  that  does  not  explicitly 
involve  V-  v„  and  may  be  diagnosed  from  the  observ- 
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able  flow.  The  Lagrangian  Rossby  number  appears 
suitable  for  this  purpose. 

When  applied  to  synoptic  data,  Koch  and  Dorian 
suggest  an  indicator  of  unbalanced  flow  is  the  Lagran¬ 
gian  Rossby  number  R(VL)  ^  0,5  in  a  region  where  |  V  | 
>  lOms**1  (to  avoid  spuriously  large  values  of  R$,L) 
from  the  denominator  being  small).  In  practice,  when 
calculating  d\/dt  from  synoptic  data,  the  local  wind 
tendency  dV/dt  must  often  be  neglected  relative  to  the 
advective  tendency  V*VV  because  of  the  crude  time 
resolution  of  the  observations.  For  this  simulation  we 
calculated  both  terms  of  dS Idt  and  found  that,  indeed, 
the  local  wind  tendency  is  usually  small  compared  with 
the  advective  tendency  V-  VV,  although  locally  it  can 
amount  to  ~30%  of  the  latter. 

Figure  10  shows  R{0L)  at  195  mb  on  day  11  of  the 
baroclinic  lifecycle  with  horizontal  velocity  divergence 
overplotted  at  130  mb,  showing  the  main  IGW  packet 
at  50°-60°N  and  a  second  wavepacket  west  of  the  cut¬ 
off  low  near  30°N.  Nearer  to  the  jet  core,  R},L)  was 
calculated  at  a  slightly  lower  level,  since  the  wave 
source  lies  below  where  the  waves  are  prominently 
seen.  In  this  simulation  large  values  of  R„L)  are  found 
extending  down  to  the  level  of  the  jet  maximum,  near 
250  mb.  In  Fig.  10,  these  large  values  (heavy  lines)  are 
seen  to  coincide  with  the  location  of  IGW  activity.  Two 
regions  of  large  R{0L)  occur:  one  between  45°  and  60°N 
that  overlies  the  main  IGW  packet  in  the  jetstream  exit 
region  and  a  second  region  lying  southwest  of  the  cut¬ 
off  low,  suggesting  that  two  separate  regions  of  imbal¬ 
ance  exist.  To  be  sure,  large  R{aL)  suggests  where  an 
imbalance  exists  but  says  nothing  about  the  intensity 
of  resulting  IGWs.  Since  IGWs  radiate  unbalanced  en¬ 
ergy  away  from  the  jet,  their  intensity  should  be  ap¬ 
proximately  proportional  to  |  V  | 2 .  This  is  consistent 
with  results  of  mechanistic  numerical  experiments  by 
Van  Tuyl  and  Young  (1982),  who  found  that  the 
strength  of  unbalanced  motions  (IGWs)  increases 
roughly  as  the  square  of  jet  core  strength.  Thus, 
stronger  IGW  radiation  is  to  be  expected  from  regions 
of  imbalance  (large  R{0L))  where  windspeeds  are 
strong.  This  may  explain  why  the  IGWs  are  much 
stronger  in  the  high-latitude  region  of  imbalance  than 
at  low  latitudes  where  windspeeds  are  weaker. 

A  combination  of  Lagrangian  Rossby  number  and 
mean  flow  kinetic  energy,  then,  provides  a  useful  di¬ 
agnostic  of  IGW  excitation  in  the  model.  It  remains  to 
be  seen  whether  it  is  systematically  useful  in  observa¬ 
tions. 

4.  Discussion 

The  ubiquitous  nature  of  IGWs  in  the  atmosphere 
has  long  been  recognized,  yet  their  effects  are  only  be¬ 
ginning  to  be  considered  in  the  context  of  the  general 
circulation.  This  situation  is  partly  due  to  their  sub¬ 
synoptic  scale  and  preferred  occurrence  at  upper  tro¬ 
pospheric  levels  and  above.  IGWs  are  nevertheless 


Fig.  10.  Lagrangian  Rossby  number  at  195  mb  on  day  1! 
(heavy  lines),  with  horizontal  velocity  divergence  at  130  mb  super¬ 
posed:  RlJ]  is  plotted  only  where  the  windspeed  exceeds  10  m  s“\ 
the  dotted  line  denoting  this  boundary. 


likely  to  play  an  important  role  in  several  aspects  of 
atmospheric  circulation. 

Our  simulations  show  that  IGWs  are  generated  in  the 
upper  troposphere,  propagate  energy  upward,  and  are 
absorbed  approaching  the  Jones  critical  level  (either  in 
weak  winds  or  as  the  wavevector  rotates  away  from  the 
mean  flow).  As  discussed  in  section  3c,  IGW  break¬ 
down  is  thought  to  occur  via  KH  instabilities  (unre¬ 
solved  in  the  model),  which  develop  as  amplitude  in¬ 
creases  and  vertical  wavelength  decreases.  These  insta¬ 
bilities  are  expected  in  the  region  of  strong  transverse 
vertical  shear,  at  which  point  the  flow  is  statically  stable 
(Dunkerton  1984;  Fritts  and  Rastogi  1985).  IGW 
breaking  is,  therefore,  a  potentially  efficient  means  of 
cross-isentropic  transport  and  may  be  responsible  for 
most  of  the  vertical  diffusion  found  in  the  lower  strato¬ 
sphere,  a  region  where  systematic  cross-isentropic 
transport  due  to  diabatic  heating  is  quite  weak.  Vertical 
mixing  due  to  IGW  breaking  must  vary  geographically 
and  seasonally,  just  as  IGW  generation  does. 

In  addition  to  cross-isentropic  transport  induced  by 
IGWs;  there  is  an  isentropic  component  due  to  advec- 
tion  by  the  horizontal  velocity.  One  normally  thinks  of 
gravity  wave  motions  as  fast,  reversible  undulations  su¬ 
perposed  on  the  “slow  manifold”  of  rotational  motions 
(e.g.,  nonlinear  baroclinic  instabilities  and  breaking 
planetary  Rossby  waves),  the  latter  being  mainly  re¬ 
sponsible  for  isentropic  mixing  on  a  global  scale.  Parcel 
displacements  due  to  IGWs,  nevertheless,  can  in  certain 
cases  be  large,  due  to  their  large  horizontal  scale  and 
low  intrinsic  frequency — particularly  as  the  critical 
level  is  approached  and  amplitude  increases.  In  this 
case  an  air  parcel  advected  through  the  IGWs  will  un¬ 
dergo  a  considerable  horizontal  excursion  during  a  half 
wave  period  before  the  transverse  velocity  reverses 
sign. 


/ams  v3628  1203  Mp  18  Tuesday  Jul  18  1 1:42  AM  LP:  J  At  Sci  (August  15  95)  1203 


76 


Month  1995 


O’SULLIVAN  AND  DUNKERTON 


0019 


Under  such  circumstances,  any  background  tracer 
field  with  a  nonzero  gradient  orthogonal  to  the  IGW 
wavevector  will  experience  considerable  deformation. 
This  is  true,  in  particular,  for  potential  vorticity  (PV) 
near  the  tropopause.  Such  deformations  of  PV  on  an 
isentropic  surface  are  in  fact  seen  in  these  simulations. 
Figure  1 1  shows  PV  on  the  430-K  surface,  illustrating 
large  meanders  of  the  (approximately  conserved)  PV 
contours  across  the  jet  in  a  region  of  strong  IGWs.  The 
meanders  are  rapidly  strained  to  smaller  scales  by  the 
background  flow  and  lost  because  of  limited  model  res¬ 
olution.  It  is  certain  that  they  will  lead  to  irreversible 
isentropic  transport  and  that  this  will  be  especially 
strong  in  the  cross-jet  direction  if  real  IGWs  behave 
like  those  simulated  here.  In  ongoing  work  we  are  in¬ 
vestigating  the  effect  of  IGWs  on  isentropic  tracer 
transport  in  the  lower  stratosphere  by  comparing  trans¬ 
port  and  mixing  in  simulations  with  and  without  IGWs. 
This  is  relevant  in  light  of  recent  contour  advection 
studies  (Plumb  et  al.  1994;  Waugh  et  al.  1994;  Waugh 
and  Plumb  1994;  Norton  1994)  in  which  the  isentropic 
tracer  redistribution  is  computed  from  advecting  winds 
known  only  at  low  resolution  and  as  infrequently  as 
daily.  Recognizing  the  importance  of  unbalanced  mo¬ 
tions,  Pierce  and  Fairlie  (1993)  speculate  on  the  role 
of  IGWs  in  enhancing  cross-jet  tracer  dispersion,  per¬ 
haps  enhancing  the  permeability  of  the  lower  strato¬ 
spheric  vortex  at  its  edge. 

Figure  12  shows  a  latitude-height  cross  section 
through  the  IGW-induced  meanders  shown  in  Fig.  1 1 . 
This  illustrates  the  effect  of  quasi-horizontal  isentropic 
IGW  motion  on  the  PV  field,  which  is  representative 
of  a  conserved,  vertically  stratified  tracer.  An  important 
implication  of  such  isentropic  transport  in  general  is 
that  IGWs  may  directly  cause  stratosphere -tropo¬ 
sphere  exchange.  If  the  IGW  breaking  level  occurred  a 


few  kilometers  lower  than  seen  in  Fig.  12,  the  quasi¬ 
horizontal  motions  would  involve  the  tropopause 
[commonly  defined  by  the  2-PVU  (potential  vorticity 
units)  contour].  For  instance,  if  the  IGW  motions  of 
Fig.  1 1  were  to  occur  100  K  lower  at  approximately 
330  K  where  the  isentropic  surface  intersected  a  slop¬ 
ing  tropopause,  then  the  2-PVU  contour  would  be  in¬ 
volved  in  lateral  meanders.  The  corresponding  version 
of  Fig.  1 1  would  then  show  filaments  of  tropospheric 
air  being  advected  into  the  higher  PV  stratosphere  and 
filaments  of  high  PV  stratospheric  air  similarly  being 
advected  into  the  troposphere.  The  alternating  filaments 
of  stratospheric  and  tropospheric  air  appear  like  inter¬ 
leaved  fingers  on  an  isentropic  surface.  Much  of  this 
exchange  would  be  irreversible,  given  the  rapid  strain¬ 
ing  of  features  by  the  background  wind  field.  There  is 
observational  evidence  for  stratosphere- troposphere 
exchange  by  waves  of  short  vertical  wavelength  that 
produce  interleaving  laminae  of  tropospheric  and 
stratospheric  air  at  and  above  the  tropospheric  jet  (Sha¬ 
piro  et  al.  1980).  Danielsen  et  al.  (1991 )  and  Kritz  et 
al.  ( 1991 )  show  evidence  of  such  filaments  mixing  air 
along  isentropes  in  the  stratosphere  above  the  jet 
stream. 

Aside  from  tracer  transport,  IGWs  are  important  in 
that  they  constitute  a  subsynoptic  signal  having  large 
amplitude  in  the  vicinity  of  jet  streaks  and  jet  stream 
exit  regions.  Van  Tuyl  and  Young  (1982)  emphasized 
this  aspect  of  IGWs  and  the  importance  of  not  sup¬ 
pressing  them  during  the  process  of  assimilating  new 
data  into  numerical  forecast  models.  IGWs  can  also 
play  a  role  in  organizing  convection  and  precipitation. 
Koch  and  Dorian  (1988)  observed  IGWs  from  a  jet 
streak  exit  region,  organizing  a  sequence  of  thunder¬ 
storms  as  they  propagated  through  a  region  of  weak 
conditional  convective  stability. 


Fig.  11.  Isentropic  potential  vorticity  on  the  430-K  surface  on  day  11.  Contour  interval  is  1  PVU 
(=  10~6  nr  s' 1  K  kg-1).  This  isentropic  surface  is  near  90  mb  at  the  location  of  the  IGW  packet. 
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Fig.  12.  Latitude-height  cross  section  of  PV  at  30°W  intersecting  the  lateral 
meanders  of  Fig.  11.  Contour  interval  is  1  PVU. 


Lastly,  the  upward  momentum  flux  due  to  IGWs  is 
considered.  Given  the  widespread  nature  of  IGWs  and 
their  preferred  location  in  jetstreams,  it  might  be 
thought  that  their  upward  flux  of  momentum  could 
be  important.  However,  the  integrated  momentum 
flux  is  proportional  to  energy  times  a  factor  6lJ2 
=  Vl  -  (//u>)2,  showing  that  smaller  scale  high- 
frequency  gravity  waves  are  more  efficient  in  fluxing 
momentum  vertically  than  IGWs  (Fritts  and  van  Zandt 
1993).  (The  factor  6+  in  their  formula  is  extraneous  to 
the  momentum  flux  and  may  be  safely  omitted.)  It  is 
of  interest  to  note  that  the  Eliassen-Palm  flux  ap¬ 
proaches  zero  more  rapidly  as  i  ->/,  by  an  additional 
factor  of  <5_ .  At  the  center  of  the  wavepacket  on  day 
1 1,  u  'w'  «  -0.08  m2  s'2;  a  time-  or  area-average  flux 
would  be  even  smaller.  Both  are  small  numbers  com¬ 
pared  to  the  momentum  flux  near  the  tropopause  esti¬ 
mated  by  Fritts  and  van  Zandt. 

IGWs  could,  nevertheless,  have  some  impact  on  mo¬ 
mentum  transport  by  modulation  of  higher-frequency 
gravity  waves  (Broutman  and  Young  1986).  Such 
waves  might  be  absorbed  at  IGW-induced  critical  lev¬ 
els  or,  conversely,  escape  to  higher  altitudes.  Partial 
reflection  due  to  steep  gradients  created  by  IGW  mix¬ 
ing  is  another  possibility.  Neither  would  be  very  im¬ 
portant  in  the  simulation  where  IGWs  account  for  a 
small  fraction  of  the  total  wind. 


5.  Conclusions 

The  excitation  and  propagation  of  inertia-gravity 
waves  (IGWs)  generated  by  an  unstable  baroclinic 
wave  was  examined  with  a  high-resolution  3D  nonlin¬ 
ear  numerical  model.  IGWs  arose  spontaneously  as  the 
tropospheric  jet  stream  was  distorted  by  baroclinic  in¬ 
stability  and  strong  parcel  accelerations  took  place,  pri¬ 
marily  in  the  jet  stream  exit  region  of  the  upper  tro¬ 
posphere.  Subsequent  propagation  of  IGWs  occurred 
in  regions  of  strong  windspeed,  in  the  tropospheric  and 
stratospheric  jets,  and  in  a  cutoff  low  formed  during  the 
baroclinic  lifecycle.  IGWs  on  the  flanks  of  these  jets 
were  rotated  inward  by  differential  advection  and  sub¬ 
sequently  absorbed  by  the  model’s  hyperdiffusion.  This 
was  the  principal  mechanism  of  IGW  absorption  in  the 
model,  the  results  being  insensitive  to  vertical  resolu¬ 
tion.  IGWs  in  reality  are  thought  to  break  via  Kelvin- 
Helmholtz  (KH)  instability,  a  process  unresolved  in 
the  3D  model.  Although  the  model’s  absorption  of 
IGWs  at  the  sidewalls  of  the  jet  is  artificial,  IGW  prop¬ 
agation  was  for  the  most  part  confined  to  regions  of 
Cj  >  f.  Only  a  few  IGWs  were  able  to  penetrate  the 
middle  stratosphere,  due  to  weak  winds  or  an  unfavor¬ 
able  alignment  of  wavevector  with  respect  to  the  mean 
flow  (except  beneath  the  polar  night  jet). 

IGWs  are  important  both  as  a  synoptic  signal  in  the 
jet  stream  that  may  influence  subsequent  tropospheric 
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developments  and  as  a  source  of  isentropic  or  cross- 
isentropic  mixing  in  the  lower  stratosphere.  Our  results 
demonstrated  for  the  first  time  numerically  a  significant 
isentropic  displacement  of  potential  vorticity  isopleths 
due  to  IGWs  above  the  tropopause.  Since  conditions 
for  IGW  propagation  are  favorable  within  a  jet,  a  region 
of  strong  isentropic  PV  gradient,  it  is  likely  that  inertia- 
gravity  waves  affect  the  permeability  of  the  lower 
stratospheric  vortex  and  may  in  some  instances  lead 
to  stratosphere-troposphere  exchange,  particularly  if 
IGW-induced  parcel  displacements  are  large  and  irre¬ 
versible.  Although  IGW  breaking  may  help  ensure  that 
the  mixing  is  permanent,  it  is  interesting  to  note  from 
the  simulations  that  straining  by  the  large-scale  flow, 
without  IGW  breaking,  was  sufficient  to  guarantee  ir¬ 
reversible  mixing  of  PV. 

Observations  of  tracer  laminae  in  the  lower  strato¬ 
sphere  sometimes  suggest  a  concurrent  IGW  signal, 
while  others  do  not;  it  is  possible  that  “fossil”  laminae 
were  due  to  an  earlier  IGW  event  or,  perhaps,  the  baro- 
clinic  instability  itself.  The  large-scale  shear  distends 
tracer  isopleths  into  many  thin  filaments  (e.g.,  Plumb 
et  al.  1994).  It  will  be  interesting  to  compare  mixing 
rates  with  and  without  IGWs  to  assess  their  role  in  is¬ 
entropic  transport. 

Considerably  higher  vertical  (and  horizontal)  reso¬ 
lution  is  required  to  determine  the  mechanisms  of  IGW 
breaking  and  the  role  of  IGWs  in  cross-isentropic  trans¬ 
port.  This  issue  is  important  in  the  lower  stratosphere 
where  other  mechanisms  of  cross-isentropic  transport 
are  quite  weak.  The  chemistry  of  this  region  is  sensitive 
to  the  rate  at  which  constituents  are  diffused  vertically 
and  exchanged  between  the  troposphere  and  strato¬ 
sphere. 

IGW-induced  mixing  is  expected  to  vary  as  a  func¬ 
tion  of  location  and  season.  The  climatology  of  IGWs 
and  their  relation  to  tropospheric  sources  are  essentially 
unknown.  If  progress  is  to  made,  it  will  begin  with  an 
analysis  of  IGWs  in  relation  to  the  synoptic-scale  flow 
at  locations  where  continuous  radar  or  lidar  observa¬ 
tions  of  high-frequency  motions  are  routinely  available. 
Relationships  established  locally  may  then,  as  a  first 
approximation,  be  extrapolated  to  global  models. 
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